Wie der Titel schon verrät muss ich Grid Raster miteinander verschneiden. Das wäre an sich kein Problem wenn es sich nur um einige Dateien handeln würde. Ich muss allerdings einige Hundert miteinander verschneiden. Im hab eigentlich recht wenig Praxis was sowas angeht aber ich denke da komm ich um ein kleines C++ Programm nicht drum herum das ich dann für ArcInfo verwenden kann, weil das manuel nicht machbar ist.

Es handelt sich im Großen und Ganzen um Niederschlagsdaten die als Grid Raster vorliegen. Diese Daten sollen auf ein bestimmtes Gebiet (Einzugsgebiet, was auch als Grid Raster vorhanden ist) zugeschnitten werden.

Die Niederscglagsdaten sind in einzelene sogenannte "event"-Ordner unterteilt,dh es geht los bei "event4" dann "event9","event10" (die zahlen sind nicht durchgehend) usw.
In diesen "event"-Ordnern sind die Grid Raster Ordner untergebracht die nach "cMonatTagStundeMinute" benannt sind.

Wie gesagt nun gilt es diese Grid Raster auf den einen Gebiets Grid Raster zuzuschneiden.

Wenn einer vielleicht eine Idee hat wie ich an dieses Problem herangehn sollte und wie der Quellcode eines solchen Programms auschauen sollte wäre ich sehr froh.

danke

Mfg

marlouwe schrieb:

Hallo zusammen !

Wie der Titel schon verrät muss ich Grid Raster miteinander verschneiden. Das wäre an sich kein Problem wenn es sich nur um einige Dateien handeln würde...

ok, wenn das Verschneiden von zwei Grids kein Problem ist, so erläutere uns doch mal, wie das geht (Link ist auch ok) - und wir überlegen dann wie es mit vielen Grids funktionieren könnte.

Gruß
Werner

noch irgend ein problem vorhanden?

Also unter Grid Raster (oder Raster Grid) verstehe ich eine räumliche Verteiltung von (Niederschlags)daten die man visuell in ArcView , ArcGis, etc darstellen kann. Bei mir sind das Niederschlagssummen für 4 oder 5 min über ein grosses Gebiet. Ich brauch allerdings nur die Daten für ein kleineres Gebiet (dessen Räumliche Darstellung auch als Grid vorliegt). Ich möchte nun die Daten des kleineren Gebietes aus dem grösseren herausschneiden.
Ich hab seid längerem keine Praxis mehr in GIS also lasst es mich wissen wenn ich was falsch verstehe.

Werner Salomon schrieb:

ok, wenn das Verschneiden von zwei Grids kein Problem ist, so erläutere uns doch mal, wie das geht (Link ist auch ok) - und wir überlegen dann wie es mit vielen Grids funktionieren könnte.

Wie gesagt ich hab das schon länger nicht mehr gemacht, allerdings weiss ich dass 2 einzelne Grids in ArcGis zu verschneiden kein so grosses Problem ist. Der befehl dazu lautet:Clip:
ArcToolBox
Analysis Tools >> Extract >> Clip

siehe Link:http://ifgivor.uni-muenster.de/vorlesungen/Geoinformatik/uebungen/ArcGIS_Beispiele/Beispiel_09/AG_B_09.html

rapso schrieb:

ueber wieviel MB an daten sprechen wir hier?

Ich spreche hier von 50GB Daten, also jede Menge Grids oder Raster wie man das nun nennen mag.

http://support.esri.com/index.cfm?fa=knowledgebase.techarticles.articleShow&d=22136

und zwar unter: Clip a grid to a specific polygon shape.-show me-

das Problem ist halt dass ich das 5000 mal machen muss

