Dear list
I'm doing some analysis of orthophotos with i.gensigset and i.smap, but
suddenly when running i.gensigset i get this:
GRASS 6.1.cvs (ortho35):~ > i.gensigset trainingmap=2616trening
group=2616 subgroup=2616 signature=2616sign1
Finding training classes ... 100%
ERROR: G_calloc: out of memory
GRASS 6.1.cvs (ortho35):~ >
I'm running it on cygwin with 6.0.0., and I have tried to reinstall both
cygwin and GRASS 3 times but the same thing happens again.
In my prosesses window of the windows task manager the i.gensigset
appears two times counting different amounts of memory when runnig the
module. I'm in short of time so if anyone has an answer I'll be very
happy!!!!
Knut
<mailto:d@nina.no>
Dear Knut,
I have had correspondence with Charlie Bouman, the author of i.gensigset, about this issue. There is a problem with the code. As I understand it, for certain large data sets, it can't solve the fitting of the variance-covariance matrix. i.gensigset is prone to result in singular matrices..this is where I show my ignorance on exactly how it works....but as the dimensionality of the data increases there is a greater chance this will occur. On Charlie's website (http://dynamo.ecn.purdue.edu/~bouman/ ) is a new version, called CLUSTER, which has an algorithm that detects and deals with singularities. You must compile this code on your computer (there is a makefile; I am sure others on this list can recommend the best way to compile if this is new to you). There is also a user's manual and some example files. I dump the data out of grass using r.stats -1n input=grid1,grid2,grid3 output=cluster_input.txt
where grid1, grid2, grid3 etc are the 'bands' in the imagery subgroup that you wish to use.
The above r.stats command must be run for each class (I use a mask to subset the map into classes). CLUSTER takes the above text file and an information file (stating the file names for each class.) as input. The file created by CLUSTER can be put into the correct folder (in the proper group, subgroup, sigset folder) and be used with i.smap. You will need to edit this file, however. It requires different numbers at the end of the lines that start with
classnum:
Otherwise i.smap will classify everything as the same class. You also need to be sure that the file named REF in the subgroup lists your rasters in the same order they are used in CLUSTER.
Charlie Bouman also said he has made minor improvements to i.smap (called smap in the version on his website), but the changes to i.gensigset were substantial. I think it would be nice if someone who knows how can update i.gensigset using CLUSTER, assuming Dr. Bouman agrees.
Cheers,
Dan
At 1:45 PM +0200 6/6/05, knut.langeland@nina.no wrote:
Dear list
I'm doing some analysis of orthophotos with i.gensigset and i.smap, but suddenly when running i.gensigset i get this:
GRASS 6.1.cvs (ortho35):~ > i.gensigset trainingmap=2616trening group=2616 subgroup=2616 signature=2616sign1 Finding training classes ... 100%
ERROR: G_calloc: out of memory
GRASS 6.1.cvs (ortho35):~ >
I'm running it on cygwin with 6.0.0., and I have tried to reinstall both cygwin and GRASS 3 times but the same thing happens again. In my prosesses window of the windows task manager the i.gensigset appears two times counting different amounts of memory when runnig the module. I'm in short of time so if anyone has an answer I'll be very happy!!!!
Knut
<mailto:d@nina.no>
--
Dan Gavin
Research Associate
University of Vermont
Botany and Agricultural Biochemistry
120B Marsh Life Science
109 Carrigan Dr.
Burlington, VT 05405-0086
802.656.0297
hi,
does r.stats+CLUSTER produce somehow "better" results, than i.gensigset?
Jachym
On Mon, Jun 06, 2005 at 08:43:11AM -0400, Daniel Gavin wrote:
Dear Knut,
I have had correspondence with Charlie Bouman, the author of
i.gensigset, about this issue. There is a problem with the code. As
I understand it, for certain large data sets, it can't solve the
fitting of the variance-covariance matrix. i.gensigset is prone to
result in singular matrices..this is where I show my ignorance on
exactly how it works....but as the dimensionality of the data
increases there is a greater chance this will occur. On Charlie's
website (http://dynamo.ecn.purdue.edu/~bouman/ ) is a new version,
called CLUSTER, which has an algorithm that detects and deals with
singularities. You must compile this code on your computer (there is
a makefile; I am sure others on this list can recommend the best way
to compile if this is new to you). There is also a user's manual and
some example files. I dump the data out of grass using r.stats -1n
input=grid1,grid2,grid3 output=cluster_input.txt
where grid1, grid2, grid3 etc are the 'bands' in the imagery
subgroup that you wish to use.
The above r.stats command must be run for each class (I use a mask to
subset the map into classes). CLUSTER takes the above text file and
an information file (stating the file names for each class.) as
input. The file created by CLUSTER can be put into the correct
folder (in the proper group, subgroup, sigset folder) and be used
with i.smap. You will need to edit this file, however. It requires
different numbers at the end of the lines that start with
classnum:
Otherwise i.smap will classify everything as the same class. You
also need to be sure that the file named REF in the subgroup lists
your rasters in the same order they are used in CLUSTER.
Charlie Bouman also said he has made minor improvements to i.smap
(called smap in the version on his website), but the changes to
i.gensigset were substantial. I think it would be nice if someone
who knows how can update i.gensigset using CLUSTER, assuming Dr.
Bouman agrees.
Cheers,
Dan
At 1:45 PM +0200 6/6/05, knut.langeland@nina.no wrote:
>Dear list
>
>I'm doing some analysis of orthophotos with i.gensigset and i.smap,
>but suddenly when running i.gensigset i get this:
>
>GRASS 6.1.cvs (ortho35):~ > i.gensigset trainingmap=2616trening
>group=2616 subgroup=2616 signature=2616sign1
>Finding training classes ... 100%
>ERROR: G_calloc: out of memory
>GRASS 6.1.cvs (ortho35):~ >
>
>I'm running it on cygwin with 6.0.0., and I have tried to reinstall
>both cygwin and GRASS 3 times but the same thing happens again.
>In my prosesses window of the windows task manager the i.gensigset
>appears two times counting different amounts of memory when runnig
>the module. I'm in short of time so if anyone has an answer I'll
>be very happy!!!!
>
>Knut
>
><mailto:d@nina.no>
>
--
Dan Gavin
Research Associate
University of Vermont
Botany and Agricultural Biochemistry
120B Marsh Life Science
109 Carrigan Dr.
Burlington, VT 05405-0086
802.656.0297
--
Jachym Cepicky
e-mail: jachym.cepicky@centrum.cz
URL: http://les-ejk.cz
GPG: http://les-ejk.cz/gnupg_public_key/
Hi everybody... i need some help.
im going to rectify some ortho photos but i dont have the calibration report of the camera. Clifford J. Mugnier (thanks) said me that the camera that was used to take the photos is a "WILD Heerbrugg camera made in Switzerlan". I really dont have the time to ask for the report that i need but i want to try to rectify the photos anyway. Somebody said me that on the Markus's book there is a section about the matter but i neither have the time to find that book (i need to finish the work im doing as soon as possible). I only have the focus leght of the camera. I also have the DTM to use for rectify. Particullary i dont have coordinates of the fiducials marks that are required by the command i.ortho.photo, i want to try to rectify the images anyway. Is there a way to do that? Is there the possibility to use a similar report? is there anybody that can send me a similar camera calibration report? Is there a "place" where i can find sample information to use?
At the moment for me is not really important the accuracy of the work, i only need to finish it. How can i do?
thaks everybody for the help
Luca P.
Email.it, the professional e-mail, gratis per te: http://www.email.it/f
Sponsor:
Iscriviti Gratis al primo corso on-line di Web Marketing in Italia!
* Impara subito come promuoverti su internet - clicca qui
Clicca qui: http://adv.email.it/cgi-bin/foclick.cgi?mid=3210&d=6-6
I had a similar problem. Here's what I did:
1 my photos where in written in tiff with correct DPI, so I just
2 opened one in Gimp
3 measured distances between fiducials
4 asked Prof. Pitagoras for help, because I had only 4 fiducials (in corners) and the photo was fairly rotated
5 image center is 0,0; coords to the left are negative, positive to the right; [mm] is the unit required in i.ortho.photo
If your DPI isn't properly stored in the image but you know it, measure distances in pixels, translate to inch -> mm.
Maciek
----- Original Message ----- From: "Luca Penasa" <luca.penasa@email.it>
To: "GRASS Users list" <grasslist@baylor.edu>
Sent: Monday, June 06, 2005 8:11 PM
Subject: [GRASSLIST:7029] i.ortho photo help
Hi everybody... i need some help.
im going to rectify some ortho photos but i dont have the calibration report of the camera. Clifford J. Mugnier (thanks) said me that the camera that was used to take the photos is a "WILD Heerbrugg camera made in Switzerlan". I really dont have the time to ask for the report that i need but i want to try to rectify the photos anyway. Somebody said me that on the Markus's book there is a section about the matter but i neither have the time to find that book (i need to finish the work im doing as soon as possible). I only have the focus leght of the camera. I also have the DTM to use for rectify. Particullary i dont have coordinates of the fiducials marks that are required by the command i.ortho.photo, i want to try to rectify the images anyway. Is there a way to do that? Is there the possibility to use a similar report? is there anybody that can send me a similar camera calibration report? Is there a "place" where i can find sample information to use?
At the moment for me is not really important the accuracy of the work, i only need to finish it. How can i do?
thaks everybody for the help
Luca P.
--
Email.it, the professional e-mail, gratis per te: http://www.email.it/f
Sponsor:
Iscriviti Gratis al primo corso on-line di Web Marketing in Italia!
* Impara subito come promuoverti su internet - clicca qui
Clicca qui: http://adv.email.it/cgi-bin/foclick.cgi?mid=3210&d=6-6
On Mon, Jun 06, 2005 at 08:11:33PM +0200, Luca Penasa wrote:
Hi everybody... i need some help.
im going to rectify some ortho photos but i dont have the calibration
report of the camera. Clifford J. Mugnier (thanks) said me that the
camera that was used to take the photos is a "WILD Heerbrugg camera made
in Switzerlan". I really dont have the time to ask for the report that i
need but i want to try to rectify the photos anyway.
To obtain a good results it's recommended to get this report.
Even google/whatever can help...
Somebody said me
that on the Markus's book there is a section about the matter but i
neither have the time to find that book
Here you can find it:
- any (online) bookshop
- http://mpa.itc.it/grassbook2/
(i need to finish the work im
doing as soon as possible). I only have the focus leght of the camera. I
also have the DTM to use for rectify. Particullary i dont have
coordinates of the fiducials marks that are required by the command
i.ortho.photo, i want to try to rectify the images anyway. Is there a
way to do that? Is there the possibility to use a similar report? is
there anybody that can send me a similar camera calibration report? Is
there a "place" where i can find sample information to use?
What's the name/number of the camera?
At the moment for me is not really important the accuracy of the work, i
only need to finish it. How can i do?
thaks everybody for the help
Luca P.
Markus
Markus Neteler wrote:
What's the name/number of the camera?
This is the number i see on the photo:
/UAiG 6025 153.26/
Clifford J. Mugnier said me that:
"The camera is a WILD Heerbrugg camera made in Switzerland. The company is
now owned by Leica. The lens type apparently is a Universal Aviogon wide
angle lens. The calibrated focal length is as you stated."
Iv obtained a good work with the metod suggested by Maciek:
"I had a similar problem. Here's what I did:
1 my photos where in written in tiff with correct DPI, so I just
2 opened one in Gimp
3 measured distances between fiducials
4 asked Prof. Pitagoras for help, because I had only 4 fiducials (in corners) and the photo was fairly rotated
5 image center is 0,0; coords to the left are negative, positive to the right; [mm] is the unit required in i.ortho.photo
If your DPI isn't properly stored in the image but you know it, measure distances in pixels, translate to inch -> mm."
Thanks for your help.
Luca P.
Email.it, the professional e-mail, gratis per te: http://www.email.it/f
Sponsor:
Hai una fotocamera digitale e vuoi stampare le tue immagini
* su vera carta fotografica professionale?
Clicca qui: http://adv.email.it/cgi-bin/foclick.cgi?mid=1533&d=7-6