Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  791 1486 (N35E04 L=175)   735    153        1 ( 728/1440- 868/1519)  -5.7 E08-W02 N33-N36
      2 1708 1394 (N08W54 L=233)   458     84        3 (1676/1376-1740/1414) -18.4 W51-W57 N06-N09
      3 1744 1325 (N03W55 L=234)   393    135        1 (1703/1303-1787/1347) -35.5 W50-W58 N00-N05
      4  199 1034 (N20E57 L=122)   377    138        1 ( 169/ 998- 240/1058) -39.8 E60-E53 N17-N22
      5 1073  861 (S07E01 L=178)  3271    898        9 ( 878/ 605-1199/1054)  72.3 E12-W06 S21-N04
      6  208  889 (N12E59 L=120)   696    102        2 ( 191/ 836- 224/ 945)  69.0 E61-E57 N08-N15
      7  160  903 (N13E65 L=114)   506    113        2 ( 130/ 868- 188/ 929) -45.0 E68-E61 N10-N15
      8  871  711 (S11E16 L=163)   190     56        0 ( 850/ 679- 894/ 747)  36.9 E17-E13 S13-S09
      9 1259  637 (S25W05 L=184)   516    146        1 (1191/ 592-1337/ 677)   5.7 W00-W10 S26-S23
     10  620  354 (S28E46 L=133)   935    240        3 ( 583/ 273- 651/ 447)  65.2 E50-E42 S34-S21
  

    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