childwin.pas

在delphi7 下开发医学图象浏览器,对医学图象进行编辑,分析的

//**************************************************************************
  VectorGraph := TVectorGraph.Create;
  ScrollBox1.DoubleBuffered := True; //防止画图时闪烁
  //设置为浏览状态
  PicState:=True;
  GRID_SPACE:=1;
//**************************************************************************
end;

procedure TMDIChild.DetermineZoom;
//Work out scale of image.
//Can scale image to fit the form size
var
  lHZoom: single;
  lZoom, lZoomPct: integer;
begin
  if (not MainForm.BestFitItem.checked) then
    exit;
  lHZoom := (ClientWidth) / g100pctImageWid;
  if ((ClientHeight) / g100pctImageHt) < lHZoom then
    lHZoom := ((ClientHeight) / g100pctImageHt);
  lZoomPct := trunc(100 * lHZoom);
  if lZoomPct < 11 then
    lZoom := 10 //.5 zoom
  else if lZoomPct > 500 then
    lZoom := 500
  else
    lZoom := lZoomPct;
  gZoomPct := lZoom;
end;

procedure TMDIChild.AutoMaximise;
//Rescales image to fit form
var
  lZoom: integer;
begin
  if (not MainForm.BestFitItem.checked) or (g100pctImageHt < 1) or (g100pctImageWid < 1) then
    exit;
  lZoom := gZoomPct;
  DetermineZoom;
  if lZoom <> gZoomPct then
  begin
    RefreshZoom;
    MainForm.ZoomSlider.Position := lZoom;
  end;
end;

function FSize(lFName: string): longint;
var
  infp: file;
begin
  assign(infp, lFName);
  FileMode := 0; //Read only
  Reset(infp, 1);

  result := FileSize(infp);
  closefile(infp);
  Filemode := 2;
end;

function TMDIChild.LoadData(lFileName: string; lAnalyze, lECAT, l2dImage, lRaw: boolean): Boolean;
//Loads and displays a medical image, also searches for other DICOM images in the same folder as the image you have opened
//  This latter feature allows the user to quickly cycle between successive images
var
  lHdrOK: boolean;
  lAllocSLiceSz, lS: integer;
  lExt: string;
  JPG: TJPEGImage;
  Stream: TmemoryStream;
  BMP: TBitmap;
  lImage: TImage;
