Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1377 1782 (N37W53 L= 97)  1281    226        3 (1302/1751-1450/1804) -20.0 W46-W58 N33-N40
      2 1737 1494 (N10W63 L=108)  1543    248        5 (1672/1411-1784/1587)  85.2 W59-W66 N04-N16
      3 1177 1516 (N28W24 L= 68)   475    134        0 (1119/1483-1249/1533)  -6.3 W19-W28 N25-N29
      4 1269 1333 (N15W22 L= 67)   254     70        0 (1242/1309-1292/1361)  26.6 W20-W24 N13-N16
      5 1049 1289 (N19W08 L= 53)   365     55        1 (1019/1262-1087/1313) -45.0 W06-W10 N17-N20
      6  212 1280 (N39E54 L=351)   787    232        0 ( 161/1263- 265/1306) -20.0 E59-E47 N37-N42
      7 1877 1077 (S17W59 L=104)   545    100        3 (1849/1045-1900/1108) -59.0 W56-W62 S18-S13
      8 1798  951 (S21W48 L= 92)   485     82        0 (1775/ 908-1818/ 989)  78.7 W46-W49 S23-S18
      9 1473  701 (S26W18 L= 63)   999    249        2 (1439/ 583-1527/ 814)  55.3 W15-W24 S33-S20
     10 1149  596 (S23E04 L= 40)   928    291        1 (1055/ 492-1198/ 727)  35.5 E12-W02 S27-S16
     11  591  303 (S27E53 L=352)   764    155        2 ( 566/ 260- 616/ 344)  36.9 E56-E48 S30-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