Wednesday, September 15, 2010

station lat/lon to 3 decimal places

On the IRIS site the seismic station locations are only given to 2 decimal places lat/lon, which is +/- l km. They don't want to real the precise position of the instruments.

If you need more info, location to 3 decimal places is available in the .sac files.

(from rob)

In the sac header files you can find the latitude and longitude to at least 3 decimal places. You can access this either via sac:
read (file)
lh stla stlo

or a couple matlab options based on which sac reader you're using
with load_sac:
[header, data] = load_sac('file');
header.stla, header.stlo

(In the above two cases stla is short for station latitude and stlo is short for station longitude)

or with readsac
[time, amp, header] = readsac('file');
header(17)
header(18)

in this case header(17) is the latitude and header(18) is longitude.

No comments:

Post a Comment