We have: ETOPO5 and Predicted Bathy (Smith and Sandwell) To use our code you should make: ~jhc/local/sgi/bin in your path. on solomon we have : /drives/solomon/disk2/Irish_ACQ/ where there are two files: -rw-r--r-- 1 jhc omg 37325824 Nov 3 15:58 world_etopo.r4 -rw-r--r-- 1 jhc omg 273716224 Oct 10 13:10 world_predicted.r4 make copies of these files in your own space..... These are the two data sets stored as floating point arrays in our image format. You should get to know this format. It is described in : ~jhc/code/subr/array.h You can make a blank mapsheet using make_blank: prompt>, make_blank Warning : no outfile name input the blank array file prefix please... NZ_region_predicted do you want an 8bit array (1) or floating point array (2) 2 input Map Sheet Characteristics: ------ which Ellipsoid ? ---- WGS84 : (1) ---- NAD27 : (2) ---- CLK66 : (3) ---- AGD84 : (4) ---- SPHER : (5) 1 what type of projection ? ---- Mercator (orthogonal) : (1) ---- UTM (orthogonal) : (2) ---- UTM (rotational) : (4) ---- Mercator (rotational) : (5) 1 input L/L of top LH corner of area (deg min deg min) (-VE in southern and western hemispheres!) -20 0 160 0 -20 0.000000 --- 160 0.000000 input L/L of bottom RH corner of area (deg min deg min) (-VE in southern and western hemispheres!) -40 0 180 0 -40 0.000000 --- 180 0.000000 input projection latitude: -30.0 xlen 1928110.139266 , ylen 2234386.893581 input desired metres per pixel 2000 mpp 2000.000000 array size 964 1117 prompt>, mv NZ_region_predicted.blank NZ_region_predicted.r4 prompt>, cp NZ_region_predicted.r4 NZ_region_etopo.r4 A better set of Geobounds for Bay of Plenty is ~~~~` -30 0 -170 0 -40 0 -180 0 choose your own bounds .... ----------------------- now you have two blank flaoting pointing arrays. Now you need to fill them with ETOPO or prediscted bathy:. >, patchWORLD -etopo -into NZ_region_etopo.r4 -source world_etopo.r4 >, patchWORLD -predict -into NZ_region_predicted.r4 -source world_predicted.r4 This takes a while, let it plug away.... You can look at what you made by typing: >, jview NZ_region_predicted.r4 and use the left mouse button to pull iout a zoomed subarea. Arrow kys allow you to panm around. Click with the middle mouse button to get topo profiles from click to click. Noe try sun-iluminating the data: >, addSUN NZ_region_etopo.r4 -azimuth 315 NZ_region_etopo.sun_315 >, jview NZ_region_etopo.sun_315 Make an 8bit depth encoded image of the data: >, r4to8bit -low -3000 -high 1000 NZ_region_etopo.r4 NZ_region_etopo.8bit >, jview NZ_region_etopo.8bit Can make a gif of the image by: >, stretchacres -low 0 -high 255 -pgm -in NZ_region_etopo.sun_315 -out NZ_region_etopo.sun_315.pgm >, xv NZ_region_etopo.sun_315.pgm to convert from a PGM to gif or tif or jpg or whatever... another way to make a color picture is: >, mix_ci -c NZ_region_etopo.8bit -i NZ_region_etopo.sun_315 -ppm -m NZ_region_etopo.ppm >, xv NZ_region_etopo.ppm to convert from a PPM to whatever... These pgms/ppms wil lbe what you draw into Coreldraw to have an the underlying image for your maps.... Now contour the data: prompt >, contgrid NZ_region_predicted.r4 NZ_region_predicted.cont_500m float NZ_region_predicted.r4, plot NZ_region_predicted.cont_500m low 0 high 0 OMG1 OMG1 array header data ---- 964 1117 - 160.000000 180.000000 -40.000000 -20.000000 - 2000.000000 2000.000000 floating point data type -- undefined rgb_bw -- Mercator Projection -- WGS84 Ellipsoid merid -- -30.000000- sun - 0.000000 0.000000 0.000000 - view - 0.000000 0.000000 0.000000 max 0.000000 min 0.000000 read in depth ranges: -5000.000000 -1000.000000 input minlevel,maxlevel,interval -8000 1000 500 number of levels: 19 frequency of contour labels (eg.100(m)) : 1000 0 1117 labeled : -1000.000000 labeled : -1000.000000 ....... take a look at the contours: jview NZ_region_predicted.sun_315 -vectbin NZ_region_predicted.cont_500m now make a vector plot of the bounds of the mapsheet: >, lineplot -area NZ_region_predicted.r4 -grid 600 100 enter projection meridian -30.0 enter scale 1000000 and look at this : >, viewmeta plot.meta convert to hpgl for importing into Coreldraw (to act as the base for the image map sheets). >, hpgl plot.meta basemap_border.plt if you want plot the contours: >, lineplot -area NZ_region_predicted.r4 -grid 600 100 -vectbin NZ_region_predicted.cont_500m 1 1 etc..... Sometimes you get too many contour vectors. It is often worth averaging down the grid before contouring: >, decr4 -count 4 -in NZ_region_predicted.r4 -out temp.dec.r4 >, contgrid temp.dec.r4 smaller.cont_500m etc.... All this should keep you occupied for a while... Source.


Яндекс.Метрика Рейтинг@Mail.ru Free Web Counter
page counter
Last Modified: January 31, 2014 @ 12:00 am