Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  842 1828 (N39E41 L=221)   360     64        0 ( 820/1810- 866/1845)  45.0 E43-E38 N37-N41
      2 1019 1802 (N42E26 L=236)   389     75        0 (1000/1774-1051/1830)  90.0 E27-E24 N39-N45
      3  677 1765 (N30E48 L=214)  1525    166        8 ( 611/1713- 745/1822)  73.3 E51-E44 N25-N36
      4 1504 1795 (N64W19 L=281)   620     53        0 (1486/1767-1522/1821)  47.5 W15-W27 N58-N70
      5  804 1691 (N28E33 L=228)  1130    245        3 ( 740/1649- 877/1758)  32.0 E38-E28 N24-N32
      6  411 1572 (N12E57 L=205)   578    136        0 ( 377/1549- 465/1607)   6.3 E60-E51 N10-N14
      7 1538 1449 (N32W22 L=284)  1660    445        0 (1409/1304-1692/1605) -24.8 W10-W36 N26-N39
      8  492 1534 (N11E47 L=214)   394     67        2 ( 458/1512- 523/1553)  36.9 E50-E45 N09-N12
      9 1642 1469 (N37W31 L=293)   336     84        0 (1605/1433-1673/1502) -15.9 W27-W34 N35-N38
     10 1802 1329 (N33W48 L=310)   535     56        0 (1771/1305-1824/1354)  -8.1 W44-W51 N31-N34
     11 1205 1161 (N05W07 L=269)   345     97        1 (1158/1132-1243/1189)  51.3 W04-W08 N02-N07
     12  458 1055 (S18E36 L=226)  1274    390        2 ( 326/ 912- 568/1170)  51.7 E45-E30 S28-S09
     13  158 1007 (S25E66 L=196)  2175    167        4 ( 130/ 922- 190/1088)  79.7 E71-E61 S30-S19
     14 1196  879 (S11W13 L=275)   362    112        0 (1127/ 858-1239/ 897)  15.9 W09-W16 S12-S09
     15  538  775 (S33E26 L=236)  2470    564        4 ( 321/ 638- 770/ 862)  27.1 E48-E09 S45-S24
     16 1307  525 (S28W33 L=294)   530    180        0 (1226/ 507-1387/ 543)  14.0 W26-W38 S29-S25
     17 1573  492 (S21W53 L=315)   487     87        0 (1533/ 465-1617/ 524)  59.0 W51-W55 S23-S18
     18 1485  358 (S31W57 L=319)  2390    287        2 (1347/ 333-1607/ 395)  26.6 W45-W68 S36-S25
  

    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