Oliver,
Would it just be a matter of counting the vertices of each polygon,
selecting the polygons with 3 vertices, and comparing the value of the
first and last vertex as you are building wkt value string below?
Doug
>The problem with your rings isn't that they are not closed, it is that
Derek,
You can do this in GRASS ( using gdal in the background) by doing
the following:
1) create a workspace in the same projection as your existing raster (
when you first start GRASS it will ask you to set up a workspace with
projection information)
2) import the existing raster la
Looks good! Thanks for the link!
Doug
Doug Newcomb
USFWS
Raleigh, NC
919-856-4520 ext. 14 doug_newc...@fws.gov
-
The opinions I express are my own and are not representative of the
official
>Hi,
>I did mean that with "contains vector graphics" and "zoom in and enjoy".
>-Jukka-
Does that also mean extract vector layers from a geopdf and save out as
shapefile ? For those of us trying to translate government-distributed
data into something useable
Doug
Doug Newcomb
USFWS
For what it's worth,
What do you consider large and what kind of method for downscaling you
wish to use?
I do this in GRASS using the r.neighbors command,
http://grass.fbk.eu/grass64/manuals/html64_user/r.neighbors.html. I've
done an averaging method to convert a 20 ft elevation grid ( st
Thanks! Sorry for the noise.
Doug
Doug Newcomb
USFWS
Raleigh, NC
919-856-4520 ext. 14 doug_newc...@fws.gov
-
The opinions I express are my own and are not representative of the
official poli
Hi Folks,
I had seen libjpeg-turbo ,
http://www.libjpeg-turbo.org/Main/HomePage, in the announcement for Fedora
14 bu had forgotten about it until Franks' FOSS4G talk last week. Has
anyone tried using it in GDAL?
Doug
Doug Newcomb
USFWS
Raleigh, NC
919-856-4520 ext. 14 doug_newc...@fw
Stephen,
You can script the import, classification, and display with python
in GRASS in a couple of ways. I generally just run python from the GRASS
command window to echo the GRASS commands and variables I need, but there
are other ways of doing it.
Doug
Doug Newcomb
USFWS
Raleigh,
Jorge,
> Hello,
>
> Using GDAL library is posible to read and write ArcGIS binary
GRID files?
I generally work with GRID files in GRASS. As mentioned before, gdal (
r.in.gdal in this case) imports the GRID layers fine. To export back to
ArcGIS, I use the r.out.arc command,
http://g
I probably do it the hard way. The USDA flies summer color imagery in
North Carolina every couple of years for crop verification. They will put
draft versions of the imagery up in single image per county mosaics in
MrSid format for download about 6-8 months before they make the final
quarter
>As we are running out of 1.x numbers (1.8.0 hopfully coming out soon),
I'm
>thinking about a 2.0 release for something like the fall of 2012 - just
in
>time for the world to end.
That should simplify keeping track of the EPSG codes...
Doug
Doug Newcomb
USFWS
Raleigh, NC
919-856-4520 ext. 1
Hi folks,
I'm working with someone trying to merge and reproject 1108 4 band tiff
images into 1 BIGTIFF image. The command line options are as follows:
gdalwarp -s_srs EPSG:26918 -t_srs EPSG:32119 -wo 'SOURCE_EXTRA=3' -wo
'SKIP_NOSOURCE=YES' -wm 4000 --config 'GDAL_CACHEMAX=2000' -tps -rcs
-
Jason,
If you're working with vector data, why not throw the data into
Postgresql/Postgis, http://postgis.refractions.net, and use the spatial
operators there to select/buffer/intersect the vector geometries as you
describe. http://postgis.refractions.net/documentation/manual-1.4
/ch07.html
Carolyn,
I may be running into something similar. Try setting the origin
projection of the source tiff using the -s_srs parameter and see what
result you get.
Doug
Doug Newcomb
USFWS
Raleigh, NC
919-856-4520 ext. 14 doug_newc...@fws.gov
---
>Folks,
>I would note that large -wm values can be very counter productive when
>used in combination with SKIP_NOSOURCE. The problem is that the larger
>the chunk size, you run into the chance that a large window will intersect
>a small amount of data and the whole window ends up being processed
>Hi Doug,
>I finally tried your parameters and they did work fine for me also. I
>had something like hundred geotiffs, 400 MB each, and I was pushing
>them to bigtiff mosaic. I tried first with your *.tif selection and then
>again by using a virtual raster file as source, created from Mapser
>> Shaun Kolomeitz wrote:
>> > could push beyond 1GB/s. Currently to process (mosaic) an 80GB image
it
>> > takes several days to complete. This is also on 32bit hardware, and I
>>
>> Shaun,
>
>> I suspect that there is a gross issue with how the warping is being
done,
>> and that it could be spe
Bruce,
If you're on a linux system with perl:
perl -pi -e 's/\t/,/g' filename
is a perl one-liner that should replace all of the tabs with commas in the
file, working on the file in place. I've used similar commands on text
files up to 379 GB in size, although it does take a while for fil
I'm confused, I thought spatialite support was being added to ogr in 1.7 .
http://www.gdal.org/ogr/drv_sqlite.html
As to ESRI, I asked at the last DOI meeting prior to the ESRI conference
this year about spatialite support and they said they had no plans to
support it. Of course, they used to
David,
zip format files for 32bit WindowsXP maxes out at 4GB. You can
create larger zip files on Vista. I'm wildly guessing that you are
running into the 4GB limit as well on reading your zip files with . 7zip,
http://www.7-zip.org/, might be able to handle zip format files larger
tha
Even,
Thanks! That did the trick!
Doug
Doug Newcomb
USFWS
Raleigh, NC
919-856-4520 ext. 14 doug_newc...@fws.gov
-
The opinions I express are my own and are not representative of the
off
Hi,
I'm trying to compile gdal 1.6.1 on a Centos 5.3 x86_64 computer. The make
portion works fine, but the make install bombs with :
/usr/bin/ld: skipping incompatible /usr/lib/libcom_err.so when searching
for -lcom_err
/usr/bin/ld: skipping incompatible /usr/lib/libcom_err.a when searching for
There seems to be an issue with the python scripts in fwtools 2.3.0 and
2.3.1
Example:
C:\Program Files\FWTools2.3.1>gdal_polygonize
Traceback (most recent call last):
File "C:\PROGRA~1\FWTOOL~1.1\bin\gdal_polygonize.py", line 34, in
import gdal, ogr, osr
File "C:\PROGRA~1\FWTOOL~1.1\p
Hi Folks,
Before I try to clumsily reinvent the wheel, I'd like to know if
someone has a method to :
1) use ogr2ogr to reproject a shape file ( easy enough ).
2) after the above projection copy the .xml file (if it exists)
from the old shapefile to the new shapefile (easy as w
24 matches
Mail list logo