Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1391 1524 (N29W37 L=229)   811    311        2 (1273/1479-1532/1569) -33.7 W30-W46 N23-N34
      2 1689 1428 (N16W57 L=249)   356     45        1 (1676/1398-1703/1457)  63.4 W55-W59 N13-N17
      3 1019 1326 (N25W06 L=197)   928    260        1 ( 935/1228-1082/1422) -71.6 W01-W07 N17-N32
      4  811 1042 (N12E13 L=179)   260    161        1 ( 777/ 995- 883/1069)  14.0 E15-E07 N09-N13
      5 1760  948 (S13W49 L=240)   310     62        1 (1737/ 922-1784/ 977)  26.6 W46-W51 S14-S11
      6  937  803 (S05E09 L=182)   397    151        1 ( 886/ 770-1010/ 843) -45.0 E11-E05 S07-S01
      7 1262  591 (S24W07 L=199)  1094    476        1 (1112/ 457-1426/ 744)  28.4 E05-W19 S30-S17
      8 1054  631 (S18E05 L=186)   223     90        1 (1017/ 601-1080/ 666)  31.0 E08-E03 S18-S15
      9 1454  611 (S27W21 L=212)   420    116        1 (1411/ 574-1510/ 655)  29.7 W17-W25 S29-S25
     10  332  534 (S14E63 L=129)   980    100        2 ( 285/ 500- 377/ 561) -10.3 E68-E57 S16-S11
     11  796  471 (S24E26 L=165)   350    104        0 ( 763/ 443- 823/ 522)  39.8 E29-E23 S25-S20
  

    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