On 06-09-2011 17:48, Frank Warmerdam wrote:
On Tue, Sep 6, 2011 at 8:48 AM, Joaquim Luis wrote:
OGR_G_* functions are intended to operate on OGRGeometry objects while
OSR* functions operate on OGRSpatialReferences. Your hSRS object is an
OGRSpatialReference and not suitable to pass to
Hi
The doc of OGR_G_GetCoordinateDimension states at
http://www.gdal.org/ogr/ogr__api_8h.html#8a5e78b0753339d1cdd282b0e151d28f
that
...
Returns:
in practice this always returns 2 indicating that coordinates
are specified within a two dimensional space.
but that's not true as I just
Hi,
Mirone has those same block algorithms (and more, like slope by fitting
a bilinear surface) and let change the the box size.
Joaquim
Evan,
Thanks for the info.
Jay
On Sat, Oct 8, 2011 at 12:29 PM, Even Rouault
mailto:even.roua...@mines-paris.org>>
wrote:
Le samedi 08 octobre 2
Hi,
I'm having troubles reading a SeaWifs L2 hdf file due to a bug in
the hdf driver. The problem here is that not all Subdatasets are
read.
Below is the list givem by gdalinfo, which lacks 2 subdatasets: the
cntl_pt_cols & cntl_pt_rows
To prove I attach a
DF4 driver will
open it pick it up.
Etienne
On Sat, Oct 22, 2011 at 9:57 AM, Joaquim Luis <jl...@ualg.pt>
wrote:
Hi,
I'm having troubles reading a
ou can
read lat/lon from GCPs. Or do you also feel that only ~100 GCPs are
not enough ?
I'm working with MODIS data and in the Arctic warping results look
ugly sometimes. I was thinking if I should try to get lat/lon grids
from an HDF file and increase number of GCPs.
Anton
On 10/23/201
On 24-10-2011 14:04, Jose Gomez-Dans wrote:
Hi Joaquim,
On 24 October 2011 13:56, Joaquim Luis <mailto:jl...@ualg.pt>> wrote:
Anton,
I don't remember the details because I programmed that some time
ago, but from what I recall that's the most accurate way of
SUBDATASET_20_DESC=[391x6] scan_ell (32-bit floating-point)
SUBDATASET_21_NAME=NETCDF:"S1998031140424.L2_MLAC_OC.x.hdf":nflag
SUBDATASET_21_DESC=[391x8] nflag (16-bit integer)
SUBDATASET_22_NAME=NETCDF:"S1998031140424.L2_MLAC_OC.x.hdf":tilt_ranges
SUBDATASET_22_DESC=[20x2] t
Hi,
I made a nice little tool for Mirone where one can very easily mosaic
either several flavors of SRTM grids or satellite images from Bing servers.
My question is about the SRTM30 grids (30 arc minutes) that are
available here (SRTM plus)
ftp://topex.ucsd.edu/pub/srtm30_plus/srtm30/data/
whi
Evan, Frank
Thanks very much. I keep insisting on the esri .hdr because I have a
function to do all the work but it's definitively time to write one to
do the same for .vrt (and learn more about it a same time).
Joaquim
Le lundi 14 novembre 2011 18:23:23, Joaquim Luis a écrit :
H
Derek,
These are the type of operations that I recommend GMT (built with GDAL
off course). With GMT5 ii as simple as
grd2xyz yourGeotiff > outfile
this saves x,y,z triplets. Add -ZTLa to write only the z values (many
other options available)
Joaquim
Hello,
I am new to GDAL and raster p
On 27-11-2011 00:07, Martin Chapman wrote:
All,
I am using a Windows 7 x64 machine with VS 2010 with both x86 and x64
compile tools. When I try to build gdal for x64 using the make file
from the Win64 command prompt and WIN64 = YES in the nmake.opt file
everything compiles fine. (using gda
Hi,
It's been a while (well, a couple of months) since I built on OSX but I
only did a svn update and I now get these linking errors
stall_name /usr/local/lib/libgdal.1.dylib -compatibility_version 17
-current_version 17.0 -Wl,-single_module
ld: warning: path '/System/Library/Frameworks/
On 14-12-2011 07:01, Even Rouault wrote:
Le mercredi 14 décembre 2011 01:48:36, Joaquim Luis a écrit :
Hi,
It's been a while (well, a couple of months) since I built on OSX but I
only did a svn update and I now get these linking errors
stall_name /usr/local/lib/libgdal.1.
Sorry for the noise that I introduced by mixing the GDAL and GMT's
lists. I already replied to Andreas only because I though my message
went only to him, but since it didn't here is a copy of it.
Joaquim
"
Andreas,
I'm very sorry but I sent that message to you by mistake as it was meant
as a
float pafLocalMin = afWin[0];
float pafLocalMax = afWin[0];
for ( int k = 1; k < 9; k++)
{
if (afWin[k] < pafLocalMin)
{
pafLocalMin=afWin[k];
}
if (afWin[k] > pafLocalMax)
{
pafLocalMax=afWin[k];
}
On 16-12-2011 03:31, Stephen Mather wrote:
Hmm. I don't follow? The revision would be this:
float GDALSlopeHydroAlg (float* afWin, float fDstNoDataValue, void* pData)
{
// Hydrologic Slope is the max
// local slope btw center cell and adjacent cells
const double radiansToDegrees
On 16-12-2011 15:32, Jan Hartmann wrote:
Hi all,
Is it possible to do a first order GCP-warp *without* shearing? I have
a rectangular scan with very imprecise control points, and I want it
georeferenced without any distortion. With a first order gdalwarp I
get back a map in the form of a para
On 16-12-2011 18:50, Jan Hartmann wrote:
On 16-12-2011 18:25, Joaquim Luis wrote:
Jan,
What you want is to apply a transformation called (at least on the
Matlab parlance) "Linear conformal" which an afine transform but
without the the shear terms. Mirone has it implemented and if y
On 21-12-2011 20:43, Etienne Tourigny wrote:
Yves,
processing of netcdf files like this may be much easier using NCO
(http://nco.sourceforge.net/)
example:
ncap -s "WND=sqrt(u^2+v^2)" in.nc out.nc
Or with GMT
grdmath U 2 POW V 2 POW ADD SQRT = out.nc
this, with GDAL's help, had the further
Etienne,
For what it worth, I repeated the tests as for Windows in
http://trac.osgeo.org/gdal/ticket/3166 and got the expected behavior.
Joaquim
Folks,
I would appreciate minimal testing for netcdf-4 in osx before
releasing rc2 (which should take a few days at the most).
Wiliam:
The drive
Hi,
I made my first attempt to build the python bindings on Windows and
like others have complained before ... nothing happened
I did set PYDIR and PY_INST_DIR but at the end no errors nor anything
else python related.
Next I went to swig\python and did
python.exe setup.py build
which resul
Generally, to build Python extensions on Widnows, you need to use the same
version of MSVC that was used to compile your Python binary. For python 2.6, I
believe that the provided binaries are compiled with MSVC 2008 for example.
See http://bugs.python.org/issue5235
If you have not that version
On 11-01-2012 08:11, Timur Chabayev wrote:
Hi,
I'd like to use GDAL in my program to open raster images (MFC), but I have one problem - when I try to open images with one raster band, these image are drawn gray (or, in few cases, almost black). I' ve looked for a solution, but haven't found any.
Hi,
I'm puzzled with this behavior of the vsi* functions
The first call that uses vsigzip because file is gzipped
/vsigzip/vsicurl/http:.../...30S.ACE2.gz gives an error. However, if I
remove the the '/vsigzip' than the command runs fine and I'm able to
access the file.
Is this intended (it
On 21-01-2012 08:35, Even Rouault wrote:
Le samedi 21 janvier 2012 03:21:43, Joaquim Luis a écrit :
Hi,
I'm puzzled with this behavior of the vsi* functions
The first call that uses vsigzip because file is gzipped
/vsigzip/vsicurl/http:.../...30S.ACE2.gz gives an error. However, if I
r
On 21-01-2012 21:15, Even Rouault wrote:
Even,
Thanks for the explanation + fix but I was not aware that the correct
syntax is to use a double slash. The docs does not sate that and in fact
the examples there http://trac.osgeo.org/gdal/wiki/UserDocs/ReadInZip
show (for example)
ogrinfo -ro -al
Also, I forgot ask before. The ACE2 driver says that elevation data is
return as float but what is the 'NoDataValue' for this?
I don't know if there's an explicit nodata, but you can perhaps derive it from
the source/quality/confidence layer datasets ?
For example :
gdalinfo
/vsigzip//vsicurl
Hi,
A couple of months ago I posted a message here about subdatsets not
being read by the HDF4 driver and opened this ticked.
http://trac.osgeo.org/gdal/ticket/4305
Jose Gomez-Dans reminded me that this was a even older issue and pointed
me to
http://osgeo-org.1560.n6.nabble.com/gdal-dev-s
On 21-02-2012 14:11, Rutger wrote:
Hey,
TJMartin wrote
...
Is there a way to combine the hillshade and topo map together - ie add a
transparency or opacity to the hillshade.
...
A common way is to convert your (i assume) RGB topo map to HSV color space
and replace the intensity band (v) with
5) Looked at GMT and I cannot find the options to load two rasters on
top of each other.
Tim,
You don't load the two grids on top of one-another. Illumination is NOT
like transparency. It uses the gradients of the grid, projects it into
the direction of light and than change the color in t
On 22-02-2012 00:54, Jay L. wrote:
List,
I am attempting to reproject 8 gtiffs, each in a local equirectangular
projection, back to gcs for use in geoserver. I immediately went to
gdalwarp to perform this task and am having an issue. The input files
are stored in positive longitude, while t
On 10-04-2012 13:14, Robert Zermeno wrote:
Ladies and Gents,
I have yet to properly configure GDAL to compile with a repository
OpenJpeg V2, but I wanted to see how GDAL uses openjpeg tile request.
I created my own project and modified IReadBlock() to be standalone to
view a tile request.
Hi,
I made some false reports in the past about the build failures that
ended up being caused by no starting by a "make clean", but that is not
the case this time.
I'm getting these errors on linking the current head.
Joaquim
LIBCMT.lib(dosmap.obj) : error LNK2005: _errno already defined in
On 12-04-2012 13:07, Even Rouault wrote:
Selon Joaquim Luis:
Hi,
I made some false reports in the past about the build failures that
ended up being caused by no starting by a "make clean", but that is not
the case this time.
I'm getting these errors on linking the current he
I made some false reports in the past about the build failures that
ended up being caused by no starting by a "make clean", but that is not
the case this time.
.^^
Clean and rebuild GDAL.
As I first said ...
___
gdal-dev m
On 12-04-2012 16:24, Jeff McKenna wrote:
Hi Joaquim,
I've definitely been in a similar situation before. Most of the time,
when I get those 'LIBCMT' or 'MSVCRT' conflict messages, it is because I
built a dependent library with the /MT option by mistake, causing the
GDAL built to try to include
On 12-04-2012 16:00, Ivan Lucena wrote:
Olá Joaquim,
Sometimes a SVN update bring things that require a complete clean up before we
can rebuild it again, so you need to run:
nmake -f makefile.vc clean
nmake -f makefile.vc
nmake -f makefile.vc install
Ivan, obrigado but ... I wish it was that
On 12-04-2012 22:32, Frank Warmerdam wrote:
On Thu, Apr 12, 2012 at 2:31 PM, Frank Warmerdam wrote:
On Thu, Apr 12, 2012 at 2:09 PM, Joaquim Luis wrote:
While at this, could I make a small request that VERSION variable in
nmake.opt be wrapped in a IFNDEF so that we can optionally set it
It also built fine with VS2010 (32 bits) but one need to put
continuation characters on OPTFLAG= ...
Joaquim
Folks,
I'm looking for a but of assistance from GDAL hackers
who use Visual Studio 2008 or later.
I'd be thankful if someone could conduct the following test for me
against the current
Hi again,
A couple of weeks ago I referred to a problem building the Win version
and at the end thought the problem was due to the WEBP driver. But now
I'm getting the same type of errors with another driver and, ... one nuance
To remember, the errors are are of the type
LIBCMT.lib(dosmap.ob
On 23-04-2012 17:28, Even Rouault wrote:
Selon Joaquim Luis:
Hi again,
A couple of weeks ago I referred to a problem building the Win version
and at the end thought the problem was due to the WEBP driver. But now
I'm getting the same type of errors with another driver and, ... one nuanc
Indeed that solved the linking issue, but while I am at this I have to
refer to still another problem.
For quite some time that build the mrsid driver (Geo_DSDK-7.0.0.2167)
but now it errors with.
However, shouldn't it be trying to link against lti_dsdk_dll.lib
instead of lti_dsdk.lib?
Agai
Hi,
I am having some troubles in writing a sub-region on an array with
GDALRasterIO.
The story is this. GMT internal grid arrays have a padding of 2 rows and
2 columns around the 'true data region' that is used for boundary
conditions but that need to be striped off upon saving on file.
Th
On 13-05-2012 20:28, Even Rouault wrote:
Le mardi 08 mai 2012 15:52:11, Joaquim Luis a écrit :
Hi,
I am having some troubles in writing a sub-region on an array with
GDALRasterIO.
The story is this. GMT internal grid arrays have a padding of 2 rows and
2 columns around the 'true data r
On 16-05-2012 02:50, William Kyngesburye wrote:
forget the directory problem - nothing should be linking the Carbon framework
in the first place. Especially on OSX 10.7, since carbon 64bit does not exist,
only 32bit. But then, that's just a warning, and may have nothing to do with
the actual
From it's placement in the above and a vague memory of the old ecw builds, I
suspect the carbon part comes from ECW. those separate libNCEcw* libraries are
an old way of building ECW, the newer method builds as a single libecwj2. ECW
is a PITA to compile on OSX (my notes still leave some u
On 16-05-2012 15:12, William Kyngesburye wrote:
On May 16, 2012, at 9:08 AM, William Kyngesburye wrote:
On May 16, 2012, at 8:50 AM, William Kyngesburye wrote:
I didn't mentioned explicitly but I did a svn update so the code should be very
close to what's in 1.9.1rc
Anyone building the relea
On 25-05-2012 14:40, Jeff McKenna wrote:
Hopefully this helps someone else down the road...
- the latest version of Poppler that I am able to get working with GDAL
(1.9.1 now) on MSVC (2008) is 0.16.4
- 0.20.0 or git master gives errors
- example of errors: http://pastebin.com/4YSnSJkn
So, stic
0.18.X and 0.20.X are indeed supported, but 0.19.X not. Odd numbered versionss
of poppler are developement versions, whereas even ones are considered stables.
I initially supported 0.19 in the hope of being reading for 0.20 release, but
there was a API break later in the 0.19 development, and beca
David,
Would this help?
https://sites.google.com/site/mironehowtos/satellite/import-a-modis-l2-sst-hdf-grid
I do the projection using the geolocation arrays and an internal
interpolation using a minimum curvature algorithm. Works pretty nicely
with MODIS L2 SSTs and chlorophyll
Joaquim
Ru
One other thing called my attention. If, as we can read on the wiki page
of OFGT
"... Most of the stand-alone programs use GDAL libraries and many of the
scripts rely heavily on GDAL command-line utilities."
than why the OFGT's licence is GPL?
I know MIT's license allows it, but it clearly l
intention to abuse any of the great work of the
gdal/ogr developers. Rather the opposite.
The license is being discussed as we speak and we're considering
various options. What would you suggest? Should it be MIT or would ,
for example, Modified BSD License do?
Best,
Anssi
On Wed, Jun 20,
Thanks again to Etienne for helping with this.
After closing the filed ticket [1] he has updated the wiki page [2] as
well with this, so most likely no one can be in my same shoes
p
[1] http://trac.osgeo.org/gdal/ticket/4717#comment:7
[2] http://trac.osgeo.org/gdal/wiki/NetCDF
Thanks for upda
ogram grdproject but I warn you already that it might not
be obvious.
Also you can try my front end to grdproject in Mirone. It should be easier but
again there will
probably be fine details that one need to adjust.
Joaquim Luis
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
ram is lightning
fast
in interpolating a grid into a series of x,y positions.
Joaquim Luis
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
Kenn Sebesta wrote:
If you release the shapefile constraint, GMT's grdtrack program is lightning
fast
in interpolating a grid into a series of x,y positions.
Well, the only reason I'm using shapefiles is because that's what the
mapping department gave me. If a shapefile can be converted into a
rm,
hTransformArg,
adfDstGeoTransform, &nPixels, &nLines,
adfExtent, 0 );
and the warping fails (the output is equal to input)
Is this a MEM driver issue, or an error of mine?
Thanks
Joaquim Luis
___
gdal
Even Rouault wrote:
Le Sunday 07 December 2008 21:30:17 Joaquim Luis, vous avez écrit :
Hi,
I have my gdalwarp MEX that works reasonably into converting between
projection systems. And now I would like to have it working to do warping
with GPS, but after hours of struggling with it, I don
Sven Geggus wrote:
Joaquim Luis <[EMAIL PROTECTED]> wrote:
I have my gdalwarp MEX that works reasonably into converting between
projection systems.
Why do you post this a a reply to my Question?
There is zero relation between my gdal2tiles Question and your stuff :(
Sven
Sorry,
Yo
trunk\gdal\ogr\ogrsf_frmts\mitab>cd gml && nmake /NOLOGO /f ma
kefile.vc && cd ..
The system cannot find the path specified.
...
Joaquim Luis
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
I think it doesn't matter
Joaquim Luis
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
s a side note, I also have MEXs that link with the OpenCV library, but those
do not sensibly
fragment the memory.
Thanks
Joaquim Luis
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
Frank Warmerdam wrote:
Joaquim Luis wrote:
So one question is, why warping with GCPs takes such large amount of
memory (apparently a chunk
of 500 Mb was not enough to process 5000 GCPs)?
Joaquim,
Is the memory fragmentation really related to the use of many GCPs? Does
the same problem
Joaquim,
OK. Well, generally speaking GDAL should allocate relatively little
memory just on loading the DLL (and perhaps calling GDALAllRegister()).
So I don't know why you are seeing the problem you see.
Frank,
It's ok. I was not expecting an easy answer but I think this issue worth
rep
rping memory problem was due to a
stupid error of over allocation - nGCPCount*nGCPCount instead of the needed
nGCPCount
It works fine now.
But the memory fragmentation issue is another story, and that still remains.
Joaquim Luis
___
gdal-dev mailing list
alled
NE2_land_only.jpg
than the world file name must be either NE2_land_only.jgw or NE2_land_only.wld
And, humm, there is much sea water on that "..._land_only" image
Joaquim Luis
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
Alex Bernstein wrote:
On Sun, Jan 4, 2009 at 11:00 PM, Joaquim Luis wrote:
Joaquim, renaming the world file to NE2_land_only.jgw was one of first
things I've done, sorry I forgot to mention that. So it's not reason
for the problems I'm having with reprojection.
The coordinate
Alex Bernstein wrote:
Joaquim,
Thanks for trying to help. I'm actually just interested in the map of
Africa from this global map. May be I'll try to figure out gdal
commands to crop and transform just the relevant region.
That will work fine. I just tried with Mirone and checked against GMT
co
Joaquim Luis wrote:
After that I wrote a MEX version of gdaltransform to be able to
reproject points using GCPs.
This MEX works but it revealed serious problems. At least under Matlab.
Matlab requires not only LOTS of memory but also that the available
memory is continuous.
For example, if
alinfo output of the two files they do not agree.
Why are the LR coords from the geotiff file different from the tif + .tfw?
And where do they come from?
Thanks
Joaquim Luis
# -- THE .TFW CONTENTS
677.678633953751# x_inc
0.0
0.0
-677.201365837256 # y_inc
-1389588.37368302
.
The fact that the pixels in the TFW are not exactly square is beyond our
control.
Recall that we start with a postscript image and want to end up with an
as-good-as-possible geotiff image.
Thanks
Joaquim
Joaquim Luis wrote:
(which executes this command)
gdalwarp -s_srs "+proj=stere
or a file containing the
WKT.
I'm not complaining (well on the contrary), but we may need to document this.
Since what version is it possible to provide a Proj4 string in a file for
-a_srs?
(I'm using the GDAL trunk)
Thanks
Joaquim Luis
___
g
I'm not complaining (well on the contrary), but we may need to
document this.
Since what version is it possible to provide a Proj4 string in a file
for -a_srs?
Joaquim,
Underlying the -a_srs function is the
OGRSpatialReference::SetFromUserInput()
method which does mention that proj.4 string
't know it yet have a look at the OpenCV library. It is amazingly
fast, but doesn't know anything about coordinates.
Joaquim Luis
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
5)
when dealing with global files
Joaquim Luis
I tried georeferencing one of these images ( http://www.oera.net/How2/TextureMaps2.htm ) using
MapWindow GIS; the pogram generated a .jgw file, which then I fixed by replacing "," by
"." and changing something like 179.12345 to
Brian,
The left side of your image contains a ... legend, so it's not part of the
map. How can any program guess what is legend and where the map starts?
Joaquim Luis
I am still battleing this conversion. I was just informed, via IRC,
that gdal won't overwrite existing files, so
d a
legend left on it. In my last post, I describe how I finally got it to
work. :)
Good for you,
... but are you sure it worked? I mean, are you sure that it is correctly
referenced?
At half way of the description I just quit following it (sorry, too complicated
for me).
Joaquim Luis
On S
on proximity to each other?
A job for GMT's blockmean or blockmedian.
Joaquim Luis
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
Okay, I downloaded the Manpages.
My initial reading is that blockmean and blockmedian work by averaging
the values.
blockmedian does not do point average. It does take the median point inside
each bin.
However, if the number of points inside a bin is even than it takes the average
of the two
Mateusz Loskot wrote:
for(Roadway::RoadWayArray::iterator itrw = _roadwayArr.begin(); itrw !=
_roadwayArr.end();itrw++)
++itrw;
if you care about performance.
Mateusz,
Just curious. Why should that impact on performance?
Joaquim Luis
___
gdal
Thankx
I'll keep this in mind if I ever one day start practicing C++
For the time being, only good old C.
Joaquim
Joaquim Luis wrote:
Mateusz Loskot wrote:
for(Roadway::RoadWayArray::iterator itrw = _roadwayArr.begin(); itrw !=
_roadwayArr.end();itrw++)
++itrw;
if you care
tool "ps2raster" (CVS version) is now able to
create KML code from any -- GMT-complicate-as-much-as-you-want -- images.
Joaquim Luis
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
t instance?
Thanks
Joaquim Luis
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
F to make it
a GeoPDF. Is it a Ghostscript command line?
Brent,
One only has to convert it to pdf using ghostscript.
I did it with GMT's ps2raster, like that
ps2raster quad.ps -Tf -A
Joaquim Luis
Brent
Klokan Petr Přidal wrote:
Hi,
There is a great blog post (and the linked "
we had a way to translate the GMT
projection syntax into
WKT, which we don't.
Joaquim
Joaquim,
Doing the format conversion from .ps to pdf is one thing (and there
are several ways to do it), but embedding the georeferencing in the PDF
to make it a GeoPDF is the interesting bit. I h
Peter J Halls wrote:
A bit of googling determines that the 'Adobe way' structure, together
with its elements, is defined in the ISO PDF definition - ISO 32000-1
(2008).
Looking at this with a PostScript Developer hat, the /LPTS has to be
associated with scaling / transforming the GPTS coordin
201 - 287 of 287 matches
Mail list logo