Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  646 1728 (N38E46 L=155)   489    135        2 ( 626/1697- 663/1771) -56.3 E49-E42 N35-N41
      2  226 1362 (N07E67 L=133)  1331    242        4 ( 189/1282- 255/1413) -35.0 E76-E63 N02-N09
      3 1385 1125 (N11W20 L=221)   228     82        0 (1364/1087-1403/1167)  90.0 W19-W21 N08-N14
      4  152 1064 (S12E67 L=133)  1951    286        1 ( 120/ 992- 203/1163) -31.0 E75-E59 S16-S06
      5  657 1125 (S01E24 L=177)   331    120        1 ( 597/1102- 716/1149)  15.9 E27-E20 S02-N01
      6 1533  783 (S06W37 L=237)  1001    246        1 (1494/ 671-1560/ 889) -77.0 W32-W39 S13-N00
      7 1127  688 (S19W12 L=213)   211    111        0 (1072/ 654-1174/ 734)  50.2 W09-W14 S21-S15
      8 1118  486 (S32W17 L=218)   602    237        0 ( 991/ 467-1225/ 530)   3.4 W07-W24 S33-S30
      9 1480  481 (S26W43 L=244)   412     50        0 (1449/ 461-1510/ 504)  -9.5 W39-W45 S27-S24
  

    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