Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  869 1712 (N38E31 L=143)   439     89        1 ( 837/1691- 909/1741) -12.5 E35-E26 N36-N39
      2  920 1574 (N29E20 L=154)   643    182        1 ( 823/1534- 980/1613)  45.0 E25-E17 N24-N33
      3 1662 1279 (N28W38 L=212)  2759    758        2 (1453/ 983-1798/1550) -38.9 W18-W54 N12-N41
      4  593 1433 (N13E37 L=137)   683    138        1 ( 547/1386- 649/1491)  76.0 E39-E34 N09-N17
      5 1157 1377 (N22W00 L=174)   244    137        0 (1106/1329-1204/1429) -26.6 E04-W04 N19-N24
      6 1325 1356 (N24W12 L=186)   488    230        0 (1240/1321-1401/1424)   0.0 W04-W16 N22-N27
      7 1368 1292 (N21W15 L=189)   236     45        0 (1349/1274-1389/1316) -33.7 W14-W17 N20-N23
      8  264 1225 (S05E57 L=117)   491     56        1 ( 243/1201- 291/1254)  76.0 E59-E54 S07-S03
      9  502 1023 (S13E33 L=141)   709    178        1 ( 458/ 931- 566/1066)  66.0 E36-E28 S18-S09
     10  258  856 (S28E54 L=120)   828    255        0 ( 217/ 782- 292/ 960)  56.3 E60-E49 S33-S21
     11 1000  869 (S12W02 L=176)   412    119        0 ( 939/ 851-1054/ 884)   9.5 E01-W05 S12-S09
     12  723  802 (S21E14 L=160)   258     86        0 ( 684/ 778- 756/ 840)  51.3 E16-E12 S23-S18
     13 1633  586 (S14W53 L=227)  3054    491        8 (1542/ 434-1733/ 746) -84.8 W47-W59 S24-S02
     14 1534  669 (S11W41 L=215)   342     63        2 (1503/ 642-1563/ 696)  26.6 W38-W43 S12-S10
     15 1097  601 (S25W15 L=189)   252     82        1 (1056/ 577-1134/ 620)  45.0 W12-W16 S27-S23
     16 1503  558 (S19W43 L=217)   426    108        0 (1448/ 539-1547/ 585)  33.7 W40-W46 S20-S16
     17  439  535 (S47E37 L=137)   657    100        1 ( 402/ 504- 475/ 572)  28.6 E42-E31 S50-S43
  

    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