Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1442 1780 (N55W57 L=165)  1261    115        1 (1409/1735-1495/1819)   9.0 W46-W84 N51-N57
      2  529 1680 (N49E51 L= 58)   470     66        1 ( 502/1652- 555/1710)  90.0 E52-E48 N46-N51
      3 1472 1580 (N38W40 L=148)   311     35        0 (1450/1559-1500/1594)   9.5 W37-W43 N37-N39
      4  740 1436 (N31E20 L= 88)  2311    564        1 ( 536/1290- 939/1566) -26.6 E40-E04 N20-N40
      5 1284 1270 (N19W18 L=126)   478    194        0 (1250/1197-1350/1331)  53.1 W15-W22 N14-N22
      6  989 1007 (N03E02 L=106)   468    130        3 ( 942/ 966-1046/1050)  33.7 E05-W01 N00-N05
      7 1691  931 (S04W45 L=153)   270     57        1 (1664/ 912-1712/ 950) -26.6 W43-W47 S05-S03
      8  864  697 (S16E11 L= 97)  1768    486        2 ( 726/ 491- 938/ 870)  55.8 E23-E05 S30-S05
      9 1490  726 (S16W30 L=139)   289    103        0 (1463/ 685-1525/ 775)  59.0 W28-W32 S18-S13
     10 1583  525 (S31W42 L=150)  4675    656        4 (1334/ 364-1771/ 718)  35.7 W23-W56 S42-S18
     11 1136  446 (S34W07 L=115)  1733    373        1 (1002/ 350-1268/ 542)  37.9 E03-W15 S42-S27
     12  676  467 (S32E28 L= 80)   276     47        0 ( 661/ 439- 691/ 493) -90.0 E29-E26 S34-S30
  

    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