Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1493 1599 (N35W33 L=307)   328     28        1 (1472/1572-1515/1619)  53.1 W30-W35 N33-N37
      2 1339 1579 (N33W20 L=294)   331     62        0 (1325/1541-1353/1618) -90.0 W18-W20 N29-N35
      3  204 1425 (N19E68 L=206)  2236    191        4 ( 151/1353- 266/1503)  51.3 E75-E63 N14-N24
      4 1700 1449 (N26W48 L=323)  1166    184        2 (1629/1393-1775/1502) -21.0 W42-W55 N23-N29
      5  276 1347 (N14E55 L=219)   606    118        2 ( 229/1318- 332/1384)  11.3 E60-E50 N11-N15
      6  440 1275 (N09E39 L=235)   715    177        1 ( 362/1249- 521/1297)  -8.7 E45-E32 N07-N10
      7 1864 1258 (N16W62 L=336)   572    110        0 (1809/1225-1900/1283) -10.3 W55-W67 N14-N16
      8 1883 1108 (N07W62 L=337)  2003    339        5 (1813/ 961-1942/1197)  70.3 W54-W72 S02-N12
      9  174  980 (S08E62 L=212)   817     85        2 ( 127/ 949- 211/1003) -10.3 E68-E57 S10-S06
     10  282  925 (S12E51 L=224)   483     85        0 ( 259/ 881- 302/ 966)  51.3 E53-E48 S14-S09
     11 1876  765 (S13W68 L=342)   839    129        3 (1836/ 736-1922/ 805)  24.4 W63-W74 S15-S10
     12 1359  598 (S28W25 L=300)  2238    527        4 (1270/ 462-1529/ 723)  58.0 W19-W35 S37-S19
     13  433  577 (S34E44 L=231)   920    237        1 ( 343/ 518- 531/ 612)   9.9 E57-E34 S37-S31
     14 1490  535 (S32W37 L=312)   422     86        1 (1451/ 504-1533/ 574)  51.3 W34-W40 S33-S28
     15 1809  537 (S28W74 L=348)   922     56        0 (1777/ 521-1843/ 553)   3.4 W67-W85 S28-S26
     16  878  337 (S50E09 L=265)  1795    356        3 ( 685/ 316-1022/ 371)   5.4 E28-W04 S52-S47
     17 1584  342 (S44W61 L=335)  2091    191        4 (1499/ 314-1673/ 370)   0.0 W48-W83 S45-S41
  

    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