Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1153 1779 (N50E13 L= 51)  4195    395       12 (1070/1609-1239/1923) -39.5 E36-E01 N37-N66
      2  389 1541 (N15E58 L=  6)  1991    237        6 ( 352/1443- 434/1647) -42.9 E68-E51 N08-N21
      3  616 1218 (N00E28 L= 36)  1170    242        2 ( 553/1122- 695/1315)  67.4 E30-E25 S05-N07
      4  934 1101 (N00E07 L= 57)   239    110        0 ( 921/1047- 945/1152) -71.6 E08-E05 S03-N03
      5 1184 1008 (N00W09 L= 74)   219    116        0 (1174/ 949-1197/1057) -80.5 W08-W11 S03-N03
      6  832 1001 (S08E11 L= 54)   322    128        0 ( 780/ 953- 881/1032)  56.3 E12-E08 S11-S05
      7 1785  774 (N00W58 L=122)  3552    306        3 (1721/ 640-1842/ 918) -28.6 W48-W72 S05-N07
      8 1445  859 (S03W29 L= 93)   562    111        1 (1402/ 817-1487/ 897)  68.2 W25-W30 S05-N00
      9  437  692 (S35E34 L= 30)   307     92        0 ( 395/ 658- 462/ 712) -18.4 E38-E31 S36-S33
     10 1178  171 (S54W52 L=116)   474     68        0 (1154/ 142-1199/ 201) -28.6 W47-W58 S56-S50
  

    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