Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  555 1709 (N41E55 L=122)   490     80        0 ( 520/1687- 593/1734)  -7.6 E63-E48 N39-N42
      2 1159 1597 (N40W03 L=180)   651    129        1 (1102/1564-1223/1637) -26.6 E01-W08 N38-N42
      3  512 1431 (N21E41 L=135)   421     88        1 ( 494/1397- 538/1474)  90.0 E42-E38 N18-N23
      4 1329 1233 (N17W17 L=194)   707    303        1 (1264/1179-1386/1316) -69.4 W12-W21 N14-N22
      5  195  952 (S13E63 L=114)  3930    309        4 ( 125/ 824- 279/1055)  29.7 E81-E52 S21-S05
      6 1704  699 (S12W53 L=229)   688    228        4 (1621/ 634-1778/ 777)  42.3 W47-W58 S16-S06
      7 1571  735 (S11W40 L=217)   370     71        1 (1555/ 698-1593/ 775)  90.0 W38-W41 S13-S08
      8 1489  702 (S14W35 L=211)   329     81        1 (1451/ 682-1535/ 720)  31.0 W32-W37 S15-S12
      9 1103  611 (S24W10 L=187)   348    155        0 (1047/ 566-1140/ 673)  58.0 W06-W11 S27-S19
     10  788  613 (S27E11 L=165)   217     12        0 ( 768/ 596- 808/ 630) -45.0 E12-E09 S28-S26
     11 1537  516 (S26W45 L=221)  1064    123        4 (1477/ 486-1595/ 547)  18.4 W39-W50 S28-S23
     12 1104  461 (S34W13 L=190)   453    129        0 (1044/ 444-1166/ 478)   6.3 W08-W17 S35-S32
  

    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