Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  751 1856 (N72E36 L=151)   835     82        1 ( 731/1820- 777/1884) -26.6 E52-E27 N66-N76
      2 1785 1459 (N23W71 L=259)   602     49        1 (1768/1436-1806/1482) -56.3 W67-W74 N21-N24
      3  636 1463 (N36E26 L=162)   303     71        1 ( 597/1446- 671/1480)  -9.5 E29-E23 N34-N36
      4 1135 1469 (N32W12 L=200)   335    120        0 (1072/1461-1194/1480)   0.0 W07-W16 N31-N33
      5 1464 1428 (N26W35 L=223)   322    161        0 (1429/1397-1534/1476) -39.8 W33-W40 N24-N29
      6 1034 1365 (N26W03 L=191)   586    166        2 (1005/1289-1066/1444) -84.3 W01-W05 N21-N31
      7  334 1370 (N30E53 L=134)   423     80        1 ( 320/1324- 350/1412) -71.6 E55-E51 N27-N33
      8 1427 1370 (N23W31 L=219)   444    168        1 (1358/1312-1498/1398) -41.2 W26-W35 N18-N25
      9 1746 1264 (N12W55 L=243)  1336    291        5 (1718/1170-1773/1397)  54.5 W51-W62 N06-N20
     10  688 1161 (N16E21 L=167)   667    204        3 ( 631/1097- 743/1210) -35.5 E24-E17 N11-N19
     11 1903  989 (S07W70 L=258)   597     49        1 (1876/ 966-1928/1009)   5.7 W65-W75 S07-S05
     12 1746  854 (S13W50 L=238)   810     71        1 (1701/ 825-1778/ 887)   7.1 W45-W53 S14-S10
     13 1705  703 (S22W48 L=235)   483    126        0 (1653/ 669-1743/ 747)  29.7 W43-W51 S23-S19
     14  231  661 (S15E67 L=120)  1973    150        4 ( 200/ 597- 268/ 743)  37.6 E76-E61 S19-S09
     15 1520  632 (S24W32 L=219)  1519    304        1 (1438/ 508-1575/ 743)  59.0 W26-W35 S32-S17
     16 1753  661 (S26W54 L=242)  1004    169        0 (1693/ 599-1809/ 722)  35.0 W48-W60 S29-S22
     17  888  607 (S20E12 L=175)   309    103        0 ( 845/ 579- 932/ 632)  18.4 E15-E09 S21-S18
     18  954  449 (S32E10 L=178)   316    112        1 ( 935/ 400- 971/ 506)  63.4 E12-E08 S35-S27
     19  860  350 (S40E20 L=168)   278     71        0 ( 828/ 333- 902/ 362)   9.5 E22-E16 S40-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