Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1078 1630 (N38W24 L=160)   867    217        4 (1039/1555-1143/1703)  19.7 W18-W32 N33-N41
      2 1491 1605 (N23W49 L=185)  1178    187        4 (1423/1548-1569/1634) -56.3 W43-W53 N17-N26
      3  548 1254 (N31E24 L=113)   973    401        1 ( 395/1226- 728/1294) -19.8 E35-E10 N26-N36
      4 1546 1173 (S01W34 L=170)  1250    279        5 (1444/1130-1648/1216) -19.7 W26-W40 S05-N02
      5  808 1081 (N14E11 L=126)   265     64        0 ( 780/1059- 842/1116) -63.4 E11-E08 N12-N16
      6 1709  992 (S16W41 L=177)   601     94        2 (1664/ 964-1744/1022)  -8.1 W37-W44 S17-S14
      7 1449  818 (S18W19 L=155)   556    151        1 (1423/ 759-1493/ 908) -84.3 W17-W20 S21-S11
      8 1109  700 (S15E04 L=132)   745    226        1 ( 984/ 655-1200/ 727) -42.0 E10-W00 S19-S10
      9  740  505 (S17E32 L=105)   258     59        1 ( 715/ 488- 779/ 524) -14.0 E33-E29 S17-S15
     10  679  440 (S19E38 L= 98)   382    120        1 ( 640/ 412- 706/ 475) -90.0 E41-E36 S21-S16
     11 1049  226 (S46E28 L=108)  1018    155        2 (1015/ 180-1074/ 271)  31.0 E33-E23 S49-S41
  

    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