Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  887 1233 (N01E13 L=164)   259    118        1 ( 836/1219- 947/1248)  26.6 E16-E10 N00-N03
      2  195 1056 (S24E61 L=116)   415    113        0 ( 151/1017- 235/1096)  28.6 E67-E56 S27-S21
      3  361  856 (S33E41 L=136)   573    162        0 ( 332/ 781- 405/ 923)  63.4 E43-E37 S38-S28
      4  257  854 (S36E52 L=124)   380    113        0 ( 247/ 791- 271/ 901)  74.1 E53-E50 S39-S32
      5  470  836 (S32E31 L=146)   262     86        1 ( 443/ 807- 506/ 875) -31.0 E33-E27 S32-S29
      6  561  788 (S32E23 L=154)   467    120        0 ( 506/ 770- 609/ 819)  40.6 E26-E19 S34-S28
      7 1143  670 (S22W18 L=194)   695    294        0 (1108/ 589-1217/ 750) -79.7 W14-W23 S26-S15
      8  799  177 (S71W35 L=211)   683    106        1 ( 751/ 149- 846/ 211)   0.0 W22-W47 S72-S67
  

    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