Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  839 1824 (N62W20 L=117)   530    103        0 ( 804/1788- 884/1867)  13.2 W13-W30 N59-N64
      2 1145 1788 (N46W40 L=136)   454     68        1 (1121/1747-1166/1818)  35.5 W36-W44 N43-N48
      3 1380 1733 (N35W51 L=148)  1141    114        3 (1311/1701-1434/1759) -24.0 W45-W55 N31-N37
      4 1467 1646 (N27W51 L=147)   581     95        1 (1440/1605-1492/1689)  63.4 W47-W53 N24-N30
      5  778 1599 (N45W01 L= 97)   234     64        1 ( 756/1568- 799/1627) -90.0 W00-W01 N42-N47
      6  527 1394 (N38E24 L= 73)   435    103        0 ( 488/1356- 579/1419) -40.6 E27-E20 N34-N41
      7  963 1364 (N25W05 L=102)   288    100        2 ( 935/1316- 992/1409) -90.0 W04-W06 N21-N28
      8  252 1384 (N45E56 L= 41)   655     65        1 ( 217/1360- 280/1401) -14.0 E63-E51 N42-N46
      9  223 1053 (N24E56 L= 41)  3268    584        5 ( 173/ 841- 291/1280) -79.5 E63-E49 N11-N38
     10  483 1228 (N29E30 L= 67)   225     84        1 ( 461/1206- 520/1250)   0.0 E31-E26 N27-N30
     11  381 1113 (N24E40 L= 57)   345    129        3 ( 330/1087- 448/1138)  -5.2 E44-E33 N23-N26
     12  437  946 (N14E38 L= 59)   243     86        0 ( 410/ 918- 458/ 988)  31.0 E40-E35 N12-N15
     13 1139  834 (S09W02 L= 99)   476    166        0 (1113/ 755-1175/ 898)  45.0 E01-W05 S12-S05
     14 1756  647 (S39W44 L=141)   722    101        2 (1708/ 621-1808/ 667)  -4.4 W38-W51 S40-S37
     15 1192  540 (S28E03 L= 94)   215     36        0 (1164/ 526-1213/ 558)   0.0 E04-E01 S28-S26
     16  911  468 (S25E22 L= 75)   371    116        0 ( 845/ 453- 965/ 484) -14.0 E26-E18 S26-S23
     17 1201  220 (S54E17 L= 80)  1543    245        2 (1101/ 150-1296/ 274)   8.7 E32-E06 S58-S49
  

    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