Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1106 1703 (N40W01 L= 73)  1802    398        4 ( 940/1638-1277/1768) -13.0 E12-W14 N35-N45
      2  727 1623 (N32E25 L= 46)   746    145        0 ( 685/1562- 767/1675) -45.0 E29-E21 N27-N35
      3 1294 1624 (N36W16 L= 87)   366     43        0 (1274/1597-1316/1649)  53.1 W14-W17 N33-N37
      4 1355 1597 (N34W20 L= 92)   444     71        1 (1330/1567-1385/1622)  38.7 W17-W22 N31-N36
      5 1528 1528 (N31W33 L=104)  1407    300        3 (1401/1469-1658/1575)  -5.7 W23-W43 N27-N33
      6 1230 1285 (N12W11 L= 83)   479    168        0 (1210/1189-1257/1365) -79.7 W09-W13 N06-N17
      7  588 1313 (N11E29 L= 43)   302     66        0 ( 548/1292- 617/1327) -21.8 E31-E26 N09-N11
      8 1817 1221 (N14W55 L=126)   444    111        0 (1784/1192-1859/1254) -23.2 W51-W59 N12-N15
      9  276 1189 (N03E52 L= 19)   288     74        0 ( 247/1161- 304/1219)  31.0 E54-E49 N01-N04
     10  652  686 (S27E23 L= 49)   248     68        0 ( 626/ 651- 677/ 720)  59.0 E25-E21 S29-S24
     11  259  583 (S34E59 L= 13)   639     92        1 ( 240/ 556- 295/ 619) -68.2 E60-E55 S36-S31
     12 1695  542 (S28W54 L=126)   418     69        1 (1667/ 518-1729/ 564)  36.9 W52-W56 S30-S26
  

    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