Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  541 1576 (N36E33 L=327)   740    137        1 ( 505/1518- 601/1640)  58.0 E35-E29 N32-N40
      2 1556 1336 (N14W36 L= 36)   264     37        2 (1528/1315-1579/1353) -14.0 W34-W38 N12-N15
      3  836 1145 (N06E11 L=350)   794    225        1 ( 746/1036- 892/1232) -57.5 E15-E07 N00-N11
      4  445 1098 (N06E36 L=324)   290     71        1 ( 412/1076- 476/1123) -14.0 E38-E34 N04-N07
      5  447 1034 (N02E37 L=324)   525    121        1 ( 410/1006- 489/1071) -33.7 E39-E33 N00-N04
      6 1777  680 (S26W54 L= 54)   701    109        0 (1742/ 623-1816/ 739)  60.3 W51-W56 S29-S22
      7 1409  509 (S36W25 L= 26)   789    199        1 (1316/ 483-1506/ 547)   8.1 W18-W32 S37-S34
      8 1192  384 (S45W09 L=  9)   329     81        1 (1153/ 372-1240/ 401)   0.0 W05-W13 S45-S43
  

    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