Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  575 1598 (N15E45 L=136)   308     65        0 ( 542/1563- 602/1626)  80.5 E46-E43 N11-N17
      2  672 1399 (N05E30 L=151)   245     71        0 ( 634/1373- 709/1419)  45.0 E32-E28 N02-N06
      3 1253 1332 (N16W04 L=186)   502    211        0 (1190/1251-1292/1412) -41.2 E01-W07 N11-N19
      4  223 1240 (S13E61 L=121)   619    197        2 ( 189/1203- 268/1286) -18.4 E66-E54 S14-S10
      5 1079 1186 (N03E01 L=180)   232     61        1 (1052/1155-1099/1215) -33.7 E03-W00 N01-N04
      6  365  939 (S28E40 L=141)  1061    254        0 ( 275/ 873- 447/1005)  42.5 E47-E34 S33-S22
      7 1364  737 (S13W28 L=209)   204    109        0 (1326/ 696-1383/ 793) -33.7 W23-W30 S14-S10
      8 1514  769 (S06W36 L=217)   507    173        2 (1432/ 744-1583/ 791)  26.6 W30-W40 S09-S04
      9  883  580 (S36W05 L=186)   875    243        1 ( 761/ 511- 991/ 629)   0.0 E04-W14 S37-S32
     10 1183  562 (S28W24 L=205)   743    237        0 (1130/ 486-1267/ 620) -84.3 W20-W29 S31-S21
     11 1647  387 (S21W67 L=248)  2009    313        4 (1583/ 352-1734/ 425)  28.1 W59-W76 S23-S15
  

    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