Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1123 1586 (N40W03 L=340)   441     80        2 (1078/1566-1165/1604)  18.4 W00-W06 N38-N41
      2  649 1476 (N28E31 L=306)   730    220        0 ( 567/1435- 760/1518)  -6.7 E38-E21 N26-N30
      3 1639 1387 (N29W45 L= 21)  1113    257        0 (1566/1300-1709/1479) -54.0 W39-W49 N22-N34
      4 1422 1420 (N30W26 L=  3)   320     83        0 (1384/1389-1453/1441) -21.8 W23-W28 N27-N31
      5  456 1332 (N17E42 L=294)  1543    321        1 ( 363/1241- 561/1439) -26.6 E54-E32 N12-N23
      6  947 1303 (N19E07 L=330)   624    218        1 ( 854/1242-1019/1352) -26.6 E13-E01 N15-N21
      7  500  771 (S17E34 L=302)   262     72        0 ( 482/ 736- 518/ 805)  68.2 E35-E32 S19-S14
      8 1421  692 (S16W28 L=  5)   320     72        1 (1394/ 650-1452/ 723)  53.1 W26-W30 S19-S14
      9  286  545 (S35E65 L=272)   780     78        2 ( 246/ 524- 326/ 573)   9.1 E81-E56 S37-S32
  

    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