die Links, die Du angegeben hast, erklären wie man das Verschneiden mit einem Tool namens "ArcGIS" oder so erledigt. Ok - was erwartest Du jetzt als Antwort? Du befindest Dich hier einem C++-Forum. Mal angenommen wir würden Dir hier jetzt den Code posten, der die von Dir beschrieben Aufgabe macht. Was nützt er Dir dann? Wie bekommst Du das in "ArcGIS" hinein?
Wenn Du ein eigenes Tool neben "ArcGIS" erstellen willst, dann müssen wir wiederum wissen, wie das Verschneiden vom Algorithmus her funktioniert (die Erklärung: wie bedient man ArcGIS nützt wenig), oder Du gibst uns einen Link auf eine C++-Bibliothek, die über entsprechende Funktionen verfügt.
Welches Format haben denn Deine Dateien? kannst Du mal einen typischen Auszug posten.

Oder anders gefragt. Was soll das von Dir ersehnte C++-Programm genau tun? was ist der Input und was soll der Output sein?

Gruß
Werner

Für alle die diese GIS Programme nicht kennen, es handelt sich hier um Software von der Firma ESRI, dazu gehören ua ArcView, ArcGis, ArcInfo etc

In ArcInfo kann man eben mit GIS Daten operieren wobei man sich vorher die dazu benötigte aml datei schreibt und die in ArcInfo einliest.

Hier zB eine aml bei der Grids miteinander multipiziert wurden:

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <iostream.h>

int i,filenumb;
int year,month,day,hour,min;//long fileli [600];
int date;
int event;

void main()
{

// 3 variables (output path ´fi1´, event number and date of the event) have to be changed to each event /////////

FILE *fi,*fi1;

fi = fopen ("out91_92.txt", "r"); //file must be in the same folder or path need to be changed

fi1 = fopen ("aml\\corr41.aml", "w"); //file name must be changed for each event

filenumb = 19854;

for (i=0 ; i<filenumb; i++)
{
fscanf (fi, "%d %d %d %d %d", &year,&month,&day,&hour,&min);
fscanf(fi,"\n");

date = month*1000000 + day*10000 + hour*100 + min;

//event number------------------------
event = 41;
//------------------------------------

//event date--------------------------
//if (date >= 10120800 && date <= 10150759) //event4
//if (date >= 11020800 && date <= 11060759) //event9
//if (date >= 11100800 && date <= 11110759) //event10
//if (date >= 11130800 && date <= 11140759) //event12
//if (date >= 11270800 && date <= 12070759) //event15
//if (date >= 12080800 && date <= 12100759) //event16
//if (date >= 12100800 && date <= 12110759) //event17
//if (date >= 12110800 && date <= 12150759) //event18
//if (date >= 12180800 && date <= 12230759) //event19
//if (date >= 12230800 && date <= 12300759) //event20
//if (date >= 12300800 && date <= 12320759 || //event21, also next line!!
// date >= 1010000 && date <= 1040759)
//if (date >= 1120800 && date <= 1240759) //event24
//if (date >= 1290800 && date <= 2150759) //event26
//if (date >= 2150800 && date <= 2170759) //event27
//if (date >= 2190800 && date <= 3030759) //event28
//if (date >= 3040800 && date <= 3060759) //event30
//if (date >= 3120800 && date <= 3130759) //event31
//if (date >= 3140800 && date <= 3160759) //event32
//if (date >= 3200800 && date <= 3240759) //event33
//if (date >= 4120800 && date <= 4150759) //event37
if (date >= 4200800 && date <= 4210759) //event41
//------------------------------------

{
fprintf(fi1, "grid\n");
fprintf(fi1, "setwindow ..\\radar199110\\g10121338\n");
fprintf(fi1, "setcell ..\\radar199110\\g10121338\n");

fprintf(fi1, "c");
if (date < 10000000)
{ fprintf(fi1, "0");}
fprintf(fi1, "%d", date);
fprintf(fi1, " = ");

fprintf(fi1, "..\\radar");
fprintf(fi1, "%d", year);
if (month < 10)
{ fprintf(fi1, "0");}
fprintf(fi1, "%d", month);
fprintf(fi1, "\\g");
if (date < 10000000)
{ fprintf(fi1, "0");}
fprintf(fi1, "%d", date);

fprintf(fi1, " * ");

fprintf(fi1, "..\\bias_grids\\e");
fprintf(fi1, "%d", event);

fprintf(fi1, "\n");

fprintf(fi1, "q\n");

}

} // von for

fclose(fi);
fclose(fi1);

} // von main

