package producer import ( "bufio" "fmt" "os" "path/filepath" "strings" "github.com/airbusgeo/godal" log "github.com/sirupsen/logrus" "gocv.io/x/gocv" ) func (r *Registrator) SaveOriginalPanToGDALGTiff(tiffFile string) error { err := savePanToGDALGTiff(r.PanImage, 0, 0, tiffFile) if err != nil { return err } return nil } func savePanToGDALGTiff(pan gocv.Mat, topLeftX, topLeftY float64, tiffFile string) error { // log.Println("Saving PAN image to TIFF file:", tiffFile) width := pan.Cols() height := pan.Rows() ds, err := godal.Create(godal.GTiff, tiffFile, 1, godal.UInt16, width, height) if err != nil { log.Error("Error creating TIFF file: ", err) return err } defer ds.Close() setGeoTransform(ds, topLeftX, topLeftY, 1.25) ds.SetMetadata("NBITS", "16") // 将通道的数据转换为uint16数组 data := make([]uint16, width*height) for y := 0; y < height; y++ { for x := 0; x < width; x++ { data[y*width+x] = uint16(pan.GetShortAt(y, x)) } } band := ds.Bands()[0] band.SetColorInterp(godal.CIGray) err = band.IO(godal.IOWrite, 0, 0, data, width, height, godal.PixelSpacing(2), godal.LineSpacing(width*2)) if err != nil { log.Error("Failed to write data to band:", err) return err } log.Info("Saved pan image to ", tiffFile) return nil } func (r *Registrator) SaveRegisteredMssToGDALGTiff(tiffFile string) error { r.rgbirImage = gocv.NewMat() gocv.Merge(r.registeredMssImages[:], &r.rgbirImage) err := SaveBGRToGDALGTiff(r.rgbirImage, 4, 0, 0, 5, []godal.ColorInterp{godal.CIBlue, godal.CIGreen, godal.CIRed, godal.CIUndefined}, tiffFile) if err != nil { return err } return nil } func SaveBGRToGDALGTiff(bgr gocv.Mat, bands int, topLeftX, topLeftY float64, resolution float64, colorInterps []godal.ColorInterp, tiffFile string) error { width := bgr.Cols() height := bgr.Rows() // 创建一个二维切片来存储图像数据 data := make([][]uint16, bands) for i := range data { data[i] = make([]uint16, width*height) } // 从gocv.Mat中提取数据 for y := 0; y < height; y++ { for x := 0; x < width; x++ { for b := 0; b < bands; b++ { data[b][y*width+x] = uint16(bgr.GetShortAt(y, x*bands+b)) } } } ds, err := godal.Create(godal.GTiff, tiffFile, bands, godal.UInt16, width, height) if err != nil { log.Error("Error creating TIFF file: ", err) return err } defer ds.Close() // ds.SetMetadata("NBITS", "16") setGeoTransform(ds, topLeftX, topLeftY, resolution) for b := 0; b < bands; b++ { band := ds.Bands()[b] band.SetColorInterp(colorInterps[b]) err := band.IO(godal.IOWrite, 0, 0, data[b], width, height, godal.PixelSpacing(2), godal.LineSpacing(width*2)) if err != nil { log.Error("Failed to write data to band:", err) return err } } log.Info("Saved BGR MSS to ", tiffFile) return nil } func (r *Registrator) BytesToRaw(mssData []byte, filePath string) error { f, err := os.OpenFile(filePath, os.O_CREATE|os.O_TRUNC|os.O_WRONLY, 0777) if err != nil { return err } w := bufio.NewWriter(f) w.Write(mssData) return nil } func (r *Registrator) SaveRegisteredMssToRaw(raw string) error { f, err := os.OpenFile(raw, os.O_CREATE|os.O_TRUNC|os.O_WRONLY, 0777) if err != nil { return err } var mssData [4][]byte for i := 0; i < MssBands; i++ { mssData[i] = r.registeredMssImages[i].ToBytes() } width := r.registeredMssImages[0].Cols() * PixelBytes height := r.registeredMssImages[0].Rows() log.Println("Writing registered MSS to RAW file...", len(mssData[0])*4) log.Println("width:", width) log.Println("height:", height) w := bufio.NewWriter(f) for row := 0; row < height; row++ { w.Write(mssData[0][row*width : (row+1)*width]) w.Write(mssData[1][row*width : (row+1)*width]) w.Write(mssData[2][row*width : (row+1)*width]) w.Write(mssData[3][row*width : (row+1)*width]) } hdr := EnviHdr{ Samples: 4 * width / PixelBytes, Lines: height, Bands: 1, } hdrFile := filepath.Join(filepath.Dir(raw), fmt.Sprintf("%s.HDR", strings.TrimSuffix(filepath.Base(raw), filepath.Ext(raw))), ) hdr.Write(hdrFile) return nil } func (r *Registrator) Report() error { return WriteReport(&r.report, r.Params.ReportFile) }