Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  928 1878 (N56E43 L=355)   551     73        0 ( 894/1868- 978/1892)  18.4 E49-E37 N53-N57
      2 1251 1783 (N56E03 L= 35)  1801    229        3 (1130/1746-1335/1813)  30.1 E13-W06 N50-N63
      3 1191 1389 (N24W03 L= 41)   285     96        0 (1153/1369-1229/1417) -11.3 W00-W05 N22-N24
      4 1777 1098 (N18W51 L= 89)  2714    317        6 (1698/ 997-1840/1227) -58.4 W43-W59 N11-N26
      5 1871 1049 (N18W63 L=101)  1835    228        6 (1836/ 945-1917/1120) -84.8 W58-W71 N10-N21
      6  334 1057 (S13E47 L=351)   444    118        0 ( 302/1012- 379/1098)  50.2 E49-E43 S15-S09
      7  932  759 (S19E00 L= 38)   534    176        0 ( 859/ 735- 988/ 794)  26.6 E04-W03 S20-S16
      8  196  737 (S36E66 L=332)   678     91        3 ( 178/ 708- 220/ 772) -38.7 E69-E61 S37-S33
      9 1396  560 (S21W35 L= 73)   755    214        2 (1326/ 511-1466/ 654)  70.0 W31-W37 S25-S14
     10 1137  471 (S33W21 L= 59)   284     49        0 (1120/ 445-1156/ 499) -36.9 W18-W23 S34-S31
     11 1527  378 (S30W57 L= 95)  2528    349        4 (1366/ 313-1657/ 443)  36.4 W46-W67 S37-S22
  

    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