Help: 2 problemes to solve...

[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      **       
*************************************************