Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  449 1669 (N36E54 L=357)  2227    201        3 ( 392/1573- 499/1721) -24.8 E70-E44 N28-N40
      2  916 1682 (N38E13 L= 38)   411    103        1 ( 875/1662- 967/1714) -14.0 E16-E08 N36-N40
      3  626 1392 (N16E28 L= 23)   664    162        0 ( 591/1325- 671/1464)  71.6 E29-E25 N12-N21
      4 1650 1427 (N25W43 L= 94)   603     89        0 (1602/1390-1690/1457) -23.2 W38-W46 N22-N26
      5 1355 1399 (N21W19 L= 70)   278     87        0 (1311/1371-1394/1431) -31.0 W16-W21 N18-N22
      6  788 1357 (N15E16 L= 35)   406    139        0 ( 725/1320- 842/1379) -14.0 E20-E12 N12-N16
      7  736 1268 (N09E19 L= 32)   227     96        0 ( 699/1229- 759/1292) -36.9 E21-E17 N06-N10
      8  187 1210 (N05E63 L=348)   957    128        4 ( 168/1154- 209/1265)  90.0 E65-E60 N01-N08
      9  595 1179 (N03E27 L= 24)   318     93        0 ( 573/1138- 624/1227) -63.4 E28-E25 N00-N06
     10  195 1094 (S02E60 L=351)  1096    176        1 ( 154/1034- 224/1186)  48.8 E65-E57 S05-N03
     11  309  847 (S17E50 L=  1)   601    122        2 ( 266/ 796- 345/ 893)  45.0 E54-E46 S20-S14
     12  553  864 (S16E30 L= 21)   599    218        0 ( 452/ 838- 650/ 885)   7.6 E37-E22 S17-S14
     13 1386  809 (S15W24 L= 75)   264     69        0 (1370/ 776-1406/ 844) -63.4 W22-W25 S16-S12
     14 1497  744 (S18W33 L= 84)   341     68        0 (1463/ 725-1529/ 767)  36.9 W30-W34 S19-S16
     15 1533  676 (S22W37 L= 88)   994    282        2 (1465/ 576-1617/ 756)  63.4 W32-W42 S28-S16
     16 1282  706 (S22W19 L= 70)   360    142        0 (1222/ 683-1352/ 741) -12.5 W14-W23 S23-S20
     17 1439  571 (S29W33 L= 84)   314     76        0 (1418/ 550-1469/ 605) -45.0 W30-W35 S30-S27
     18 1587  480 (S34W50 L=101)   837     99        3 (1551/ 453-1627/ 524) -12.5 W45-W54 S35-S30
     19 1221  421 (S42W20 L= 71)   311     54        0 (1197/ 402-1255/ 447)   0.0 W17-W23 S43-S39
  

    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