Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  541 1705 (N46E49 L=129)  3514    478        7 ( 415/1631- 679/1762)  20.9 E58-E35 N39-N51
      2 1132 1647 (N45W06 L=184)   292     62        1 (1108/1616-1160/1667)   0.0 W04-W09 N42-N46
      3  838 1629 (N42E18 L=160)   431     58        0 ( 806/1604- 869/1658)  45.0 E20-E16 N40-N44
      4 1241 1547 (N38W15 L=193)   295     51        1 (1219/1531-1271/1572) -45.0 W12-W17 N36-N39
      5  434 1360 (N21E43 L=135)   446     91        1 ( 406/1332- 462/1387) -26.6 E46-E40 N19-N22
      6  370 1253 (N14E47 L=131)  1169    277        3 ( 292/1191- 424/1290) -24.4 E54-E42 N10-N16
      7  363 1165 (N08E46 L=132)   986    309        3 ( 293/1110- 455/1226) -18.4 E53-E38 N05-N11
      8  202 1156 (N07E63 L=115)   545    128        1 ( 183/1104- 220/1211) -63.4 E66-E60 N03-N09
      9  626  941 (S04E25 L=153)   261     44        3 ( 607/ 895- 649/ 985)  78.7 E26-E23 S06-S01
     10 1866  845 (S07W66 L=244)   817    113        2 (1843/ 803-1887/ 886) -45.0 W62-W69 S09-S04
     11  391  776 (S15E43 L=135)  1623    315        3 ( 352/ 677- 438/ 868)  85.6 E45-E39 S22-S09
     12  972  775 (S13E02 L=176)   225     73        1 ( 950/ 735-1002/ 812)  68.2 E03-W00 S15-S10
     13 1429  580 (S24W30 L=208)   291     84        0 (1393/ 567-1471/ 599)  -8.1 W26-W33 S25-S23
     14 1233  430 (S36W18 L=196)   825    152        0 (1160/ 404-1302/ 460) -14.0 W12-W24 S37-S33
     15  834  401 (S40E13 L=165)  1296    309        3 ( 693/ 359- 969/ 437)  11.8 E25-E01 S44-S36
     16 1438  374 (S40W38 L=216)   350     57        0 (1414/ 344-1461/ 397)  76.0 W36-W39 S42-S38
  

    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