Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1141 1830 (N49E20 L=153)   496     82        0 (1095/1807-1179/1849)   0.0 E24-E16 N47-N50
      2  623 1688 (N23E46 L=127)   852    172        0 ( 564/1632- 654/1760) -78.7 E50-E43 N18-N28
      3 1343 1354 (N20W10 L=183)   207     76        0 (1315/1324-1374/1385) -26.6 W07-W12 N18-N20
      4  926 1332 (N08E13 L=160)   245     36        1 ( 901/1313- 948/1350)  71.6 E14-E11 N06-N09
      5 1271 1209 (N09W09 L=182)   375    120        2 (1250/1154-1291/1272) -56.3 W07-W11 N06-N12
      6 1758 1001 (N12W47 L=220)   379     66        1 (1737/ 964-1771/1038) -63.4 W45-W48 N09-N13
      7  185  931 (S30E62 L=111)  1114    144        2 ( 128/ 894- 253/ 956)  11.9 E72-E53 S32-S27
      8  537  762 (S34E25 L=147)   373    130        0 ( 495/ 723- 574/ 822)  54.5 E28-E23 S36-S29
      9  244  794 (S38E55 L=117)   550     86        1 ( 203/ 774- 291/ 814)  13.0 E62-E49 S39-S36
     10  456  763 (S36E32 L=140)   273     74        1 ( 441/ 726- 477/ 801) -53.1 E33-E30 S37-S33
     11  804  375 (S53W06 L=179)   322     56        0 ( 777/ 354- 830/ 399)  -9.5 W02-W08 S54-S51
  

    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