Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1311 1769 (N39W51 L= 25)   936    145        0 (1232/1751-1384/1782) -21.0 W44-W57 N36-N42
      2  632 1707 (N61E08 L=326)   310     23        0 ( 614/1687- 648/1729)  63.4 E10-E05 N59-N63
      3 1657 1653 (N20W70 L= 43)  1086     86        2 (1615/1616-1692/1683) -80.5 W67-W74 N17-N23
      4 1313 1628 (N31W39 L= 13)   442    109        0 (1266/1592-1363/1662)   5.2 W34-W45 N30-N32
      5 1359 1537 (N25W37 L= 11)   707    190        0 (1267/1516-1451/1555) -21.0 W30-W43 N21-N27
      6  848 1524 (N40W04 L=338)   449    113        0 ( 811/1497- 888/1543)  -9.5 W01-W07 N37-N41
      7 1824 1472 (N06W74 L= 48)  1703    123        2 (1786/1421-1874/1512) -41.2 W68-W81 N01-N08
      8  412 1382 (N43E35 L=299)   709    115        0 ( 363/1346- 464/1423) -32.5 E40-E29 N39-N46
      9 1060 1287 (N20W09 L=343)   479    173        1 (1004/1202-1089/1343)  83.7 W07-W10 N14-N23
     10  968 1317 (N24W05 L=339)   198     48        0 ( 949/1287- 991/1342)  33.7 W03-W06 N22-N24
     11 1805 1256 (S04W58 L= 31)   767     94        1 (1776/1202-1841/1297)  18.4 W52-W62 S06-S02
     12  134 1032 (N26E69 L=264)  1018    138        0 ( 106/ 981- 178/1094) -31.6 E75-E62 N22-N30
     13 1664  883 (S20W35 L=  9)   299     81        0 (1627/ 867-1705/ 903) -33.7 W32-W38 S21-S17
     14 1313  705 (S19W08 L=342)   646    207        1 (1256/ 606-1349/ 786)  48.4 W02-W11 S23-S14
     15  752  635 (S09E26 L=308)   748    198        1 ( 724/ 543- 793/ 728)  42.0 E30-E20 S13-S03
     16  849  271 (S34E37 L=296)  3093    493        2 ( 654/ 163-1029/ 405)   8.9 E69-E18 S38-S29
  

    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