begin
  ReleaseDICOMmemory;
  RescaleClear;
  gFilePath := extractfilepath(lFileName);
  FDICOM := true;
  gScheme := 1;
  gSlice := 1;
  LoadColorScheme('', gScheme); //load Black and white
  Result := TRUE;
  gImgOK := false;
  FFileName := lFileName;
  gAbort := true;
  if not fileexists(lFilename) then
  begin
    result := false;
    showmessage('Unable to find the file: ' + lFilename);
    exit;
  end;
  Self.caption := extractfilename(lFilename);
  lExt := UpperCase(ExtractFileExt(FFileName));
  if (l2DImage) or ('.JPG' = lExt) or ('.JPEG' = lExt) or ('.BMP' = lExt) then
  begin
    FDICOM := false;
    if ('.JPG' = lExt) or ('.JPEG' = lExt) then
    begin
      {JPEGOriginal := TJPEGImage.Create;
      TRY
         JPEGOriginal.LoadFromFile(FFilename);
         Image.Picture.Graphic := JPEGOriginal
      FINALLY
          JPEGOriginal.Free
      END;}
      //the following longer method makes sure the user can save the JPEG file...
      {Stream := TMemoryStream.Create;
      try
         Stream.LoadFromFile(FFilename);
         Stream.Seek(0, soFromBeginning);
         Jpg := TJPEGImage.Create;
         try
            Jpg.LoadFromStream(Stream);
            BMP := TBitmap.create;
            try
               BMP.Height := JPG.Height;
               BMP.Width := JPG.Width;
               BMP.PixelFormat := pf24bit;
               BMP.Canvas.Draw(0,0, JPG);
               Image.Picture.Graphic := BMP;
            finally
                   BMP.Free;
            end;
         finally
                JPG.Free;
         end;
      finally
             Stream.Free;
      end;{}
//next bit allows contrast adjustable JPEG images....
      Stream := TMemoryStream.Create;
      try
        Stream.LoadFromFile(FFilename);
        Stream.Seek(0, soFromBeginning);
        Jpg := TJPEGImage.Create;
        try
          Jpg.LoadFromStream(Stream);
          gDICOMData.XYZdim[1] := JPG.Width;
          gDICOMData.XYZdim[2] := JPG.Height;
          Image.Width := JPG.Width;
          Image.Height := JPG.Height;

          {                 BMP := TBitmap.create;
                           try
                              BMP.Height := JPG.Height;
                              BMP.Width := JPG.Width;
                              BMP.PixelFormat := pf24bit;
                              BMP.Canvas.Draw(0,0, JPG);
                              Image.Picture.Graphic := BMP;
                           finally
                                  BMP.Free;
                           end;}
        finally
          JPG.Free;
        end;
      finally
        Stream.Free;
      end;

      gDICOMData.SamplesPerPixel := 3;
      gDICOMData.Storedbits_per_pixel := 8;
      gDICOMData.Allocbits_per_pixel := 8;
      gDICOMData.ImageStart := 0;
      g100PctImageWid := gDICOMData.XYZdim[1];
      g100PctImageHt := gDICOMData.XYZdim[2];
      gDICOMData.XYZdim[3] := 1;
      gECATposra[1] := 0;
      //CloseFile(infp); //we will read this file directly
      lAllocSLiceSz := (gDICOMdata.XYZdim[1] * gDICOMdata.XYZdim[2]);
        //24bits per pixel: number of voxels in each colour plane
      gBuff24Sz := lAllocSliceSz * 3;
      GetMem(gBuff24, lAllocSliceSz * 3);
      decompressJPEG24(FFilename, gBuff24, lAllocSliceSz, gECATposra[1], Image);
      MainForm.ColUpdate;
      DetermineZoom;
      SetDimension(gDIcomData.XYZdim[2], gDIcomData.XYZdim[1], 24, gBuff24, false);
      DICOMImageRefreshAndSize;
      Image.Refresh;
      FDICOM := true;
      gImgMin := 0;
      gImgMax := 255;
      gWinMin := 0;
      gWinMax := 255;

      {}
      //?? what if gDICOMdata.monochrome = 4 -> is YcBcR photometric interpretation dealt with by the JPEG comrpession or not? I have never seen such an image, so I guess this is an impossible combination
      //Reset(infp, 1); //other routines expect this to be left open

    end
    else
      Image.Picture.Bitmap.LoadFromFile(FFilename);
    //if Image.Picture.Bitmap.PixelFormat = pf8 then
    gDICOMData.SamplesPerPixel := 3;
    gDICOMData.Storedbits_per_pixel := 8;
    gDICOMData.Allocbits_per_pixel := 8;
    gDICOMData.ImageStart := 54;
    gDICOMData.XYZdim[1] := Image.Picture.Width;
    gDICOMData.XYZdim[2] := Image.Picture.Height;
    g100PctImageWid := gDICOMData.XYZdim[1];
    g100PctImageHt := gDICOMData.XYZdim[2];
    Image.Width := Image.Picture.Width;
    Image.Height := Image.Picture.Height;
    gDICOMData.XYZdim[3] := 1;
    if self.WindowState <> wsMaximized then
    begin
      self.ClientHeight := gDICOMdata.XYZdim[2];
      self.ClientWidth := (gDICOMData.XYZdim[1]);
    end;
    MainForm.ColUpdate;
    ContrastAutobalance1.enabled := false;
    OptionsImgInfoItem.enabled := false;
    gImgOK := true;
    automaximise;
    Image.Refresh;
    exit;
  end;
  FDICOM := true;
  if lRaw then
  begin
    lHdrOK := true;
    gImgOK := true;
  end
  else if lAnalyze then
    OpenAnalyze(lHdrOK, gImgOK, gDynStr, FFileName, gDicomData)
  else if lECAT then
    read_ecat_data(gDICOMdata, true {verbose}, true {offset tables supported}, lHdrOK, gImgOK, gDynStr, FFileName)
  else
    read_dicom_data(true, true, true, true, true, true, true, gDICOMdata, lHdrOK, gImgOK, gDynStr, FFileName);
  rescaleInit;

  if gDICOMdata.ElscintCompress then
  begin
    showmessage('Unable to descode Elscint compressed images.');
    gImgOK := false;
  end;
  if (lHdrOK) and (gImgOK) and (not fileexists(FFileName)) then
  begin
    MainForm.OpenDialog.Title := 'Select Interfile image file...';
    if MainForm.OpenDialog.Execute and Fileexists(MainForm.OpenDialog.Filename) then
      FFilename := MainForm.OpenDialog.Filename
    else
      gImgOK := false;
    MainForm.OpenDialog.Title := 'Open';
  end; //1.33
  HdrShow;
  if gECATJPEG_table_entries > 0 then
  begin
    //showmessage('ecatabba'+inttostr(gDICOMdata.CompressOffset));
    //gDicomData.ImageStart := gECATJPEG_pos_table[1];
    {if (gECATJPEG_table_entries = 1) then begin
            gECATposra[1]:=gDICOMdata.CompressOffset;
            gECATszra[1]:=gDICOMdata.CompressSz;
    end else} if (gECATJPEG_table_entries > kMaxECAT) then
    begin
      gImgOK := false;
      Showmessage('This ECAT file has too many slices (' + inttostr(gECATJPEG_table_entries) + ').');
    end
    else
    begin
      gECATslices := gECATJPEG_table_entries;
      for lS := 1 to gECATslices do
      begin
        //showmessage(inttostr(gECATJPEG_pos_table[lS]));
        gECATposra[lS] := gECATJPEG_pos_table[lS];
        gECATszra[lS] := gECATJPEG_size_table[lS];
      end;
    end;
    freemem(gECATJPEG_pos_table);
    freemem(gECATJPEG_size_table);
    gECATJPEG_table_entries := 0;
  end;
  gBlack := 1;
  gScale := 1;
  gPro := 0;
  gCustomPalette := 0;
  if red_table_size > 0 then
  begin
    //gCustomPalette := 0;
  end
  else
  begin
    if gDICOMdata.monochrome = 1 then
      gScheme := 0
    else
      gScheme := 1;
    LoadColorScheme('', gScheme); //load Black and white
  end;
  gWinCen := 0;
  gWinWid := 0;
  if (gDICOMdata.XYZdim[2] < 1) or (gDICOMdata.XYZdim[1] < 1) or (not lHdrOK) or (not gImgOK) then
  begin
    showmessage('LoadData: Error reading image.');
    ReleaseDICOMmemory;
    OptionsImgInfoItemClick(nil);
    exit;
  end;
  LowerSlice1.enabled := gDicomdata.XYZdim[3] > 1;
  HigherSlice1.enabled := gDicomdata.XYZdim[3] > 1;
  Mosaic1.enabled := gDicomdata.XYZdim[3] > 1;
  if self.WindowState <> wsMaximized then
  begin
    self.ClientHeight := gDICOMdata.XYZdim[2];
    self.ClientWidth := (gDICOMData.XYZdim[1]);
  end;
  if (gDICOMdata.RLERedSz > 0) or (gDICOMdata.SamplesPerPixel > 1) or (gCustomPalette > 0) then
  begin
    gDICOMdata.WindowCenter := 127;
    gDICOMdata.WindowWidth := 255;
    gImgMin := 0;
    gImgMax := 255;
    gWinCen := gDICOMdata.WindowCenter;
    gWinWid := 0; //gDICOMdata.WindowWidth;
    gFastCen := 128;
    gFastSlope := 128;
  end;
  gAbort := false;
  Overlay1.enabled := true;
  gSlice := 0; {force a new image to be displayed - so gSlice should be different from displayimage requested slice}
  DisplayImage(True, True, 1, -1, 0);
  Screen.Cursor := crDefault;
  if Self.Active then
    MainForm.ColUpdate;
  MainForm.StatusBar.Panels[5].text := 'id:s:a:i ' + inttostr(gDicomData.PatientIDInt) + ':' +
    inttostr(gDicomData.SeriesNum) + ':' + inttostr(gDicomData.AcquNum) + ':' + inttostr(gDicomData.ImageNum);

end;

procedure TMDIChild.FileOpenItemClick(Sender: TObject);
//File menu, OPEN
begin
  MainForm.FileOpenItemClick(Sender);
end;

procedure TMDIChild.FileExitItemClick(Sender: TObject);
begin
  MainForm.FileExitItemClick(Sender);
end;

procedure TMDIChild.HdrShow;
var
  lLen, lI: integer;
  lStr: string;
begin
  if not FDICOM then
  begin
    //showmessage('Unable to show DICOM header information. This is not a DICOM file.');
    EXIT;
  end;
  Memo1.Lines.Clear;
  //Memo1.lines.add(inttostr(gDicomData.ImageStart));
  lLen := Length(gDynStr);
  if lLen > 0 then
  begin
    lStr := '';
    for lI := 1 to lLen do
    begin
      if gDynStr[lI] <> kCR then
        lStr := lStr + gDynStr[lI]
      else
      begin
        Memo1.Lines.add(lStr);
        lStr := '';
      end;
    end;
    Memo1.Lines.Add(lStr);
  end; //lLen > 0
end;

procedure TMDIChild.OptionsImgInfoItemClick(Sender: TObject);
begin
  MainForm.HdrBtn.Down := not MainForm.HdrBtn.Down;
  MainForm.HdrBtn.Click;
end;
(*
procedure TMDIChild.decompressJPEG24x (lFilename: string; var lOutputBuff: ByteP0; lImageVoxels,lImageStart{gECATposra[lSlice]}: integer);
var
Stream: Tmemorystream;
Jpg: TJPEGImage;

TmpBmp: TPicture;
lImage: Timage;
lRow:  pRGBTripleArray;
lHt0,lWid0,lInc,i,j: integer;
begin
try
Stream := TMemoryStream.Create;
Stream.LoadFromFile(lFilename);
Stream.Seek(lImageStart, soFromBeginning);
try
Jpg := TJPEGImage.Create;
Jpg.LoadFromStream(Stream);
//lImage.Create(Image);

Image.Height := JPG.Height;
Image.Width := JPG.Width;

//Image.Picture.Graphic:=JPG;
//Image.Picture.Assign(jpg);
Image.Picture.Bitmap.Assign(jpg);
{Image.Picture.Bitmap.Height := JPG.Height;
Image.Picture.Bitmap.Width := JPG.Width;
Image.Picture.Bitmap.PixelFormat := pf24bit;
     {}
//lImageVoxels = (JPG.Height*JPG.Width);
lWid0 := JPG.Width-1;
lHt0 := JPG.Height-1;
lInc := (3*lImageVoxels)-1; //*3 because 24-bit, -1 since index is from 0
//showmessage(inttostr(lWid0)+'@'+inttostr(lHt0));
FOR j := lHt0