Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1072 1695 (N34E16 L=307)  1175    191        3 ( 984/1637-1156/1737)  -8.1 E22-E08 N32-N36
      2 1330 1415 (N23W08 L=331)   273     69        0 (1311/1378-1347/1454) -78.7 W07-W09 N20-N25
      3 1183 1408 (N18E00 L=322)   366     98        1 (1144/1383-1238/1435)   0.0 E02-W03 N16-N19
      4  290 1391 (S01E58 L=264)   904    124        1 ( 238/1374- 360/1419)   0.0 E65-E51 S01-N00
      5 1838 1229 (N27W52 L= 14)   513     82        0 (1822/1181-1865/1260)  54.5 W49-W55 N23-N30
      6 1648 1109 (N13W35 L=358)   560    145        1 (1624/1051-1686/1172) -39.8 W33-W39 N10-N16
      7 1860 1067 (N18W55 L= 18)  2098    236        5 (1824/1003-1943/1154) -13.2 W50-W69 N14-N22
      8 1373 1071 (N04W18 L=341)   540    164        0 (1288/1044-1452/1107)  41.2 W13-W22 N00-N07
      9  497  893 (S26E30 L=292)   269     32        0 ( 471/ 873- 524/ 913) -11.3 E32-E27 S27-S24
     10 1800  765 (S00W58 L= 20)   681     98        0 (1761/ 749-1856/ 790)   5.2 W53-W64 S01-N01
     11 1719  642 (S09W56 L= 18)   593     83        3 (1691/ 614-1752/ 685) -14.0 W52-W60 S10-S07
     12  284  644 (S46E53 L=269)  2451    313        1 ( 153/ 605- 432/ 679)   1.1 E85-E34 S47-S44
     13 1086  634 (S27W14 L=337)   198     84        0 (1048/ 609-1122/ 660)  -8.1 W10-W17 S28-S26
     14 1702  547 (S14W60 L= 23)   935    146        4 (1671/ 487-1732/ 602) -60.3 W56-W64 S18-S11
     15 1637  503 (S19W57 L= 20)   464     96        1 (1617/ 449-1654/ 538) -51.3 W53-W60 S21-S16
     16 1461  487 (S25W44 L=  7)   681    110        1 (1406/ 463-1510/ 512)   6.3 W39-W48 S26-S22
     17 1466  334 (S33W57 L= 19)  2026    290        3 (1330/ 304-1598/ 374)  33.7 W47-W66 S38-S26
  

    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