2 Problemes: Export-Map-Format; "Overlay Maps"

[andreas.hoppe@gsf.de]

Hello GSTAT-Users,

first I want to excuse me for the following faults in the english grammer and 
spelling (I hope there aren't too much).

For some months I am working on a GIS-project, forthat we bought "ArcView 
3.0a" (unfortunately without "Spatial Data Analyst"-Modul). 
In learning more about GIS and kriging I found GSTAT 2.0f (an GNUPLOT 3.6/347 
and  PCRaster 2.0) on the Internet (It's good!), but I'm an engineer with low 
geographical knowledge.

O.k., here my task:

We have some measuring sites (only a few - about 20) were the precipitation 
(rain) is analysed. Our idea to get a more detailed map of the deposition of 
the chemicals was to use a more narrow net of precipitation amount sites 
(over 150) and  multiply the (kriged map of the) concentration of the 
chemicals (of the few sites) with the (kriged map of the) raining amounts (of 
the many sites).

I think the "right" method is "Inverse Distance weighted Interploation", if 
not: Please tell me!

O.k. - so far so good.


1) How do I multiply one pixel_x1_y1 of map1 with pixel_x2_y2 of map2 and 
receive finally a map3 ?
   (My Input-map-format is taken from the PCRaster distribution; '*.map')
   (Both maps are equal in size and locations)


2) By the way: All attemps to generate a different output-map in GSTAT with 
"set mapformat=(-4 to +3)" failed.
   The most interessting map format for me would be ArcInfo/binary (because 
of the use of ArcView).
   The Error message was:

         NULL argument in function 'write_csf()'
         this is a bug, please send the next line to e.pebesma@geog.uu.nl
         mapio.c, line 1092, 2.0f (July 1998) dos/dpmi

Any ideas ??

3) As I mentioned I use the free PCRaster (professional) package (because 
nothing other worked) for displaying my maps on the screen. This works fine: 
display map1.map

But: All attemps failed here, too, to get a print of it.

How can I print '*.map' maps? No matter if with PCRaster or with GNUPLOT or 
something else ? 
Or how can I tell GSTAT to generate another map that is printable?


To complete this mail with all information I can give, here is my comand file 
of GSTAT:
(modelled on example 8)

#
# ordinary block kriging on a mask
#
set gnuplot='gpt36';
set gpterm='svga';
set mapformat=0; 
#    0=input,          1=PCR,            2=ArcInfo,      3=Idrisi,
#   -2=ArcInfobinary, -3=Idrisi binary, -4=Gnuplot binary
data(rain): 'rain_ges.txt', x=1, y=2, v=13, 
max=12, radius=5000; # force;    # within a distance of 5000
# set plotfile='rain.gp';
mask: 'c:\programm\pcrast2\test2.map'; 
blocks;
blocksize: dx=40, dy=40;        # 40 x 40 block averages
# set output = 'rain_ok.out';   # output ascii file
predictions(rain): 'c:\programm\pcrast2\rain_bk2.map';




All hints and tips to solve my problemes are welcome.
Please e-mail to: andreas.hoppe@gsf.de
Thanks a lot.


*************************************************
**********         Andreas Hoppe       **********
********        GSF-Forschungszentrum     *******
******          Tel. +49/89/3187-4618      ******
****            Fax. +49/89/3187-3313        ****
**           e-mail: andreas.hoppe@gsf.de      **       
*************************************************