Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1309 1666 (N35W43 L=158)  2365    454        3 (1127/1601-1487/1697) -12.3 W26-W58 N28-N39
      2 1107 1517 (N33W20 L=135)   655    186        1 (1071/1427-1150/1593)  70.0 W17-W24 N26-N37
      3  443 1066 (N23E36 L= 79)   283     67        0 ( 413/1047- 475/1092) -38.7 E38-E33 N21-N25
      4 1924 1011 (S21W68 L=183)   781     60        0 (1900/ 982-1950/1043) -36.9 W63-W75 S24-S18
      5 1670  857 (S21W36 L=151)   313     32        0 (1654/ 830-1686/ 884)  45.0 W34-W37 S21-S19
      6 1602  794 (S22W29 L=144)   290     76        1 (1563/ 782-1641/ 808) -18.4 W26-W32 S23-S20
      7 1692  703 (S31W36 L=151)  1649    287        3 (1616/ 598-1795/ 807)  17.5 W28-W47 S36-S26
      8 1255  658 (S20W04 L=119)   392    117        0 (1230/ 595-1275/ 715)  49.4 W00-W06 S23-S16
      9  498  397 (S18E55 L= 59)   507     75        1 ( 475/ 361- 533/ 428)  59.0 E58-E52 S20-S15
  

    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