Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1422 1837 (N40W69 L=283)   733     78        1 (1388/1815-1460/1853)  -7.6 W63-W79 N38-N41
      2 1417 1675 (N31W48 L=262)   733    203        0 (1323/1648-1509/1694) -30.3 W41-W54 N26-N34
      3 1075 1521 (N31W17 L=231)   594    241        0 ( 963/1498-1174/1548) -24.8 W09-W22 N27-N33
      4  750 1481 (N37E06 L=208)   956    297        0 ( 619/1434- 889/1529)  -5.7 E16-W04 N35-N38
      5 1460 1341 (N11W34 L=248)   716    214        0 (1405/1246-1505/1430) -90.0 W32-W35 N04-N16
      6 1889 1163 (S11W66 L=280)  1103    140        2 (1853/1121-1951/1217) -25.0 W62-W78 S14-S06
      7  369  881 (N10E45 L=169)  1828    517        1 ( 214/ 698- 537/ 998) -46.4 E58-E36 S03-N18
      8 1368  742 (S21W14 L=228)   530    130        0 (1329/ 684-1394/ 806)  50.2 W10-W16 S23-S17
      9 1737  638 (S39W43 L=257)  2639    376        4 (1635/ 519-1836/ 794)  31.8 W33-W54 S44-S30
     10  512  398 (S22E53 L=161)  1094    212        2 ( 476/ 308- 552/ 499)  45.0 E59-E45 S28-S15
     11  667  409 (S24E40 L=174)   478    160        0 ( 649/ 334- 693/ 480)  42.0 E45-E35 S29-S20
     12 1223  279 (S50E09 L=205)  1082    172        2 (1146/ 239-1310/ 327)  -3.6 E16-W00 S51-S47
  

    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