Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1806 1253 (S04W58 L=328)   343    130        0 (1761/1232-1837/1271) -23.2 W53-W62 S05-S01
      2 1441 1222 (N06W29 L=299)   240    115        0 (1389/1193-1473/1264)   8.1 W25-W32 N05-N07
      3  985 1155 (N15W01 L=271)   306    111        0 ( 938/1109-1014/1207)  33.7 E02-W04 N13-N17
      4 1167 1186 (N12W12 L=282)   725    268        0 (1049/1154-1311/1203) -17.4 W04-W20 N08-N14
      5 1701  938 (S17W39 L=309)   646    139        2 (1663/ 907-1753/ 973) -39.8 W36-W43 S19-S14
      6  468  913 (N14E37 L=233)  2120    566       10 ( 275/ 855- 686/ 968) -23.2 E51-E22 N07-N19
      7 1571  799 (S21W27 L=297)   336    107        0 (1556/ 745-1589/ 848)  63.4 W25-W28 S24-S18
      8 1071  712 (S12E06 L=264)   264    109        0 (1051/ 653-1102/ 761)  45.0 E08-E02 S14-S09
      9 1540  560 (S36W21 L=291)   341     70        0 (1512/ 542-1570/ 577) -36.9 W18-W22 S38-S34
     10  553  354 (S22E52 L=217)  2652    342        3 ( 460/ 225- 648/ 473)  81.9 E56-E49 S33-S11
  

    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