Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1563 1317 (N06W40 L=146)   462    132        1 (1537/1258-1606/1360)  26.6 W36-W44 N03-N08
      2 1734 1213 (S04W50 L=155)  1272    309        2 (1634/1120-1836/1329) -60.6 W45-W57 S12-N04
      3  827 1240 (N21E06 L=100)   249     91        1 ( 804/1193- 859/1282)  45.0 E07-E03 N18-N22
      4  245 1001 (N20E53 L= 53)  2163    544        3 ( 114/ 903- 409/1073) -31.0 E66-E40 N11-N26
      5 1532 1002 (S10W29 L=135)   489    153        0 (1513/ 915-1547/1070)  63.4 W26-W31 S14-S06
      6  762 1050 (N12E14 L= 91)   259    106        0 ( 708/1032- 816/1067)  -8.1 E17-E10 N10-N12
      7 1774  935 (S21W46 L=152)   335    118        1 (1739/ 907-1815/ 947)   0.0 W42-W50 S22-S19
      8 1196  792 (S13W04 L=110)   368    124        0 (1126/ 758-1246/ 832)  15.9 W00-W07 S13-S11
      9 1647  696 (S32W32 L=138)   309     54        0 (1627/ 663-1667/ 721)  36.9 W30-W34 S33-S30
     10  909  543 (S21E20 L= 86)   558    250        2 ( 889/ 432- 937/ 643)  52.4 E25-E15 S27-S14
     11 1366  571 (S31W09 L=115)   479     99        0 (1346/ 543-1393/ 618) -90.0 W07-W11 S33-S27
     12  338  567 (S06E59 L= 47)  1057    175        1 ( 278/ 534- 422/ 602) -26.6 E64-E51 S09-S03
  

    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