Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1431 1653 (N49W23 L= 45)   314     62        0 (1400/1642-1464/1666)   8.1 W19-W26 N47-N50
      2 1054 1540 (N33E07 L= 16)  1565    439        4 ( 879/1424-1255/1632)  -9.8 E21-W08 N27-N36
      3  640 1586 (N28E38 L=345)   481     70        2 ( 614/1545- 669/1619) -78.7 E39-E35 N25-N30
      4 1598 1463 (N37W36 L= 59)   333     99        1 (1554/1448-1638/1480)   7.1 W31-W39 N36-N38
      5  582 1389 (N15E35 L=348)   308    154        0 ( 543/1341- 639/1436) -29.1 E38-E29 N12-N17
      6 1529 1281 (N23W29 L= 52)   382     60        1 (1503/1260-1558/1308) -33.7 W27-W31 N21-N25
      7  621  898 (S14E23 L=360)   345    139        0 ( 579/ 860- 670/ 954) -29.7 E26-E19 S16-S10
      8 1365  416 (S33W37 L= 60)   402     85        0 (1332/ 399-1400/ 435)   0.0 W33-W39 S34-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