Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  643 1478 (N31E30 L= 37)   525    130        0 ( 570/1461- 701/1492)  11.3 E35-E25 N29-N32
      2 1323 1322 (N22W19 L= 86)   273    109        0 (1269/1290-1367/1340) -23.2 W15-W22 N20-N23
      3 1787 1277 (N19W59 L=125)   382     51        1 (1769/1246-1800/1300) -71.6 W57-W61 N17-N20
      4  896 1254 (N17E09 L= 58)   476    132        1 ( 835/1228- 950/1271)  15.9 E12-E05 N15-N18
      5 1650 1209 (N15W43 L=110)   337    111        0 (1600/1188-1682/1236)   8.1 W39-W46 N13-N16
      6 1779  736 (S15W58 L=125)  1284    218        2 (1743/ 651-1835/ 822)  65.6 W54-W62 S20-S09
      7  257  770 (S16E58 L=  9)   670    116        0 ( 236/ 718- 281/ 817)  54.5 E61-E54 S19-S12
      8  556  623 (S24E32 L= 34)   599    100        0 ( 513/ 582- 595/ 665) -39.8 E35-E29 S26-S21
      9 1251  443 (S35W19 L= 85)  1037    221        2 (1220/ 342-1293/ 548) -66.4 W15-W24 S43-S27
  

    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