Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  448 1634 (N18E58 L=289)   558     95        2 ( 429/1603- 468/1667) -38.7 E61-E55 N16-N20
      2  296 1564 (N11E72 L=276)  1910    187        2 ( 258/1525- 335/1604)  56.3 E76-E67 N08-N14
      3  224 1146 (S15E58 L=289)   538     99        2 ( 204/1116- 248/1180) -36.9 E60-E55 S16-S13
      4 1848 1048 (N20W56 L= 43)  1153    280        2 (1782/1025-1921/1076)  10.0 W48-W65 N17-N22
      5 1715  866 (N06W48 L= 35)   400     68        1 (1689/ 826-1738/ 903)  80.5 W46-W49 N02-N08
      6  945  638 (S28W06 L=353)   304    100        0 ( 891/ 607- 986/ 674)  -7.1 W01-W09 S29-S26
      7 1738  565 (S10W64 L= 52)  1178    206        4 (1707/ 486-1776/ 621) -77.5 W60-W68 S14-S05
      8 1639  524 (S15W56 L= 44)   627     85        3 (1614/ 489-1667/ 557) -68.2 W53-W58 S17-S12
      9 1118  374 (S40W28 L= 16)  1830    400        2 ( 986/ 317-1277/ 478)  -2.1 W14-W41 S43-S36
     10 1604  357 (S25W64 L= 52)   948    224        1 (1559/ 341-1667/ 370)  21.8 W59-W70 S26-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