Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  911 1374 (N28E05 L=157)   358     68        1 ( 870/1348- 952/1390)   0.0 E07-E01 N26-N28
      2 1540 1288 (N16W37 L=199)   354    111        0 (1479/1277-1595/1308) -20.6 W33-W41 N15-N18
      3  708 1039 (N08E20 L=142)   379    104        1 ( 674/1009- 755/1078) -21.8 E21-E16 N06-N10
      4 1903 1018 (S05W70 L=232)  1062     84        1 (1862/ 998-1947/1041) -12.5 W63-W81 S07-S03
      5 1808  829 (S15W57 L=219)   583    105        1 (1768/ 794-1843/ 853) -16.7 W52-W62 S18-S13
      6  706  703 (S13E23 L=139)   458    161        0 ( 641/ 678- 755/ 747)   7.1 E27-E19 S14-S10
      7 1386  652 (S22W21 L=183)  1401    337        3 (1261/ 555-1509/ 735)  26.6 W12-W30 S27-S16
      8 1744  632 (S28W54 L=216)   471     86        0 (1713/ 593-1777/ 659)  33.7 W51-W57 S30-S26
      9 1218  464 (S33W09 L=171)   446     55        0 (1188/ 440-1247/ 495)  36.9 W06-W11 S35-S31
     10  594  339 (S39E45 L=117)   544     80        2 ( 575/ 297- 611/ 381)  37.9 E50-E40 S42-S35
     11 1066  268 (S49E05 L=157)   508     82        0 (1036/ 246-1102/ 301) -31.0 E06-E01 S51-S45
  

    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