Der Input ist eine Datei 'out91_92.txt' mit Datums- und Zeit-Angaben. Etwa in der Form

Jahr Monat Tag Stunde Minute
2004 4 20 8 10
2004 4 20 8 11
...

Output wäre also eine aml-Datei - etwa in dieser Form

grid
setwindow ..\radar199110\g10121338
setcell ..\radar199110\g10121338
c04200810 = ..\radar200404\g04200810 * ..\bias_grids\e41
q
grid
setwindow ..\radar199110\g10121338
setcell ..\radar199110\g10121338
c04200811 = ..\radar200404\g04200811 * ..\bias_grids\e41
q

.. mit hunderten von Einträgen.

Ich habe Deinen Code (das war C-Code) mal in C++ übersetzt und ein bisschen Fehlerbehandlung eingebaut.

#include <iostream>
#include <iomanip>
#include <fstream>
#include <string>
#include <sstream>

int main()
{
    using namespace std;
    //event number------------------------
    const int event_ = 41;
    //------------------------------------

    // 3 variables (output path ´fi1´, event number and date of the event) have to be changed to each event /////////

    const char* in_fn = "out91_92.txt";
    ifstream in( in_fn );       // file must be in the same folder or path need to be changed
    if( !in.is_open() )
    {
        cerr << "Fehler beim Oeffnen von " << in_fn << endl;
        return -1;
    }
    ostringstream out_fn_buf;
    out_fn_buf << "aml\\corr" << event_ << ".aml";   // Dateiname: "aml/corr<event#>.aml"
    ofstream out( out_fn_buf.str().c_str() );
    if( !out.is_open() )
    {
        cerr << "Fehler beim Oeffnen von " << out_fn_buf.str() << endl;
        return -1;
    }

    const int filenumb = 19854;
    for( int i=0 ; i<filenumb; ++i )
    {
        int year,month,day,hour,min;    //long fileli [600];
        if( !(in >> year >> month >> day >> hour >> min).ignore( 9999, '\n' ) )
        {
            cerr << "Fehler beim Lesen des " << (i+1) << ".Datums in der Datei " << in_fn << endl;
            return -1;
        }
        string date;
        {
            // Format: MMDDhhmm
            ostringstream date_buf;
            date_buf.fill( '0' );
            date_buf << setw(2) << month << setw(2) << day << setw(2) << hour << setw(2) << min;
            date = date_buf.str();
        }
        //           MMDDhhmm             MMDDhhmm
        if( date >= "04200800" && date < "04210800" ) // event41: 20.04 8:00 bis 21.04 7:59
        {
            out << "grid" << endl;
            out << "setwindow ..\\radar199110\\g10121338" << endl;
            out << "setcell ..\\radar199110\\g10121338" << endl;

            const char fill_char = out.fill( '0' ); // setw(2) << month
            out << "c" << date << " = ..\\radar" << year << setw(2) << month << "\\g" << date 
                << " * ..\\bias_grids\\e" << event_ << endl;
            out.fill( fill_char );
            out << "q" << endl;     // quit
        }
    } // von for
    return 0;
}

Aus 'date' sollte man noch eine Klasse machen, aber das ist erstmal nicht so wichtig.

Bleiben noch 'ne Menge Fragen: Wie muss dieser Code jetzt verändert oder erweitert werden, um Deinen Anforderungen zu genügen?
Liegt Dir die passende Datei 'out91_92.txt' vor? Musst Du Directories lesen, oder hast Du eine Liste der Dateien? Benötigst Du eine Funktion, die aus einer Event-Nummer das Datums-Intervall erzeugt?
Woher kennst Du die Anzahl der Dateien (hier filenumb = 19854)? Wie muss der Ausgabe-Pfad (welcher Pfad genau ist gemeint) mit der Event-Nummer geändert werden?

Kennst Du einen Link, wo so ein Verschneiden-Algorithmus erklärt ist?

Gruß
Werner

