Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1256 1664 (N39W04 L=238)  1378    395        0 (1098/1580-1420/1755) -13.0 E09-W17 N35-N43
      2  366 1626 (N23E65 L=169)  1013     82        3 ( 330/1592- 396/1656)  26.6 E70-E59 N20-N25
      3  492 1125 (S08E35 L=199)   251     74        0 ( 469/1090- 521/1154)  63.4 E36-E33 S10-S05
      4  485  888 (S22E33 L=200)   776    242        1 ( 432/ 816- 531/ 975)  73.3 E37-E29 S26-S16
      5  139  835 (S28E73 L=161)   698     47        1 ( 108/ 822- 168/ 854)   4.1 E81-E67 S28-S26
      6  773  806 (S23E12 L=221)   470    202        0 ( 686/ 751- 855/ 850)  38.7 E17-E07 S27-S19
      7 1261  788 (S16W18 L=252)   310     94        1 (1208/ 763-1296/ 819)  45.0 W15-W20 S18-S13
      8 1249  661 (S24W21 L=255)  2716    821        4 ( 947/ 552-1623/ 810)  34.7 W03-W42 S35-S08
      9 1699  641 (S15W55 L=288)   668     91        0 (1655/ 600-1732/ 679)  56.3 W52-W57 S18-S12
     10  828  586 (S36E05 L=228)   416     98        0 ( 783/ 559- 877/ 616)   0.0 E09-E01 S37-S34
     11 1712  494 (S23W66 L=300)   651     76        0 (1665/ 472-1748/ 509)   4.4 W59-W72 S23-S21
  

    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