OpenStreetMap

h4ck3rm1k3's Diary

seextractor on osm

Posted by h4ck3rm1k3 on 21 August 2009 in English.

http://astroa.physics.metu.edu.tr/MANUALS/sextractor/

What is SExtractor?
SExtractor (Source-Extractor) is a program that builds a catalogue of objects from an astronomical image. It is particularly oriented towards reduction of large scale galaxy-survey data, but it also performs well on moderately crowded star fields.

convert landsat_4000_2000_xy_86000_168000.png landsat_4000_2000_xy_86000_168000.fits

sextractor landsat_4000_2000_xy_86000_168000.fits

nice...

Discussion

Comment from h4ck3rm1k3 on 21 August 2009 at 22:09

## copyright 2009 james michael dupont
## AGPL 3.0
## converts the CATS results into OSM format.
## example usage :
# convert landsat_4000_2000_xy_86000_168000.png -depth 16 landsat_4000_2000_xy_86000_168000.fits
# sextractor landsat_4000_2000_xy_86000_168000.fits
# perl ./convert.pl 86000 168000 test.cat > test.osm

use strict;
use warnings;

# bottom_left
my $basex = shift @ARGV;
my $basey = shift @ARGV;
my $filename = shift @ARGV;

my $tilesize = 2000;

my $topleft = $basey - $tilesize;

my $resolution = 4000;

print '';

my $id = 0;
open IN, $filename;
while ()
{
next if /\#/;
my ($test,$X_IMAGE,$Y_IMAGE,$X2,$Y2) = split (/\s+/);

print "\n";

$id --;
# print join ("!",($NUMBER,$FLUXERR_ISO,$FLUX_AUTO,$FLUXERR_AUTO,$X_IMAGE,$Y_IMAGE,$FLAGS)) . "\n";

my $newy = ($topleft + $Y_IMAGE) / $resolution;
my $newx = ($basex + $X_IMAGE) / $resolution;
my $newy2 = ($topleft + $Y2) / $resolution;
my $newx2 = ($basex + $X2) / $resolution;

# print "Data:" . $Y_IMAGE . "/" . $X_IMAGE . "\n";
print "\n";
$id --;
print "\n";
$id --;
print "\n";
$id --;
print "\n";

$id --;

my $id1 = $id + 4;
my $id2 = $id + 3;
my $id3 = $id + 2;
my $id4 = $id + 1;
print "

";

#print $newy . "/" . $newx . "\n";

}

print "";
close IN;

Comment from h4ck3rm1k3 on 21 August 2009 at 22:17

well, of course you need to swap lat and long... doh.