Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1332 1795 (N46W55 L=197)   472     51        1 (1312/1767-1355/1826)  38.7 W51-W60 N43-N47
      2 1473 1627 (N31W50 L=191)   370     25        0 (1454/1598-1490/1648)  21.8 W47-W52 N29-N32
      3 1570 1611 (N27W57 L=199)   573     77        1 (1530/1590-1611/1633) -45.0 W54-W60 N24-N29
      4  915 1482 (N37W04 L=145)   364     66        2 ( 879/1452- 947/1510)   0.0 W00-W06 N35-N38
      5  692 1422 (N38E14 L=127)   611    131        1 ( 672/1366- 717/1489)  83.7 E15-E12 N33-N42
      6  659 1176 (N23E20 L=121)   993    318        1 ( 606/1059- 731/1319) -83.7 E22-E16 N14-N32
      7 1674 1158 (S00W44 L=185)   338     90        0 (1661/1110-1688/1205)  36.9 W42-W46 S02-N02
      8 1874 1022 (S14W61 L=202)  1180    211        1 (1812/ 917-1908/1095)  29.7 W51-W67 S18-S10
      9 1899  860 (S25W64 L=205)   680     69        1 (1876/ 831-1922/ 893)  45.0 W60-W69 S26-S23
     10  717  781 (S01E23 L=118)   429    153        0 ( 644/ 763- 777/ 813)   0.0 E27-E18 S02-N00
     11 1504  686 (S24W23 L=164)  1203    287        2 (1426/ 558-1568/ 791)  45.0 W15-W29 S30-S18
     12  999  706 (S11E08 L=133)   203    102        0 ( 987/ 657-1010/ 757)  63.4 E09-E06 S13-S07
     13  384  633 (S05E51 L= 90)   398    125        0 ( 362/ 576- 410/ 685) -90.0 E52-E49 S08-S01
     14 1008  611 (S17E10 L=132)   534     95        3 ( 973/ 579-1049/ 642)  21.8 E11-E06 S18-S14
     15  868  501 (S21E21 L=120)   775    307        0 ( 801/ 424- 924/ 566) -84.3 E25-E17 S26-S16
     16  770  476 (S21E28 L=113)   467     91        1 ( 731/ 446- 803/ 515)  18.4 E31-E25 S22-S18
     17  847  296 (S36E32 L=109)  5139    692       12 ( 620/ 168-1106/ 496)  13.7 E80-E06 S44-S26
     18 1531  327 (S55W27 L=168)  1546    153        1 (1456/ 301-1602/ 364)  -9.0 W17-W37 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