Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  959 1900 (N62E44 L=111)   785     59        0 ( 939/1874- 981/1932) -14.0 E59-E35 N59-N65
      2  998 1866 (N59E32 L=123)   426     40        0 ( 974/1849-1020/1885)  63.4 E35-E29 N56-N61
      3  657 1578 (N27E37 L=118)   542    167        0 ( 611/1552- 723/1626)  -5.2 E43-E32 N25-N29
      4  324 1549 (N20E68 L= 87)  1026     95        0 ( 280/1524- 370/1570)   0.0 E78-E61 N18-N20
      5  557 1465 (N18E40 L=115)   419     60        0 ( 537/1435- 582/1491) -31.0 E42-E37 N16-N19
      6 1880 1278 (N31W70 L=225)  1080     98        1 (1839/1243-1912/1314) -12.1 W63-W78 N28-N32
      7 1570 1265 (N23W32 L=187)   413     98        0 (1555/1213-1585/1312) -71.6 W31-W33 N20-N26
      8 1821  918 (N07W59 L=214)   452    101        0 (1800/ 871-1841/ 963)  63.4 W57-W61 N03-N09
      9  239  848 (S25E57 L= 98)   805    122        0 ( 216/ 793- 271/ 900)  53.1 E60-E52 S28-S20
     10  378  800 (S25E42 L=113)   338     95        0 ( 344/ 785- 423/ 828)   0.0 E44-E37 S26-S23
     11 1044  606 (S26W09 L=164)  2511    569        6 ( 811/ 508-1258/ 732)  41.2 E05-W19 S35-S14
     12 1549  617 (S16W43 L=198)   670    140        1 (1500/ 573-1611/ 655)  54.5 W40-W47 S19-S12
     13  610  260 (S66E31 L=124)   522     64        0 ( 591/ 244- 642/ 286) -18.4 E37-E23 S68-S63
  

    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