Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  795 1748 (N58W12 L=345)  3164    523        6 ( 618/1599- 986/1884)   8.1 E11-W46 N52-N62
      2 1582 1694 (N25W67 L= 40)  3133    337        7 (1462/1627-1704/1767) -51.3 W59-W77 N17-N33
      3 1573 1612 (N21W58 L= 31)  1271    295        0 (1481/1572-1698/1656) -33.7 W50-W68 N15-N25
      4  772 1339 (N30E07 L=326)   428    147        0 ( 699/1302- 817/1388)   5.7 E12-E02 N29-N32
      5 1450  770 (S20W18 L=351)   434    112        0 (1431/ 717-1464/ 824)  63.4 W16-W19 S23-S17
      6 1284  659 (S22W05 L=338)   315    153        0 (1247/ 604-1306/ 729)  45.0 W01-W07 S24-S18
      7 1788  625 (S43W48 L= 21)   510     34        1 (1769/ 599-1806/ 650) -63.4 W45-W51 S44-S40
      8 1376  475 (S37W06 L=339)  1590    309        2 (1338/ 375-1424/ 587)  52.4 W01-W11 S46-S30
      9  966  398 (S29E23 L=310)   616    134        1 ( 912/ 357-1036/ 441) -63.4 E24-E19 S33-S25
     10  553  325 (S24E57 L=276)   916    118        2 ( 495/ 289- 599/ 353) -45.0 E62-E53 S27-S20
     11  710  302 (S29E46 L=287)   564     97        1 ( 684/ 259- 727/ 349)  33.7 E50-E41 S32-S25
     12  704  231 (S34E53 L=279)   596     89        1 ( 661/ 218- 754/ 244) -20.6 E58-E48 S35-S32
  

    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