| kriging [message #74488] |
Wed, 19 January 2011 08:00 |
ggggg
Messages: 2
Registered: January 2011
|
Junior Member |
|
|
hallo,
i have some precipation measurments from various swiss stations (over
500 points). now i try to interpolate this data with kriging. but i
don't know how to do that.
as an example i tried do do it with 5 points with the values in array
f, and the long-and latitudes in the arrays lat and lon. (see further
down)
i plotted the map of switzerland with the subprogram "maps" , but then
i don't know what to do. i found the griddata function but how i have
to use it?
hope you can help me!
thanks
ms
code:
pro grid
lat=[45.9,47.1,46.5,46,47]
lon=[6,9,8,10,7.2]
f=[100,200,11,120,100]
maps
stop
end
Pro maps
maplimit = [45.7, 5.9, 47.8, 10.6] ;NMM22 europe
; Handle TrueColor displays:
DEVICE, DECOMPOSED=0
; Load discrete color table:
tek_color
; Match color indices to colors we want to use:
black=0 & white=1 & red=2
green=3 & dk_blue=4 & lt_blue=5
; Set up an orthographic projection centered over the north
; Atlantic.Fill the hemisphere with dark blue. Specify black
; gridlines:
;MAP_SET, /ORTHO, 40, -30, 23, /ISOTROPIC, $
; /HORIZON, E_HORIZON={FILL:1, COLOR:dk_blue}, $
; /GRID, COLOR=black
MAP_SET, /ORTHOGRAPHIC, 48, 7, limit = maplimit,/iso, $
/GRID, LONDEL=15, LATDEL=15,color=0, title ='stationen',/hires
; Fill the continent boundaries with solid white:
MAP_CONTINENTS, /FILL_CONTINENTS, COLOR=black,/hires
; Overplot coastline data:
MAP_CONTINENTS, /COASTS, COLOR=white, /hires
; Add rivers, in light blue:
;MAP_CONTINENTS, /RIVERS, COLOR=lt_blue
; Show national borders:
MAP_CONTINENTS, /COUNTRIES, COLOR=white, MLINETHICK=2, /hires
end
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
