Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1198 1561 (N35W26 L=171)   375     89        1 (1162/1536-1246/1578) -18.4 W23-W29 N33-N36
      2  550 1370 (N37E26 L=119)   851    172        1 ( 501/1316- 611/1424) -45.0 E30-E21 N32-N41
      3  316 1130 (N26E47 L= 98)  1745    358        3 ( 271/ 947- 375/1249) -73.3 E51-E43 N13-N33
      4 1748 1170 (S02W51 L=196)   393     73        3 (1721/1146-1779/1193) -45.0 W48-W53 S03-N00
      5 1689 1149 (S01W45 L=191)   303     64        0 (1673/1113-1705/1186)  26.6 W43-W47 S03-N00
      6  852 1122 (N16E08 L=137)   300     90        1 ( 808/1104- 908/1136)  -8.1 E11-E04 N14-N17
      7 1066  915 (S00W00 L=145)  1111    268        3 ( 980/ 790-1118/1024) -90.0 E02-W03 S07-N07
      8  144  981 (N18E68 L= 77)   566    106        3 ( 131/ 941- 160/1017) -90.0 E69-E65 N15-N20
      9 1691  948 (S13W41 L=186)   388     53        2 (1665/ 921-1717/ 975)  45.0 W38-W43 S14-S11
     10 1030  620 (S17E08 L=137)   925    189        3 ( 953/ 535-1072/ 671) -74.7 E11-E05 S23-S12
     11  829  312 (S34E32 L=113)  1095    193        3 ( 782/ 259- 871/ 362)  26.6 E38-E26 S38-S31
  

    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