Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1361 1690 (N45W01 L= 13)  1955    320        5 (1253/1575-1484/1799) -25.6 E11-W13 N39-N51
      2  736 1696 (N25E40 L=333)   706     99        2 ( 687/1665- 773/1729)  80.5 E43-E36 N21-N27
      3 1169 1598 (N31E09 L=  4)   378     98        1 (1153/1545-1181/1646) -50.2 E11-E06 N27-N33
      4  824 1565 (N19E28 L=345)   255     96        1 ( 795/1544- 868/1598)   0.0 E30-E24 N17-N20
      5 1406 1460 (N29W10 L= 23)   262    129        0 (1352/1418-1455/1489)  -7.1 W06-W14 N28-N29
      6 1856 1182 (N29W54 L= 67)   509     90        1 (1829/1147-1886/1224) -20.6 W50-W59 N26-N30
      7  183 1141 (S20E65 L=308)   828    111        3 ( 166/1090- 201/1200) -74.1 E67-E61 S22-S15
      8  802  866 (S22E08 L=  4)   267     88        1 ( 752/ 849- 845/ 880)  26.6 E11-E05 S23-S20
      9 1405  399 (S29W48 L= 61)   950    221        0 (1355/ 317-1483/ 456) -79.7 W45-W53 S34-S23
  

    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