Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1097 1471 (N22E01 L= 49)   250     17        0 (1076/1452-1120/1494) -45.0 E02-W00 N21-N23
      2 1764 1464 (N33W50 L=100)   483     65        1 (1743/1432-1780/1491)  76.0 W47-W51 N31-N35
      3  149 1339 (N07E74 L=336)   843    121        1 ( 120/1295- 172/1385)  63.4 E79-E70 N04-N10
      4 1086 1229 (N07W01 L= 51)   181     34        0 (1068/1206-1108/1246)  56.3 W00-W02 N05-N08
      5  213 1143 (S05E59 L=351)   717    155        1 ( 175/1079- 248/1200) -41.2 E64-E54 S08-S01
      6  903  971 (S10E07 L= 43)   184     96        1 ( 880/ 945- 929/1003)  71.6 E07-E04 S11-S08
      7  753  847 (S19E15 L= 35)   197     49        0 ( 736/ 822- 778/ 877) -56.3 E15-E13 S20-S17
      8 1802  798 (S08W58 L=107)  1535    255        2 (1727/ 760-1871/ 846)   3.4 W49-W66 S09-S04
      9  571  758 (S26E27 L= 23)   830    247        0 ( 525/ 709- 665/ 825) -24.0 E30-E20 S29-S21
     10 1627  540 (S26W50 L=100)   604     77        0 (1595/ 499-1655/ 587) -80.5 W48-W52 S28-S22
  

    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