Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  366 1352 (N35E47 L= 64)   418     69        0 ( 330/1338- 401/1368)  -7.1 E51-E43 N33-N35
      2 1271 1298 (N19W20 L=132)   722    276        0 (1204/1243-1354/1339) -41.2 W15-W24 N14-N22
      3 1717 1041 (S04W47 L=158)   479     84        0 (1676/1016-1758/1067) -29.7 W43-W50 S06-S02
      4  978  615 (S18E09 L=102)   432    109        1 ( 938/ 568-1020/ 663)  35.5 E12-E05 S20-S15
      5  647  486 (S23E35 L= 77)   676    195        0 ( 573/ 434- 717/ 552)  21.8 E42-E27 S25-S18
      6 1165  472 (S31W01 L=112)   437    141        1 (1117/ 421-1210/ 529)  45.0 E03-W04 S34-S27
      7  752  394 (S31E31 L= 81)   837    211        0 ( 666/ 344- 839/ 448)  16.4 E39-E22 S33-S27
      8 1239  407 (S38W05 L=117)   456    106        1 (1187/ 375-1291/ 441)  18.4 W00-W09 S39-S35
  

    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