Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  651 1710 (N29E47 L=136)  1519    310        2 ( 532/1672- 740/1745)   6.3 E57-E39 N26-N31
      2  381 1571 (N14E61 L=122)  2648    352        9 ( 309/1477- 444/1655)  71.6 E67-E55 N07-N19
      3 1139 1576 (N32E08 L=175)   306     95        0 (1094/1536-1178/1601) -15.9 E11-E04 N30-N33
      4 1663  749 (S02W48 L=231)   706    138        4 (1633/ 690-1684/ 802) -51.3 W45-W50 S05-N00
      5 1337  712 (S13W27 L=209)   394    127        0 (1288/ 654-1398/ 767)  71.6 W25-W28 S17-S08
      6  742  544 (S40E05 L=177)   971    270        1 ( 626/ 517- 888/ 564)  17.5 E14-W05 S42-S36
      7 1506  511 (S20W46 L=229)   407    103        0 (1455/ 492-1546/ 539)  33.7 W42-W48 S22-S17
      8  327  518 (S54E48 L=135)   504     87        0 ( 294/ 500- 356/ 535)   0.0 E52-E44 S55-S53
      9 1250  401 (S34W35 L=218)   379    109        1 (1207/ 377-1308/ 431)  14.0 W30-W38 S35-S31
     10 1587  383 (S25W61 L=243)  1601    394        2 (1497/ 338-1654/ 413)   6.0 W52-W71 S27-S22
  

    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