Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1016 1791 (N52W31 L= 17)   311     45        1 ( 994/1769-1040/1811) -90.0 W29-W33 N49-N54
      2  761 1697 (N54W04 L=349)  2486    481        5 ( 588/1599- 935/1788)   1.5 E16-W23 N50-N56
      3 1471 1726 (N31W59 L= 44)   639    105        0 (1425/1691-1512/1757) -80.5 W56-W60 N28-N34
      4 1404 1611 (N27W44 L= 30)   336     83        0 (1384/1568-1427/1647)  68.2 W41-W46 N24-N29
      5  700 1464 (N39E09 L=337)   343     55        2 ( 670/1440- 729/1486) -45.0 E11-E06 N36-N40
      6 1494 1437 (N15W40 L= 26)   581    122        0 (1434/1412-1550/1471) -54.5 W37-W43 N11-N18
      7  492 1125 (N23E31 L=315)  1190    302        1 ( 355/1068- 614/1179) -27.8 E41-E21 N17-N29
      8  221 1061 (N25E56 L=289)   735    135        0 ( 203/ 994- 252/1122)  60.3 E59-E52 N21-N28
      9 1541 1010 (S10W30 L= 15)   722    197        1 (1489/ 930-1582/1108)  35.0 W24-W34 S13-S05
     10  365  980 (N17E43 L=303)  1232    310        4 ( 274/ 886- 454/1087) -49.1 E50-E37 N10-N25
     11  646 1045 (N15E22 L=324)   219     43        1 ( 629/1019- 669/1067) -90.0 E22-E20 N13-N16
     12 1224  843 (S10W07 L=352)   927    268        3 (1185/ 738-1281/ 983)  78.7 W04-W10 S17-S02
     13 1641  806 (S24W33 L= 19)  1427    276        2 (1521/ 727-1715/ 900)  13.2 W22-W40 S26-S20
     14  305  808 (N09E52 L=293)   438    134        0 ( 284/ 750- 326/ 873) -90.0 E53-E51 N05-N12
     15 1422  796 (S18W17 L=  3)   389    129        0 (1378/ 753-1461/ 853)  20.6 W13-W21 S19-S16
     16  823  396 (S27E31 L=314)   495    128        0 ( 778/ 347- 876/ 434)  15.3 E37-E26 S29-S25
     17  631  280 (S30E54 L=292)  1554    221        1 ( 588/ 207- 654/ 384)  29.7 E67-E45 S35-S23
  

    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