Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  183 1313 (N38E62 L=261)  1291    188        3 ( 155/1275- 216/1361)  78.7 E64-E57 N35-N40
      2 1030 1272 (N21W06 L=329)  1294    309        4 ( 881/1224-1170/1316) -36.3 E03-W13 N15-N26
      3  555 1234 (N28E27 L=296)   610     86        2 ( 517/1200- 594/1270) -45.0 E29-E23 N25-N30
      4 1841 1187 (S03W62 L= 25)  1103    168        6 (1810/1123-1875/1244) -74.1 W59-W65 S07-N00
      5  179 1094 (N24E62 L=262)  2488    395        3 ( 132/ 926- 246/1210) -71.6 E67-E56 N13-N31
      6 1443 1075 (N01W26 L=349)   465    106        0 (1396/1033-1488/1112)  15.9 W22-W29 N00-N02
      7 1200  901 (S04W08 L=331)  1768    532        3 (1100/ 755-1320/1062)  47.1 W01-W17 S12-N04
      8 1313  763 (S14W12 L=335)   234     79        0 (1291/ 735-1330/ 808)  53.1 W10-W13 S16-S11
      9  269  580 (S09E67 L=256)  3112    471        5 ( 210/ 518- 335/ 639) -69.4 E76-E58 S13-S04
     10  728  525 (S18E29 L=294)  1195    255        4 ( 677/ 428- 782/ 622)  45.0 E35-E23 S24-S12
     11  403  505 (S15E56 L=268)   997    163        2 ( 353/ 432- 453/ 564) -78.7 E58-E52 S20-S10
     12  615  303 (S32E49 L=275)  1149    170        1 ( 551/ 280- 657/ 341)   5.2 E54-E43 S34-S29
  

    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