Irgendwie komme ich mit den Zahlen auf keinen grünen Zweig… Zum einen hatte ich offenbar PX und PY vertauscht, aber auch wenn ich die passend habe, bringt das nicht viel. Der Punkt ist, dass ich mit den angegebenen Skalenfaktoren grob auf 3,8 * 148 = 562 Einheiten (in der gewählten Projektion) pro Pixel komme, mit deinen Koordinaten der Eckpunkte aber nur auf ganz grob 20 Einheiten pro Pixel. Da kommt irgendwas nicht richtig hin.
meine ich auch, dass da was faul sein muß. Hab hoch ein paar neue Kalibrierungen vesucht und komme immer falsch hin.
Hier nochmals ein CAL-File aber zum selben Image. Hab es halt nochmals kalibriert und jetzt die dynamische Methode mit den 3 Referenzpunkten verwendet.
#JOSM PicLayer plugin calibration data
#Wed Feb 18 23:40:41 CET 2015
POSITION_Y=3144843.921914018
POSITION_X=9519879.519414948
M12=-50.87634699525655
M11=0.022148575397414223
M10=6.719567798201829E-5
M02=-23.772573686878683
INITIAL_SCALE=81985.45583742265
M01=-0.0010414742927757161
M00=0.022211044268432425
W = 1668
H = 2338
PX = POSITION_X
PY = POSITION_Y
S = INITIAL_SCALE/100
print PX,PY,S
X1 = [-W/2,-H/2, 1]
X1 = [0,0,1]
print X1
X2 = [[S, 0, 0],
[0, S, 0],
[0, 0, 1]]
print X2
np.dot(X2,X1)
X3 = [[M00, M01, M02],
[M10, M11,-M12],
[ 0, 0, 1]
]
print X3
X4 = [[ 1, 0, PX],
[ 0, 1, PY],
[ 0, 0, 1]
]
print X4
np.dot(np.dot(np.dot(X4,X3),X2),X1)
damit krieg in in Python das raus: array([ 9.51985575e+06, 3.14489480e+06, 1.00000000e+00])
und das ist fast POSITION_Y/POSITION_X: 9519879. 3144843.9
alles sehr merkwürdig.
Ich mach mal damit Schluss für heute. Dreh mich eh im Kreis.
Gruss
walter
hab zwar keine Lib installiert sondern nur ein interaktives math programm, aber dass da irgendwo 'ne Formel oder die Matrix faul sein muss, war mir inzwischen schon klar - nur wo?
Danke für dein Beispiel. Ich werde es am DO mal in aller Ruhe analysieren und meine Auswertung daran anpassen.
Gruss
walter
EDT: wow, hab es mal überflogen und bin beeindruckt.
Wenn die Formeln im Wiki so falsch sind, ist es kein Wunder, wenn da Murks rauskommt 
Grundsätzlich müsste die Formel aber stimmen. In welcher Reihenfolge man Rotation oder Skalierung durchführt, spielt keine Rolle. Wichtig ist, dass man die Translation zuletzt ausführt, weil man ansonsten nicht um das Zentrum des Bildes drehen würde.
Vertauscht? Wie wählen die bitte ihren Index?! Vorzeichen… das mache ich meistens durch ausprobieren (je nachdem, wie rum man projiziert).
Das wiederum macht Sinn.
Gut zu wissen. Hätte ja auch auf irgendeine Art rechtwinklig sein können. 
Ich werde den Teil im Wiki mal aktualisieren.
PS: Weiß jemand, wie man schnell an ein LaTeX-Plugin im Wiki kommen kann?
PPS: So besser?
Aber sicher! Soooo verstehe ich das natürlich sofort
Ich mach mich nachher dran, die Sache umzusetzten und dann sehen wir weiter.
Gruss
walter
Wer ist denn für das Wiki zuständig? Wäre wirklich super, wenn dort das (ein) LaTeX-Plugin hinzugefügt werden könnte. Mit den Bildern ist das nervig.
Matrixmultiplikation ist einfach:
Edit: Unfug getilgt.
Ich bin von dir immer wieder überrascht. Gibs zu, du hast schon längst eine Zeitmaschine gebaut, damit du das alles schaffst und immer noch Zeit für solche Spielereien hast. 
Leider immer noch daneben - wenn ich nicht völlig schief liege:
zuerst mal die Lage in Josm - ja, es ist Nepal:
Das Image habe ich mit dem PicLayer-Pluging referenziert. Leider hat es keinen transparenten Hintergrund, aber ich hab Josm so eingestellt, dass man hoffentlich die nach diesen Image erfasste Außengrenze von Lalitpur einigermaßen erkennen kann.
Danach hab ich das Programm von Netzwolf um eine weitere Ausgabe erweitert und bekomme jetzt das:
Bildmitte : x: 9513536 y: 3179064 lon: 85.461555 lat: 27.444086
Oben links : x: 9496946 y: 3202516 lon: 85.312518 lat: 27.630889
Oben rechts: x: 9530127 y: 3202516 lon: 85.610592 lat: 27.630889
Unten links : x: 9496946 y: 3155613 lon: 85.312518 lat: 27.256967
Unten rechts: x: 9530127 y: 3155613 lon: 85.610592 lat: 27.256967
insert into refnodes values(ST_SetSRID(ST_MakePoint( 9513536, 3179064),3857),'Bildmitte ');
insert into refnodes values(ST_SetSRID(ST_MakePoint( 9496946, 3202516),3857),'Oben links ');
insert into refnodes values(ST_SetSRID(ST_MakePoint( 9530127, 3202516),3857),'Oben rechts');
insert into refnodes values(ST_SetSRID(ST_MakePoint( 9496946, 3155613),3857),'Unten links ');
insert into refnodes values(ST_SetSRID(ST_MakePoint( 9530127, 3155613),3857),'Unten rechts');
Den zweiten Block konnte ich mit Cut&Paste in Postgresql/PostGIS eingeben. Zudem hab ich mit den Werten 9496946,3202516 (Oben links) “meine” Daten in Postgis umgerechnet und danach beides (also die Rohdaten vom Script als auch mein in PostGIS umgerechtete Image) mittels QGIS dargestellt:
Und ja: Qgis ist auf 3857 eingestellt.
Wie man leicht erkennt liegt das Bild noch ein wenig falsch, deckt sich aber total mit den vom Script errechneten Basiswerten. Daher können die Basiswerte auch nicht ok sein. Es müsste alles noch (um die halbe Bildgröße?) nach Nordwest verschoben werden.
Grübelnde Grüsse
walter
ps: so nah waren wir bisher nicht dran 
Ich weiß nicht, ob das noch was mit dem aktuellen Workflow zu tun hat, aber da wird der Wert der ersten Zeile (Pixelkoordinaten Bildmitte) in der nächsten Zeile durch die Koordinaten der Ecke (0,0) überschrieben.
Nö, das ist mausetot. Ich nehme nur noch den Script von Netzwolf.
Nur nochmals zur Sicherheit:
./josmPiclayerCalDecoder.pl 1654 2338 nepal_01_Admin.jpg-move.cal
mit
#JOSM PicLayer plugin calibration data
#Wed Feb 18 02:01:41 CET 2015
POSITION_Y=3132393.459019553
POSITION_X=9534742.422238246
M12=-50.68955846622368
M11=0.021788490422173075
M10=0.0
M02=-23.031323884902726
INITIAL_SCALE=81985.45583742265
M01=0.0
M00=0.021788490422173075
ergibt obiges Ergebnis.
Nahmd,
Ich hoffe, dass war ein sachgemäßer Eingriff, denn bei unsachgemäßen Eingriffen erlischt die Garantie! 
Liefert bei mir im JOSM folgendes Ergebnis:
http://www.netzwolf.info/tmp/madhyamanchal.jpg
Und die Anzeige passt zu dem vom Skript berechneten Ergebnissen:
$ /opt/www/www.netzwolf.info/htdocs/2015/josmPiclayerCalDecoder.pl 1654 2338 <rahmen.png.cal
Bildmitte : x: 9513536 y: 3179064 lon: 85.461555 lat: 27.444086
Oben links : x: 9496946 y: 3202516 lon: 85.312518 lat: 27.630889
Oben rechts: x: 9530127 y: 3202516 lon: 85.610592 lat: 27.630889
Unten links : x: 9496946 y: 3155613 lon: 85.312518 lat: 27.256967
Unten rechts: x: 9530127 y: 3155613 lon: 85.610592 lat: 27.256967
Folgerung: wir nutzen eine unterschiedliche Projektion.
(die speichert das &$%" Piclayer nicht in die Config).
Ich hab beim JOSM eingestellt:
Merkator / EPSG:3857
Gruß Wolf
nö, nur diese zusätzliche Ausgabe:
sub postgis {
my ($label, $imageX, $imageY) = @_;
my ($x, $y) = transform ($imageX, $imageY);
my ($lon, $lat) = epsg3857toDegrees ($x, $y);
printf "insert into refnodes values(ST_SetSRID(ST_MakePoint(%8d, %8d),3857),'%s');\n",
$x, $y, $label;
}
...
#------------------------------------------------------
# und jetzt noch ernster
#------------------------------------------------------
printf"\n";
postgis ("Bildmitte ", 0 , 0 );
postgis ("Oben links ", -$imageWidth/2, +$imageHeight/2);
postgis ("Oben rechts", +$imageWidth/2, +$imageHeight/2);
postgis ("Unten links ", -$imageWidth/2, -$imageHeight/2);
postgis ("Unten rechts", +$imageWidth/2, -$imageHeight/2);
und hier das Image: http://www.mofald.gov.np/districtmap/25Lalitpur/01_Admin.jpg
So langsam versteh ich die Welt nicht mehr - oder die mich nicht?
So, und jetzt hab ich mir durch ein verd. “rm * ~” den ganzen Kram zerschossen 
Melde mich irgendwann mal wieder.
Gruss
walter
Bingo!
ich musste ja alles nochmals von vorne durchspielen und jetzt funzt es. Wird wohl ein komisches Kalibrierungsfile gewesen sein. Die jetzt noch zu erkennenden geringfügigen Abweichungenen kommen von der Neukalibrierung, sind also nicht schlimm.
Danke nochmals an alle Helfer und Gruss
walter
Ich bekomme beim Kalibrieren des Bildes vom Lalitpur-Distrikt die folgende cal-Datei
#JOSM PicLayer plugin calibration data
#Thu Feb 19 18:27:34 CET 2015
POSITION_Y=3195594.011383851
POSITION_X=9497775.91890219
M12=54.696858132803214
M11=0.29141155222927567
M10=0.0034887618851810187
M02=2.717468729812156
INITIAL_SCALE=6124.726383092037
M01=-0.0034888386184205216
M00=0.29141176398165014
Die Abweichung in den Skalierungsfaktoren ist nicht weiter tragisch, die resultierende Skalierung S*M00 bzw M11 stimmt auf Promille genau überein.
Gravierend ist die Abweichung in PX und PY (Bildmittelpunkt in projizierten Koordinaten) zu Post #24
POSITION_Y=3132393.459019553
POSITION_X=9534742.422238246
Auch die Eckpunkte des Bildes
NW 3215180/9481030
SO 3168480/9514920
stimmen im Rahmen der Kalibriergenauigkeit gut mit den Angaben in #9 überein.
Alles in Mercator/3857.
Fazit: Schon beim Erzeugen des cal-Files ging was schief. Der Bildmittelpunkt in #24 passt nicht zu den Ecken in #9 (und bei mir).
Da ist noch das beschrieben, was im Wiki als alte Version des cal-Files bezeichnet wird (nur Skalierung, Translation und Rotation).
Die neuere berechnet aus den drei Skalierungspunkten eine affine Transformation mit Scherungen und linearen Verzerrungen und schreibt z.B. kein ANGLE mehr.