Wenn ich das recht verstehe, so soll das C++-Programm nur ein aml-File erzeugen, was irgendeine Art Makro ist, das von ArcInfo gelesen wird und dann Dir die Klickerei der hunderte Dateien abnimmt.

Absolut korrekt !! output ist also eine AML Datei (Arc Macro Language) die in ArcInfo eingelesen werden kann

Werner Salomon schrieb:

Der Input ist eine Datei 'out91_92.txt' mit Datums- und Zeit-Angaben. Etwa in der Form
[code]Jahr Monat Tag Stunde Minute
2004 4 20 8 10
2004 4 20 8 11
...

Absolut korrekt !! Die datei ist vorhanden und sieht volgendermaßen aus:

von

1991 10 12 13 33
1991 10 12 13 38
1991 10 12 13 42
1991 10 12 13 47
1991 10 12 13 52

etc bis

1992 04 22 10 16
1992 04 22 10 21
1992 04 22 10 26
1992 04 22 10 31
1992 04 22 10 35

Ich hab hier noch von einem Vorgänger (der es von einem vorgäger hat, und alle nicht so richtig bescheid wissen) einen code , der was ähnliches gemacht hat.

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <iostream.h>

int i,filenumb;
int year,month,day,hour,min;
int date;

void main()
{
	FILE	*fi,*fi1;

	fi = fopen ("out02_03.txt", "r");
	fi1 = fopen ("2003\\gridclip.aml", "w");	//number must be changed to correct date; 3 times in this program

	filenumb = 20000; 

	for (i=0 ; i<filenumb; i++)
	{
		fscanf (fi, "%d %d %d %d %d", &year,&month,&day,&hour,&min);
		fscanf(fi,"\n");

		date = month*1000000 + day*10000 + hour*100 + min;

	if (year == 2003 && month == 03 && day >= 02 && day <= 30)		// if #1	//number must be changed to correct date;
	{

		fprintf(fi1, "gridclip");
		fprintf(fi1, " ");
		fprintf(fi1, "g");

		if (date < 10000000) 
			{ fprintf(fi1, "0");}
		fprintf(fi1, "%d", date);

		fprintf(fi1, " ");
		fprintf(fi1, "cg");

		if (date < 10000000) 
			{ fprintf(fi1, "0");}
		fprintf(fi1, "%d", date);

		fprintf(fi1, " ");
		fprintf(fi1, "227178");

		fprintf(fi1, " ");
		fprintf(fi1, "710761");

		fprintf(fi1, " ");
		fprintf(fi1, "249028");

		fprintf(fi1, " ");
		fprintf(fi1, "727111\n");

		fprintf(fi1, "gridascii");
		fprintf(fi1, " ");
		fprintf(fi1, "cg");
		if (date < 10000000) 
			{ fprintf(fi1, "0");}
		fprintf(fi1, "%d", date);

		fprintf(fi1, " ");
		if (date < 10000000) 
			{ fprintf(fi1, "0");}
		fprintf(fi1, "%d\n", date);

	}	// von if
	}	// von for

	fclose(fi);
    fclose(fi1);

}	// von main

Ich weiss dass er die gleiche Art von Daten wie ich benutzt hat, halt eben für ein anderes Jahr (2002/03) aber ich weiss nicht wie die genau hiessen oder in welchen Ordnern sie abgespeichert waren.

Um nochmal meine datenlage zu erklären:

Diese vielen 5min Niederschlagsdaten (Rasterdaten) sind zusammengefasst in einzelne Niederschlagsereignisse genannt: "eventx". Also sagen wir wir haben ein Niederschlagsereignis genannt "event9" , dann haben wir alle 5 min ein Bild (Raster) wärend des Ereignisses. Also im Folder "event9" gehts zB los mit Raster "c11020801", "c11020806", "c11020810", etc bis "c11051402" also die zahlen stehn für "c,Monat,Tag,Stunde,Minute". Event9 geht also vom 2 Nov um 8:01 bis 5 Nov um 14:02.

Werner Salomon schrieb:

Musst Du Directories lesen, oder hast Du eine Liste der Dateien?

