[R-sig-Geo] [STATSGRASS] Re: Unable to read GRASS6 data using gstat 2.4.5
Markus Neteler
neteler at itc.it
Wed Aug 24 17:19:08 CEST 2005
[.,, maybe we should take
statsgrass at grass.itc.it, r-sig-geo at stat.math.ethz.ch
out of the loop as it is unrelated there.
]
http://grass.itc.it/pipermail/grass-commit/2005-May/017534.html
Check for nad2bin in configure script, so that if PROJ isn't installed
properly the user hears about it before compilation fails in lib/proj
What does
which nad2bin
say? Maybe 'nad2bin' is not in path.
nad2bin needs to be in path for installing datum gridshift files.
Markus
On Wed, Aug 24, 2005 at 11:05:42AM -0400, Thomas Adams wrote:
> Markus,
>
> I have a problem with ./configure with 6.1-CVS
> (grass-6.1.cvs_src_snapshot_2005_08_20.tar.gz) — I have not had previous
> problems configuring & compiling. I get a configure error:
>
> Unable to locate nad2bin.
>
> I've been compiling against proj-4.4.9 previously for GRASS 6.0.1, etc.
> just fine.
>
> Tom
>
>
> Markus Neteler wrote:
>
> >Tom,
> >
> >I have generated the location and imported the file:
> >
> >cat temp_20050801_08 | v.in.ascii -z out=temp_20050801_08 x=3 y=2 z=4
> >col='id varchar(3), y double, x double, elevation double, temperature
> >double' skip=1
> >
> ># file: gstat.cmd.temps
> >#
> ># gstat command file
> >#
> >data(temp_20050801_08): 'temp_20050801_08', x=1, y=2, v=1;
> >variogram(temp_20050801_08): 0.2772 Nug(0) + 0.85 Sph(987521);
> >
> >Running that:
> >
> >gstat gstat.inp
> >gstat: Linux version 2.4.5 (17 June 2005)
> >Copyright (C) 1992, 2004 Edzer J. Pebesma
> >using Marsaglia's random number generator
> >data(temp_20050801_08): WARNING: Vector used instead of sites.
> >WARNING: Illegal TimeStamp string
> >GRASS site list temp_20050801_08: 0 cat, 3 dim, 1 str, 4 dbl.
> >Warning: gstat/grass: specify z coordinate if you want to use it
> >
> >gstat/grass: 117 sites read successfully.
> >temp_20050801_08 (GRASS site list)
> >attribute: col[1] [x:] x_1 : [ -91, -78]
> >n: 117 [y:] y_2 : [ 35, 43]
> >sample mean: 39.8882 sample std.: 1.84532
> >[starting interactive mode]
> >
> >cat temp_20050801_08.est
> >#gstat Linux 2.4.5 (17 June 2005) [gstat gstat.inp]
> >#sample semivariogram
> >#Wed Aug 24 15:45:39 2005
> >#data(temp_20050801_08): 'temp_20050801_08', x=1, y=2, v=1;
> >#[1] mean: 39.8882 variance: 3.4052
> >#cutoff: 5.08806 interval width: 0.339204
> >#direction: total
> ># from to n_pairs av_dist semivariogram
> > 0 0.339204 30 0.247674 0.0128483
> >0.339204 0.678408 89 0.527931 0.0741331
> >0.678408 1.01761 142 0.847027 0.174927
> >1.01761 1.35682 175 1.18733 0.294193
> >[...]
> >
> >
> >No crash so far (I used 6.1 on Linux).
> >
> >Disussion of the warnings:
> >
> >data(temp_20050801_08): WARNING: Vector used instead of sites.
> >-> that's correct since sites are treated as vectors in GRASS 6
> >
> >WARNING: Illegal TimeStamp string
> > I have fixed that (bug in siteslib) in 6.1-CVS and
> > 6.0-CVS (will appear in a potential 6.0.2)
> > This or similar messages no longer appear.
> >
> >GRASS site list temp_20050801_08: 0 cat, 3 dim, 1 str, 4 dbl.
> >Warning: gstat/grass: specify z coordinate if you want to use it
> >-> depends on you, the user.
> >
> >I suggest to update to 6.1-CVS.
> >
> >A new run with GRASS' siteslib recompiled leads to:
> >
> >Copyright (C) 1992, 2004 Edzer J. Pebesma
> >using Marsaglia's random number generator
> >data(temp_20050801_08): WARNING: Vector used instead of sites.
> >GRASS site list temp_20050801_08: 0 cat, 3 dim, 1 str, 4 dbl.
> >Warning: gstat/grass: specify z coordinate if you want to use it
> >gstat/grass: 117 sites read successfully.
> >
> >which looks good (to me).
> >
> >
> >Markus
> >
> >
> >On Wed, Aug 24, 2005 at 09:37:14AM -0400, Thomas Adams wrote:
> >
> >
> >>Markus,
> >>
> >>Not the v.in.ascii command, but the columns are:
> >>
> >>ID|latitude|longitude|elevation|temperature
> >>
> >>elevation in feet MSL
> >>temperature in deg F
> >>
> >>Tom
> >>
> >>
> >>Markus Neteler wrote:
> >>
> >>
> >>
> >>>Tom,
> >>>
> >>>[location created]
> >>>before guessing, please send the v.in.ascii command to me
> >>>(or tell me what the columns represent)
> >>>
> >>>Markus
> >>>
> >>>On Tue, Aug 23, 2005 at 08:25:10AM -0400, Thomas Adams wrote:
> >>>
> >>>
> >>>
> >>>
> >>>>Markus,
> >>>>
> >>>>Attached are the files you requested. I have also attached the ascii
> >>>>vector file as well.
> >>>>
> >>>>Thanks so much for your help!
> >>>>
> >>>>Regards,
> >>>>Tom
> >>>>
> >>>>
> >>>>Markus Neteler wrote:
> >>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>>>>On Mon, Aug 22, 2005 at 04:19:16PM -0400, Thomas Adams wrote:
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>>List:
> >>>>>>
> >>>>>>I have point temperature data in GRASS6 (as vectors); I'm using GSTAT
> >>>>>>2.4.5 which I have compiled with GRASS support (and confirmed with
> >>>>>>gstat -v). I have:
> >>>>>>
> >>>>>># file: gstat.cmd.temps
> >>>>>>#
> >>>>>># gstat command file
> >>>>>>#
> >>>>>>data(temp_20050101_08): 'temp_20050101_08', x=1, y=2, v=1;
> >>>>>>variogram(temp_20050101_08): 0.2772 Nug(0) + 0.85 Sph(987521);
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>From the GRASS6 prompt, when I try to execute:
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>>gstat gstat.cmd.temps
> >>>>>>
> >>>>>>I get:
> >>>>>>
> >>>>>>data(temp_20050801_08): WARNING: Vector used instead of sites
> >>>>>>WARNING: Invalid absolute datetime format
> >>>>>>Memory fault
> >>>>>>
> >>>>>>I understand "WARNING: Vector used instead of sites", but I am not
> >>>>>>sure of the implication since in GRASS6 sites are vectors. The 2nd
> >>>>>>message followed by the Memory Fault is very puzzling.
> >>>>>>
> >>>>>>Regards,
> >>>>>>Tom
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>Tom,
> >>>>>
> >>>>>could you send the temp_20050101_08 files to me. Since it is
> >>>>>vector, I would need the mapset/vector/temp_20050101_08/ directory
> >>>>>and the output of
> >>>>>g.proj -wef > /tmp/temp_gstat.prj
> >>>>>to recreate the location here.
> >>>>>
> >>>>>Markus
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>--
> >>>>Thomas E Adams
> >>>>National Weather Service
> >>>>Ohio River Forecast Center
> >>>>1901 South State Route 134
> >>>>Wilmington, OH 45177
> >>>>
> >>>>EMAIL: thomas.adams at noaa.gov
> >>>>
> >>>>VOICE: 937-383-0528
> >>>>FAX: 937-383-0033
> >>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>>
> >>>
> >>>
> >>>
> >>>>PROJCS["Lambert Conformal
> >>>>Conic",GEOGCS["wgs84",DATUM["D_WGS_1984",SPHEROID["wgs84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["standard_parallel_1",60],PARAMETER["standard_parallel_2",30],PARAMETER["latitude_of_origin",39],PARAMETER["central_meridian",-86],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]
> >>>>
> >>>>
> >>>>
> >>>>
> >>>
> >>>
> >>>
> >>>>ALN|38.9|-90.05|545.|64.0
> >>>>ARR|41.77|-88.47|705.|63.0
> >>>>BLV|38.55|-89.95|460.|68.0
> >>>>BMI|40.48|-88.92|880.|69.0
> >>>>CMI|40.04|-88.27|750.|68.0
> >>>>CPS|38.57|-90.15|415.|67.0
> >>>>DEC|39.87|-88.92|670.|66.0
> >>>>DPA|41.91|-88.24|760.|66.0
> >>>>LWV|38.76|-87.61|430.|71.0
> >>>>MDH|37.78|-89.25|410.|64.0
> >>>>MDW|41.79|-87.75|615.|73.0
> >>>>MLI|41.45|-90.51|590.|66.0
> >>>>MTO|39.48|-88.28|720.|67.0
> >>>>MVN|38.32|-88.87|480.|66.0
> >>>>MWA|37.75|-89.01|470.|66.0
> >>>>ORD|41.98|-87.9|665.|72.0
> >>>>PIA|40.66|-89.69|660.|70.0
> >>>>PWK|42.11|-87.9|645.|70.0
> >>>>RFD|42.2|-89.09|735.|64.0
> >>>>SPI|39.84|-89.68|600.|68.0
> >>>>SQI|41.75|-89.67|645.|63.0
> >>>>UGN|42.42|-87.87|730.|70.0
> >>>>AID|40.1|-85.62|920.|64.0
> >>>>BAK|39.27|-85.9|655.|70.0
> >>>>BMG|39.14|-86.61|840.|67.0
> >>>>EKM|41.72|-86.|780.|70.0
> >>>>EVV|37.97|-87.57|384.|71.0
> >>>>EYE|39.83|-86.3|820.|70.0
> >>>>FWA|41.|-85.2|791.|66.0
> >>>>GSH|41.52|-85.79|830.|72.0
> >>>>GYY|41.62|-87.42|590.|70.0
> >>>>HNB|38.25|-86.95|530.|68.0
> >>>>HUF|39.45|-87.3|500.|67.0
> >>>>IND|39.73|-86.28|792.|73.0
> >>>>LAF|40.42|-86.9|600.|68.0
> >>>>MIE|40.24|-85.4|935.|67.0
> >>>>SBN|41.7|-86.32|773.|68.0
> >>>>VPZ|41.45|-87.01|770.|67.0
> >>>>BWG|37.|-86.43|535.|70.0
> >>>>FFT|38.18|-84.9|800.|71.0
> >>>>JKL|37.6|-83.32|1365.|70.0
> >>>>LEX|38.03|-84.6|966.|70.0
> >>>>LNP|36.98|-82.53|2680.|66.0
> >>>>LOU|38.23|-85.66|545.|72.0
> >>>>LOZ|37.08|-84.08|1211.|68.0
> >>>>PAH|37.03|-88.58|325.|65.0
> >>>>SDF|38.18|-85.73|477.|74.0
> >>>>ARB|42.22|-83.75|835.|67.0
> >>>>AZO|42.23|-85.55|870.|71.0
> >>>>BEH|42.13|-86.42|645.|70.0
> >>>>BTL|42.3|-85.25|950.|71.0
> >>>>DET|42.42|-83.02|625.|74.0
> >>>>DTW|42.23|-83.33|640.|71.0
> >>>>IRS|41.81|-85.44|925.|64.0
> >>>>JXN|42.27|-84.47|998.|69.0
> >>>>OEB|41.93|-85.05|960.|66.0
> >>>>RMY|42.25|-84.96|940.|70.0
> >>>>TTF|41.94|-83.43|615.|68.0
> >>>>YIP|42.24|-83.53|715.|71.0
> >>>>DKK|42.49|-79.27|690.|71.0
> >>>>JHW|42.15|-79.27|1720.|63.0
> >>>>AOH|40.71|-84.03|975.|67.0
> >>>>BJJ|40.87|-81.88|1135.|66.0
> >>>>BKL|41.52|-81.68|585.|73.0
> >>>>CAK|40.92|-81.43|1208.|69.0
> >>>>CGF|41.57|-81.48|880.|72.0
> >>>>CLE|41.4|-81.85|777.|69.0
> >>>>CMH|40.05|-82.88|812.|68.0
> >>>>CVG|39.07|-84.67|869.|74.0
> >>>>DAY|39.75|-84.23|1002.|68.0
> >>>>DFI|41.34|-84.43|705.|65.0
> >>>>FDY|41.02|-83.67|797.|65.0
> >>>>ILN|39.42|-83.82|1075.|67.0
> >>>>LCK|39.82|-82.93|745.|66.0
> >>>>LHQ|39.76|-82.66|865.|65.0
> >>>>LUK|39.1|-84.43|480.|68.0
> >>>>MFD|40.82|-82.52|1295.|68.0
> >>>>MGY|39.59|-84.22|960.|72.0
> >>>>MNN|40.62|-83.06|995.|66.0
> >>>>OSU|40.08|-83.07|900.|69.0
> >>>>PHD|40.47|-81.42|890.|66.0
> >>>>TDZ|41.56|-83.48|620.|66.0
> >>>>TOL|41.6|-83.82|669.|69.0
> >>>>YNG|41.27|-80.67|1178.|67.0
> >>>>ZZV|39.95|-81.9|881.|66.0
> >>>>AGC|40.35|-79.92|1250.|68.0
> >>>>AOO|40.3|-78.32|1469.|67.0
> >>>>BFD|41.8|-78.63|2142.|60.0
> >>>>BTP|40.78|-79.95|1250.|63.0
> >>>>DUJ|41.18|-78.9|1824.|66.0
> >>>>ERI|42.08|-80.17|732.|71.0
> >>>>GKJ|41.63|-80.21|1400.|68.0
> >>>>JST|40.32|-78.83|2285.|66.0
> >>>>LBE|40.28|-79.4|1050.|70.0
> >>>>PIT|40.5|-80.22|1137.|66.0
> >>>>BNA|36.17|-86.7|590.|72.0
> >>>>CSV|35.95|-85.08|1810.|70.0
> >>>>MQY|36.02|-86.52|545.|73.0
> >>>>TRI|36.48|-82.4|1525.|67.0
> >>>>TYS|35.82|-83.98|980.|70.0
> >>>>BCB|37.22|-80.42|2130.|64.0
> >>>>HLX|36.77|-80.82|2740.|64.0
> >>>>HSP|37.95|-79.82|3790.|63.0
> >>>>MKJ|36.9|-81.35|2560.|68.0
> >>>>PSK|37.13|-80.68|2110.|64.0
> >>>>ROA|37.32|-79.97|1149.|68.0
> >>>>SHD|38.27|-78.85|1200.|63.0
> >>>>VJI|36.68|-82.03|2070.|68.0
> >>>>BKW|37.78|-81.12|2504.|64.0
> >>>>BLF|37.3|-81.2|2870.|67.0
> >>>>CKB|39.29|-80.23|1200.|70.0
> >>>>CRW|38.35|-81.63|939.|65.0
> >>>>EKN|38.88|-79.85|1992.|63.0
> >>>>HLG|40.17|-80.64|1190.|66.0
> >>>>HTS|38.37|-82.55|827.|67.0
> >>>>LWB|37.87|-80.4|2303.|61.0
> >>>>MGW|39.65|-79.92|1240.|67.0
> >>>>PKB|39.27|-81.57|615.|67.0
> >>>>
> >>>>
> >>>>
> >>>>
> >>>
> >>>
> >>>
> >>>
> >>--
> >>Thomas E Adams
> >>National Weather Service
> >>Ohio River Forecast Center
> >>1901 South State Route 134
> >>Wilmington, OH 45177
> >>
> >>EMAIL: thomas.adams at noaa.gov
> >>
> >>VOICE: 937-383-0528
> >>FAX: 937-383-0033
> >>
> >>
> >>
> >
> >
> >
>
>
> --
> Thomas E Adams
> National Weather Service
> Ohio River Forecast Center
> 1901 South State Route 134
> Wilmington, OH 45177
>
> EMAIL: thomas.adams at noaa.gov
>
> VOICE: 937-383-0528
> FAX: 937-383-0033
>
--
Markus Neteler <neteler itc it> http://mpa.itc.it
ITC-irst - Centro per la Ricerca Scientifica e Tecnologica
MPBA - Predictive Models for Biol. & Environ. Data Analysis
Via Sommarive, 18 - 38050 Povo (Trento), Italy
More information about the R-sig-Geo
mailing list