Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  409 1719 (N20E78 L=113)  1304    103        0 ( 370/1686- 445/1752)  71.6 E85-W00 N17-N23
      2 1244 1689 (N42E08 L=183)   418     67        1 (1217/1651-1268/1719) -71.6 E10-E06 N39-N45
      3  388  800 (S35E37 L=154)   285     47        1 ( 361/ 780- 411/ 826)  11.3 E40-E35 S36-S33
      4  237  612 (S53E65 L=126)  1061     89        2 ( 186/ 598- 282/ 630)   3.7 E85-E54 S53-S50
      5  964  452 (S40W16 L=208)   983    284        3 ( 934/ 329- 999/ 584) -41.4 W09-W26 S47-S32
      6 1274  503 (S27W33 L=224)   270    120        0 (1212/ 485-1327/ 520)  39.8 W29-W35 S29-S24
      7 1423  398 (S27W50 L=241)   665    138        1 (1357/ 376-1485/ 425)  45.0 W46-W53 S30-S23
  

    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