Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  607 1327 (N30E24 L=208)   379     88        1 ( 563/1305- 654/1344) -15.9 E27-E20 N28-N31
      2 1207 1201 (N14W14 L=246)   659    208        0 (1100/1145-1287/1228) -38.7 W07-W18 N09-N17
      3  175 1193 (N25E63 L=169)  1500    333        3 ( 116/1163- 231/1215)   0.0 E71-E56 N23-N26
      4  435 1059 (N15E39 L=193)  2176    459        9 ( 384/ 862- 511/1213) -75.4 E42-E34 N02-N25
      5  533 1166 (N21E30 L=201)   325     50        0 ( 505/1138- 560/1195) -45.0 E32-E28 N19-N23
      6 1145  958 (N01W06 L=238)   413    104        0 (1087/ 921-1189/ 991) -45.0 W03-W08 S02-N03
      7  240  872 (N05E59 L=173)  1122    208        4 ( 191/ 807- 282/ 944) -53.1 E63-E54 N01-N09
      8 1004  565 (S22E08 L=224)   752    248        1 ( 960/ 470-1032/ 690)  86.2 E10-E06 S28-S13
      9 1417  255 (S58W24 L=256)   593     80        0 (1379/ 239-1459/ 270)   0.0 W19-W29 S59-S56
  

    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