Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  601 1410 (N12E35 L=103)   316     86        0 ( 575/1381- 635/1445) -31.0 E37-E32 N10-N13
      2  484 1216 (S02E38 L=100)   617    204        2 ( 434/1126- 539/1279) -26.6 E43-E31 S05-N01
      3 1736 1146 (N21W45 L=183)   839    114        1 (1697/1091-1762/1198) -50.2 W41-W47 N17-N23
      4  152 1125 (S13E68 L= 70)  1105    146        2 ( 138/1044- 169/1189) -71.6 E70-E65 S17-S08
      5 1199 1151 (N09W08 L=146)   336    112        0 (1139/1117-1250/1174)   0.0 W03-W11 N07-N10
      6  463  926 (S19E34 L=104)   596    220        2 ( 380/ 866- 532/ 994)  47.7 E39-E29 S24-S13
      7  510  807 (S25E29 L=109)  1123    400        0 ( 392/ 717- 656/ 859)  36.4 E38-E18 S33-S19
      8 1324  752 (S11W24 L=162)   305    104        0 (1274/ 728-1376/ 776)  45.0 W21-W26 S13-S08
      9 1574  429 (S25W55 L=193)  1269    222        2 (1506/ 362-1656/ 520)  77.0 W52-W58 S30-S17
     10 1115  423 (S37W22 L=160)   663    165        0 (1030/ 409-1196/ 439)  24.4 W16-W27 S39-S34
  

    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