Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1228 1898 (N60W65 L= 31)   826    111        1 (1173/1877-1274/1920)  -1.9 W52-W82 N58-N61
      2  628 1700 (N59E25 L=301)   464     30        1 ( 599/1679- 661/1717)  -6.3 E29-E20 N57-N61
      3 1155 1672 (N45W26 L=353)   951    185        3 (1075/1633-1246/1697)  -3.6 W18-W34 N42-N46
      4  537 1671 (N58E37 L=289)   483     42        2 ( 513/1641- 558/1696)  76.0 E41-E33 N55-N60
      5 1300 1640 (N39W36 L=  2)   277     26        0 (1276/1621-1319/1659) -56.3 W34-W37 N37-N40
      6 1552 1615 (N31W55 L= 21)   693     61        1 (1530/1577-1585/1648)  38.7 W51-W59 N29-N33
      7 1368 1505 (N29W34 L=360)   290     84        0 (1327/1481-1406/1524) -45.0 W31-W35 N26-N30
      8  853 1387 (N31E05 L=321)  1907    533        4 ( 623/1313-1087/1479) -30.1 E22-W09 N22-N40
      9 1715 1442 (N17W59 L= 26)   375     68        1 (1680/1425-1746/1460)  -7.1 W55-W63 N15-N18
     10 1467 1266 (N13W32 L=358)   254     76        0 (1434/1233-1489/1300) -68.2 W30-W33 N10-N15
     11 1107 1085 (N09W06 L=332)   265     92        0 (1090/1036-1123/1129)  68.2 W04-W07 N06-N11
     12  658 1088 (N16E22 L=305)   206    107        0 ( 619/1058- 700/1112)  18.4 E24-E18 N14-N16
     13  432  994 (N13E39 L=287)  1070    299        4 ( 328/ 876- 528/1067) -45.0 E47-E33 N04-N18
     14 1416  859 (S10W21 L=347)   549     97        0 (1396/ 806-1434/ 920)  74.1 W19-W22 S13-S06
     15  907  614 (S16E14 L=312)   304    118        0 ( 865/ 585- 943/ 658)  23.2 E17-E10 S17-S13
     16 1078  442 (S31E07 L=319)  1285    325        2 ( 976/ 374-1188/ 575)  31.4 E16-W02 S34-S23
     17  671  221 (S43E53 L=274)   986    101        2 ( 621/ 186- 715/ 251)  10.3 E67-E45 S46-S39
  

    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