Benötigst Du eine Funktion, die aus einer Event-Nummer das Datums-Intervall erzeugt?

kann das die Lösung sein ?:

//event number------------------------
		event = 41;
//------------------------------------

//event date--------------------------
		//if (date >= 10120800 && date <= 10150759)	//event4
		//if (date >= 11020800 && date <= 11060759)	//event9
		//if (date >= 11100800 && date <= 11110759)	//event10
		//if (date >= 11130800 && date <= 11140759)	//event12
		//if (date >= 11270800 && date <= 12070759)	//event15
		//if (date >= 12080800 && date <= 12100759)	//event16
		//if (date >= 12100800 && date <= 12110759)	//event17
		//if (date >= 12110800 && date <= 12150759)	//event18
		//if (date >= 12180800 && date <= 12230759)	//event19
		//if (date >= 12230800 && date <= 12300759)	//event20
		//if (date >= 12300800 && date <= 12320759 ||	//event21, also next line!!
		//		date >= 1010000 && date <= 1040759)
		//if (date >= 1120800 && date <= 1240759)	//event24
		//if (date >= 1290800 && date <= 2150759)	//event26
		//if (date >= 2150800 && date <= 2170759)	//event27
		//if (date >= 2190800 && date <= 3030759)	//event28
		//if (date >= 3040800 && date <= 3060759)	//event30
		//if (date >= 3120800 && date <= 3130759)	//event31
		//if (date >= 3140800 && date <= 3160759)	//event32
		//if (date >= 3200800 && date <= 3240759)	//event33
		//if (date >= 4120800 && date <= 4150759)	//event37
		if (date >= 4200800 && date <= 4210759)	//event41
//------------------------------------

Werner Salomon schrieb:

Woher kennst Du die Anzahl der Dateien (hier filenumb = 19854)?

Weiss ich nicht genau, hab das einfach so übernommen, an der Anzahl der dateien hat sich ja nichts geändert

Werner Salomon schrieb:

Wie muss der Ausgabe-Pfad (welcher Pfad genau ist gemeint) mit der Event-Nummer geändert werden?

HAb ich auch keinen blassen schimmer

Gruss

.. HAb ich auch keinen blassen schimmer

das ist zuwenig um daraus ein Programm zu machen

Damit fragst Du uns hier mehr oder weniger in welchen Directories und mit welchen Filenamen die Daten untergebracht sind, die bei Dir auf der Platte (oder Server) liegen.
Das musst Du schon selber wissen, welche Dateien das sind, die zu einem Raster verschnitten werden sollen (die eines Events?). Das hat weder was mit C++ noch was mit 'wie benutzte ich ArcInfo' zu tun, sondern mit der Datenhaltung bei Dir.

Versuche noch mal selber heraus zu bekommen welche Datei-Namen in eine aml-Datei hineinkommen. Wenn Du es mit Worten vollständig erklären kannst, sollte es kein Problem sein, ein C++-Programm daraus zu machen.

Gruß
Werner

Wie muss der Ausgabe-Pfad (welcher Pfad genau ist gemeint) mit der Event-Nummer geändert werden?

Ich weiss nicht ob ich die Frage richtig verstanden habe, aber ich denke ich brauch für jedes Event eine aml Datei.

wie gesagt

marlouwe schrieb:

Also im Folder "event9" gehts zB los mit Raster "c11020801", "c11020806", "c11020810", etc bis "c11051402" also die zahlen stehn für "c,Monat,Tag,Stunde,Minute". Event9 geht also vom 2 Nov um 8:01 bis 5 Nov um 14:02.

Es sind alles in allem 20 events wenn ich mich recht erinner(bin im Moment nicht im Büro, werd das aber morgen nachprüfen). Ein Event kann locker 3 Tage dauern, und ich hab alle 5 min ein Rasterbild , dass können ganz schnell eine ganze Menge werden. Verschnitten werden müssen also alle Rasterbilder von alle Events mit jeweils immer dem selben einen Rasterbild des Einzugsgebiets. Ich versuch später noch eine Skizze anzufertigen.