Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  838 1762 (N33E37 L=264)   421     44        2 ( 813/1733- 859/1791) -51.3 E39-E35 N30-N35
      2 1203 1609 (N32E06 L=295)  1121    305        1 (1067/1567-1296/1682)  47.3 E14-E01 N27-N41
      3 1164 1389 (N17E02 L=300)   209     66        1 (1134/1359-1187/1416)  78.7 E03-E01 N14-N19
      4 1516 1262 (N20W22 L=324)   908    274        0 (1458/1166-1611/1303)  57.5 W19-W29 N12-N23
      5 1870  934 (N14W61 L=  3)  1093    131        2 (1849/ 865-1890/1003) -82.9 W57-W63 N09-N18
      6 1041  932 (S12W03 L=305)   264    110        0 ( 998/ 898-1098/ 958)   0.0 W00-W07 S12-S10
      7  980  737 (S24W06 L=307)   817    324        1 ( 881/ 677-1093/ 833)  58.4 W00-W09 S29-S16
      8 1190  699 (S20W19 L=321)   301    111        2 (1140/ 663-1230/ 727)   0.0 W15-W21 S21-S19
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y