Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  863 1794 (N40E32 L=329)   901    117        2 ( 803/1766- 923/1820)  -4.8 E38-E26 N37-N41
      2 1211 1695 (N40E02 L=359)  1024    153        2 (1160/1636-1251/1752) -41.2 E07-W01 N35-N43
      3 1099 1716 (N39E11 L=350)   618    106        1 (1048/1672-1142/1750) -18.4 E15-E06 N35-N40
      4  362 1610 (N19E64 L=297)  1172    123        2 ( 328/1549- 405/1655)  90.0 E66-E61 N14-N22
      5  534 1422 (N09E38 L=322)   286     67        0 ( 505/1390- 558/1453)  76.0 E39-E37 N06-N10
      6 1121 1361 (N15E01 L=  0)   419    158        0 (1045/1334-1183/1375)   7.1 E05-W03 N13-N16
      7 1559 1333 (N23W28 L= 29)   251     57        0 (1536/1311-1588/1361) -14.0 W26-W30 N21-N24
      8  210 1137 (S12E60 L=301)  1761    288        3 ( 125/1024- 282/1233)  37.4 E71-E53 S18-S05
      9  264 1090 (S14E53 L=307)   362     70        0 ( 242/1058- 285/1121)  53.1 E55-E51 S16-S12
     10  328  980 (S20E46 L=314)   399    123        1 ( 315/ 921- 337/1048)  90.0 E47-E45 S23-S16
     11 1554  876 (S04W35 L= 36)   222     61        0 (1541/ 845-1568/ 912) -76.0 W34-W36 S06-S02
     12  279  764 (S34E53 L=308)   578    114        2 ( 236/ 728- 309/ 816)  40.6 E58-E48 S36-S30
     13 1563  715 (S13W41 L= 41)   263     93        1 (1536/ 686-1604/ 740)  59.0 W39-W42 S15-S10
     14 1323  612 (S24W28 L= 29)   863    247        4 (1290/ 498-1365/ 726) -63.4 W25-W32 S31-S17
     15 1697  676 (S12W53 L= 54)   544    122        1 (1659/ 634-1727/ 715) -12.5 W48-W57 S14-S09
     16 1555  479 (S26W51 L= 52)  1806    336        5 (1437/ 390-1682/ 575)  64.8 W46-W58 S34-S17
     17 1606  387 (S30W65 L= 65)  1148    156        0 (1530/ 345-1672/ 432)  52.1 W59-W68 S34-S25
  

    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