Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  678 1601 (N36E33 L= 66)   415    131        0 ( 633/1554- 717/1640)  66.8 E34-E31 N32-N39
      2  297 1517 (N27E68 L= 30)  2620    213        7 ( 229/1446- 348/1590)  90.0 E78-E59 N21-N31
      3  390 1533 (N29E57 L= 42)   383     37        2 ( 367/1515- 410/1552)  45.0 E59-E54 N27-N30
      4 1260 1421 (N29W13 L=112)   873    265        0 (1179/1316-1312/1527) -60.3 W07-W16 N21-N35
      5  998 1441 (N28E05 L= 94)  1176    213        5 ( 926/1360-1066/1510) -48.0 E10-W00 N22-N32
      6  361 1452 (N23E55 L= 44)   347     51        1 ( 341/1427- 379/1473)  56.3 E57-E53 N21-N24
      7  450 1362 (N18E43 L= 56)   296     93        0 ( 410/1341- 482/1381)   9.5 E46-E40 N16-N19
      8 1571 1302 (N23W36 L=135)   381     79        1 (1548/1260-1595/1338) -68.2 W34-W38 N20-N25
      9 1844 1146 (N14W62 L=161)   968     83        2 (1821/1095-1868/1188)  80.5 W59-W65 N10-N16
     10  191  932 (S11E63 L= 36)   665     73        2 ( 169/ 895- 229/ 966)  68.2 E66-E57 S13-S08
     11  405  923 (S09E41 L= 58)   437    158        0 ( 344/ 872- 466/ 958)  31.0 E46-E35 S12-S06
     12  165  878 (S14E67 L= 31)  1126    117        4 ( 136/ 833- 207/ 926) -29.1 E72-E61 S16-S11
     13 1221  676 (S18W15 L=114)  3407    766        6 ( 928/ 512-1446/ 921)  44.0 E02-W27 S31-S01
     14  745  800 (S14E16 L= 83)   265     84        0 ( 699/ 786- 782/ 814)   9.5 E19-E13 S14-S12
     15 1782  729 (S12W59 L=158)   683     85        0 (1755/ 695-1821/ 760)  33.7 W56-W63 S13-S09
     16  320  627 (S29E53 L= 46)  1741    222        3 ( 248/ 541- 401/ 684)  20.1 E73-E43 S36-S24
     17 1718  577 (S22W59 L=157)  1997    283        4 (1656/ 533-1805/ 625)  26.6 W51-W68 S25-S18
     18  617  431 (S41E29 L= 70)  1136    141        2 ( 565/ 386- 683/ 478)  29.7 E35-E21 S45-S36
  

    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