Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1502 1602 (N30W43 L= 59)   554     94        0 (1455/1567-1543/1626) -38.7 W39-W46 N27-N32
      2 1051 1456 (N25W07 L= 22)   726    152        0 ( 977/1429-1126/1478) -16.7 W01-W11 N22-N27
      3 1755 1450 (N18W61 L= 76)   477     44        0 (1727/1432-1785/1465)   7.1 W57-W65 N16-N18
      4  749 1333 (N21E14 L=  1)   221     57        0 ( 728/1308- 774/1364) -45.0 E15-E12 N18-N22
      5  333  610 (S18E56 L=320)   594     79        1 ( 305/ 586- 370/ 638) -18.4 E59-E52 S19-S16
      6 1687  415 (S49W61 L= 77)   728     49        0 (1659/ 392-1716/ 440) -11.9 W53-W72 S51-S46
      7 1628  389 (S50W53 L= 69)   600     71        0 (1592/ 369-1660/ 406) -12.1 W46-W61 S52-S48
      8 1499  351 (S52W37 L= 52)   706     90        0 (1474/ 313-1533/ 394) -33.7 W32-W44 S55-S47
      9 1548  342 (S53W45 L= 60)   654     67        0 (1522/ 309-1573/ 374) -26.6 W38-W53 S57-S50
  

    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