Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  791 1360 (N27E06 L=335)   242     64        0 ( 757/1338- 830/1381) -45.0 E08-E04 N25-N29
      2 1842 1065 (S14W56 L= 38)  4759    792       12 (1797/ 677-1893/1342)  74.2 W50-W67 S38-N01
      3 1705 1175 (S04W46 L= 27)  2176    424        1 (1594/ 987-1765/1307) -76.6 W40-W49 S16-N05
      4  789 1250 (N21E09 L=333)   517    213        1 ( 683/1228- 869/1274) -14.0 E15-E03 N17-N22
      5 1139 1215 (N11W11 L=352)   229     83        1 (1097/1196-1184/1234) -36.9 W08-W13 N09-N12
      6 1314  876 (S12W13 L=355)   401     93        2 (1294/ 835-1343/ 919)  53.1 W11-W14 S14-S09
      7 1371  876 (S13W17 L=358)   218     22        0 (1347/ 860-1393/ 893)  14.0 W14-W18 S14-S11
      8 1063  554 (S25E09 L=333)  4343   1134        6 ( 733/ 187-1440/ 822)  18.4 E53-W19 S43-S17
      9  888  732 (S11E14 L=328)   284    167        0 ( 812/ 692- 934/ 785)  18.4 E19-E10 S11-S08
     10  460  461 (S18E53 L=289)   762    135        2 ( 419/ 410- 512/ 513) -82.9 E56-E50 S21-S13
     11  485  295 (S29E65 L=277)   705     58        2 ( 463/ 269- 520/ 324)  13.0 E72-E59 S30-S27
  

    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