Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1148 1883 (N55E30 L= 42)   645     40        1 (1128/1858-1180/1906)  90.0 E35-E27 N51-N57
      2  584 1657 (N18E49 L= 23)  1298    209        1 ( 543/1573- 659/1734) -24.2 E58-E38 N14-N23
      3 1217 1231 (N10W05 L= 77)   339    110        2 (1195/1176-1240/1273) -59.0 W03-W07 N07-N12
      4 1307 1194 (N11W11 L= 83)   226    112        1 (1292/1139-1334/1235) -59.0 W09-W12 N07-N13
      5  318 1174 (S16E50 L= 22)   467    132        1 ( 298/1111- 334/1233) -54.5 E53-E47 S19-S12
      6 1604 1120 (N16W31 L=103)   290     62        1 (1591/1080-1622/1157) -63.4 W29-W32 N13-N17
      7  169 1112 (S22E66 L=  6)   751     85        1 ( 150/1077- 190/1144)  51.3 E69-E63 S24-S19
      8  200 1042 (S26E61 L= 11)  1899    188        4 ( 156/ 966- 242/1127)  55.0 E66-E55 S30-S20
      9  374 1018 (S23E42 L= 30)   822    226        0 ( 273/ 966- 470/1061)  32.0 E50-E34 S28-S18
     10 1069  997 (S07W03 L= 75)   362    108        1 (1038/ 959-1120/1026)  45.0 W00-W06 S09-S05
  

    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