Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  635 1741 (N57E39 L=274)   446     58        1 ( 604/1721- 669/1754)   6.3 E43-E34 N54-N57
      2  834 1524 (N38E11 L=301)   551    153        1 ( 768/1484- 906/1557)  20.0 E16-E05 N35-N39
      3 1660 1418 (N24W51 L=  3)   751    179        0 (1613/1367-1721/1482) -58.0 W47-W55 N20-N28
      4 1758 1376 (N20W60 L= 12)   969    154        0 (1723/1303-1809/1433) -69.4 W57-W62 N15-N23
      5 1799 1386 (N20W65 L= 18)   589     95        0 (1770/1346-1831/1412) -38.7 W62-W69 N17-N22
      6  236 1230 (N19E60 L=252)  1453    194        5 ( 200/1147- 274/1319)  70.0 E62-E57 N14-N25
      7  931 1135 (N12E05 L=307)   413    173        0 ( 870/1094-1003/1159) -20.6 E09-W00 N09-N13
      8  812 1094 (N10E13 L=300)   359    151        1 ( 781/1032- 868/1127)  51.3 E14-E09 N06-N12
      9  347 1056 (N09E46 L=266)   522    156        0 ( 317/ 998- 396/1087)  45.0 E48-E41 N05-N10
     10  617  889 (S02E26 L=286)  1794    433        1 ( 557/ 763- 748/1070) -58.6 E29-E18 S09-N09
     11 1479  935 (S04W28 L=341)   704    238        0 (1380/ 884-1556/1009)  17.1 W21-W34 S06-N00
     12  290  739 (S11E55 L=258)   553    131        2 ( 263/ 702- 342/ 807) -39.8 E57-E49 S13-S06
     13  266  651 (S17E61 L=251)  2150    184        4 ( 220/ 589- 310/ 736) -60.9 E66-E57 S20-S11
     14  902  559 (S24E11 L=301)   302    106        0 ( 850/ 531- 943/ 586) -33.7 E14-E08 S26-S21
     15 1031  342 (S41E04 L=308)  2658    356        7 ( 926/ 247-1141/ 466)  43.4 E14-W04 S50-S32
     16  600  261 (S47E50 L=262)  1608    146        3 ( 523/ 239- 671/ 288)   7.9 E69-E39 S49-S44
  

    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