Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1606 1720 (N29W74 L= 44)  1045     78        4 (1566/1692-1638/1751) -59.0 W70-W80 N25-N30
      2 1530 1666 (N29W59 L= 28)   569     46        1 (1510/1638-1555/1694)  26.6 W55-W62 N27-N30
      3  581 1644 (N58E20 L=309)   329     64        0 ( 558/1620- 610/1662)  18.4 E23-E16 N56-N59
      4 1083 1584 (N38W21 L=350)   506    114        1 (1033/1539-1123/1610) -63.4 W18-W23 N34-N40
      5 1657 1556 (N19W61 L= 30)  2188    308        5 (1497/1519-1776/1602) -18.4 W46-W76 N13-N23
      6  296 1388 (N45E53 L=277)   983    137        0 ( 234/1360- 357/1407)  -8.5 E63-E43 N41-N46
      7  504 1324 (N36E29 L=301)   281     49        1 ( 479/1299- 533/1346)   0.0 E30-E25 N34-N37
      8 1394  837 (S12W17 L=346)  1067    280        2 (1349/ 725-1428/ 983)  55.3 W12-W21 S18-S04
      9  743  722 (S04E23 L=306)   423    129        0 ( 719/ 663- 773/ 785)  45.0 E25-E19 S06-N00
     10 1027  366 (S33E18 L=311)  2071    528        5 ( 779/ 247-1152/ 499)  15.9 E41-E06 S37-S26
     11  573  423 (S19E46 L=283)   484     67        1 ( 547/ 389- 596/ 448)  76.0 E48-E43 S21-S17
     12  786  392 (S25E32 L=297)   397    112        2 ( 754/ 352- 819/ 435)  32.0 E36-E28 S27-S22
  

    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