Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  856 1763 (N45E32 L=210)   641     73        1 ( 833/1727- 888/1796) -40.6 E36-E29 N41-N47
      2  280 1516 (N17E70 L=171)   861     84        1 ( 250/1485- 316/1540)  -6.7 E79-E62 N15-N18
      3 1743 1489 (N42W55 L=297)   614     82        1 (1702/1465-1788/1509)  -4.8 W48-W62 N41-N43
      4 1048 1218 (N11E02 L=240)  1285    441        3 ( 970/1125-1152/1287)  48.4 E06-W03 N05-N16
      5  797 1172 (N04E16 L=226)   317    168        0 ( 786/1082- 808/1253) -71.6 E17-E14 N00-N09
      6 1162 1204 (N12W05 L=247)   340    166        2 (1087/1181-1239/1229)  11.3 W00-W10 N10-N13
      7  229 1050 (S12E57 L=185)   573    139        2 ( 207/ 988- 247/1127)  90.0 E59-E55 S15-S07
      8 1052  846 (S11W05 L=246)   771    357        1 ( 981/ 716-1137/ 995)  75.3 W01-W07 S19-N00
      9  171  899 (S22E65 L=177)  1171    152        2 ( 115/ 847- 232/ 943)  20.0 E78-E55 S26-S18
     10 1487  658 (S15W37 L=278)   395     98        1 (1442/ 619-1524/ 689) -20.6 W32-W40 S16-S13
     11  526  361 (S56E35 L=207)  1553    168        2 ( 465/ 317- 600/ 406) -14.9 E43-E25 S58-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