Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  947 1545 (N40W08 L=355)  1611    299        3 ( 850/1501-1083/1600) -40.6 W00-W16 N34-N46
      2 1339 1418 (N23W29 L= 16)   361    120        0 (1323/1363-1362/1464)  51.3 W27-W32 N19-N25
      3  328 1303 (N36E48 L=299)  1060    187        2 ( 247/1257- 407/1334) -13.4 E60-E39 N32-N38
      4 1507 1299 (N12W36 L= 23)   250     41        0 (1483/1278-1526/1324)  33.7 W34-W38 N11-N13
      5  536 1206 (N27E29 L=319)   717    148        2 ( 456/1173- 605/1232) -18.4 E35-E23 N23-N29
      6 1614  892 (S14W34 L= 21)  1285    433        4 (1544/ 764-1706/1057)  47.0 W27-W43 S20-S05
      7  471 1025 (N17E36 L=312)   318    119        0 ( 450/ 977- 517/1055)  31.0 E37-E31 N14-N18
      8  288  935 (N14E53 L=294)   672    160        3 ( 261/ 878- 320/1000) -60.3 E55-E50 N10-N17
      9  900  656 (S12E15 L=332)   292    131        1 ( 846/ 629- 956/ 696)  20.6 E18-E10 S13-S10
     10  672  617 (S11E30 L=317)   357     92        0 ( 634/ 579- 698/ 646)  31.0 E33-E27 S12-S08
  

    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