Le portail francophone de la géomatique

while(fgets(s,127,fd) != NULL )
{
  REMOVE_BL(s);
  if( ! *s )
     continue;

  if (!strncmp(s,"X0 (m)",6) || !strncmp(s,"X0(m)",5))
  {
     fscanf(fd,"%f\n",&pci_trf.translation.x);
     fgets(s,127,fd);
     fscanf(fd,"%f\n",&pci_trf.translation.y);
     fgets(s,127,fd);
     fscanf(fd,"%f\n",&pci_trf.rotation);
     fgets(s,127,fd);a
     fscanf(fd,"%f\n",&pci_trf.echelle);
    break;
  }
}

fclose(fd);

// Ouverture du fichier .tfw
mdlFile_buildName(geo_filename,unite,dir,dummy,"tfw");
fd = fopen(geo_filename,"w+");
if (fd)
{
  angle = pci_trf.rotation / 200. * fc_pi;
  fprintf(fd,"%.6f\n",pci_trf.echelle * cos(angle));
  fprintf(fd,"%.6f\n",pci_trf.echelle * sin(angle));
  fprintf(fd,"%.6f\n",pci_trf.echelle * sin(angle));
  fprintf(fd,"%.6f\n",-pci_trf.echelle * cos(angle));
  fprintf(fd,"%.6f\n",pci_trf.translation.x);
  fprintf(fd,"%.6f\n",pci_trf.translation.y);
  fclose(fd);
}

Paramètre de la transformation associée :
X0 (m)=
 837545.05                        translation en X (X0)
Y0 (m)=
 55733.855                        translation en Y (Y0)
Téta (Grad)=
 327.2881                          angle de rotation en grades (Téta)
K =
 .106764                            coefficient d’échelle (K)

0.044375                 K x cos(Téta)
-0.097105                          K x sin(Téta)
-0.097105                          K x sin(Téta)
-0.044375                          -1 x K x cos(Téta) 
837545.050000                 X0
55733.855000                   Y0