[R-sig-Geo] [STATSGRASS] Re: Unable to read GRASS6 data using gstat 2.4.5
Thomas Adams
Thomas.Adams at noaa.gov
Wed Aug 24 17:05:42 CEST 2005
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
More information about the R-sig-Geo
mailing list