Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1475 1610 (N43W12 L=339)   784    113        3 (1448/1549-1500/1669) -63.4 W09-W15 N38-N47
      2 1623 1257 (N24W30 L=357)   259     90        0 (1584/1223-1658/1288)  -9.5 W26-W33 N23-N25
      3 1781 1232 (N28W44 L= 11)   326     89        0 (1763/1192-1795/1277) -45.0 W41-W46 N26-N30
      4  378 1149 (S15E43 L=284)   276     37        0 ( 356/1128- 399/1171) -14.0 E45-E41 S16-S14
      5 1574 1076 (N12W30 L=357)   343     78        1 (1530/1059-1612/1108)  18.4 W27-W33 N10-N13
      6 1031  796 (S19W07 L=334)   885    263        2 ( 935/ 706-1126/ 860)  54.0 W03-W11 S25-S13
      7  336  638 (S48E44 L=283)   648    106        0 ( 297/ 613- 383/ 673)  -5.2 E49-E37 S49-S45
      8 1679  529 (S13W60 L= 27)   812    163        1 (1636/ 487-1736/ 585)  76.0 W55-W62 S15-S07
      9 1398  469 (S25W42 L=  9)   427    110        0 (1347/ 431-1436/ 493)   7.1 W38-W46 S26-S23
     10 1191  373 (S38W35 L=  2)   322     92        1 (1146/ 360-1247/ 389)  51.3 W32-W38 S39-S34
  

    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