Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1346 1740 (N51W40 L=280)   685     85        3 (1319/1703-1373/1773)  71.6 W34-W44 N48-N54
      2 1857 1335 (N17W73 L=313)  1167    117        0 (1810/1302-1903/1371) -20.6 W67-W84 N13-N19
      3  361  933 (N01E46 L=194)   279     64        0 ( 344/ 902- 379/ 963) -45.0 E47-E44 S01-N02
      4  798  852 (S05E15 L=225)   406    129        1 ( 731/ 835- 852/ 876) -20.6 E19-E11 S06-S03
      5 1282  707 (S17W15 L=255)   243     87        0 (1263/ 680-1318/ 747) -53.1 W13-W17 S18-S14
      6  969  705 (S15E05 L=235)   202     47        0 ( 944/ 682- 994/ 727)  45.0 E06-E03 S16-S13
      7  323  599 (S21E57 L=182)  1904    282        5 ( 284/ 503- 370/ 714)  39.8 E71-E49 S28-S13
      8 1577  603 (S26W38 L=278)   507     80        2 (1550/ 564-1606/ 637)  63.4 W36-W41 S28-S23
      9  756  525 (S27E22 L=218)   406    163        0 ( 719/ 469- 828/ 569)  37.9 E25-E15 S30-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