Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1179 1673 (N41W34 L=239)   540     97        0 (1133/1657-1231/1689) -33.7 W31-W37 N38-N42
      2 1440 1515 (N23W42 L=247)   415     94        2 (1400/1482-1477/1551) -71.6 W40-W43 N20-N26
      3 1300 1397 (N21W27 L=232)   966    369        1 (1138/1312-1449/1454) -48.8 W18-W33 N12-N28
      4  987 1237 (N20W03 L=208)  2409    511        6 ( 758/1179-1205/1271) -32.0 E10-W14 N11-N26
      5  405 1179 (N30E39 L=166)   724    210        3 ( 368/1107- 455/1270) -67.4 E42-E34 N24-N36
      6 1199 1026 (N03W10 L=215)   714    232        1 (1176/ 926-1219/1147)  57.5 W06-W14 S02-N09
      7 1770  900 (S21W46 L=252)   355     51        0 (1755/ 864-1785/ 929) -90.0 W45-W48 S23-S19
      8 1804  851 (S25W50 L=255)   883    102        0 (1761/ 819-1850/ 887)   0.0 W44-W55 S27-S24
      9 1733  501 (S50W47 L=252)   819     67        1 (1691/ 483-1769/ 525) -14.0 W39-W56 S52-S46
     10 1499  462 (S42W18 L=223)  1138    181        2 (1420/ 401-1571/ 518)  12.1 W10-W24 S44-S38
     11 1561  372 (S53W25 L=230)   489     48        1 (1527/ 357-1592/ 391) -12.5 W20-W29 S55-S51
  

    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