Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1021 1798 (N52W30 L= 21)   409     62        1 ( 981/1784-1054/1818) -38.7 W25-W32 N49-N53
      2 1570 1689 (N27W63 L= 54)   724     78        1 (1531/1661-1611/1713) -26.6 W60-W67 N24-N29
      3 1511 1602 (N24W50 L= 41)   606    200        0 (1410/1582-1592/1620) -30.3 W43-W55 N20-N27
      4 1593 1513 (N17W51 L= 42)   548    158        1 (1562/1435-1618/1570)  31.0 W44-W55 N13-N20
      5 1718 1493 (N12W61 L= 52)  1759    151        5 (1648/1427-1780/1550) -26.6 W56-W68 N08-N16
      6 1042 1519 (N32W14 L=  6)   707    172        0 ( 956/1491-1132/1540) -35.0 W09-W19 N27-N35
      7 1419 1499 (N20W38 L= 29)   765    289        1 (1297/1443-1554/1538) -31.0 W29-W44 N14-N23
      8 1342 1435 (N19W30 L= 21)   403    187        0 (1266/1393-1434/1470) -45.0 W26-W34 N14-N22
      9  457 1074 (N19E34 L=317)   484    103        3 ( 425/1037- 502/1104) -38.7 E36-E30 N16-N20
     10 1848  956 (S21W56 L= 47)   536    117        0 (1809/ 913-1897/ 991)   5.2 W51-W62 S22-S19
     11  336  904 (N12E47 L=304)   353     79        1 ( 320/ 867- 355/ 941) -90.0 E48-E45 N09-N13
     12 1709  643 (S37W40 L= 31)   358     48        1 (1683/ 623-1730/ 668)   0.0 W37-W42 S38-S35
     13 1521  583 (S35W22 L= 13)   389     75        0 (1486/ 557-1554/ 613)   0.0 W19-W25 S36-S34
     14 1384  459 (S40W09 L=  0)   466     92        2 (1340/ 431-1434/ 489)   0.0 W05-W13 S41-S38
     15 1526  301 (S61W24 L= 15)   660     71        0 (1489/ 274-1567/ 329)   5.2 W18-W31 S63-S58
     16  608  278 (S32E56 L=295)  1021     60        0 ( 574/ 249- 648/ 312) -68.2 E59-E52 S34-S29
  

    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