Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  802 1773 (N34E43 L=171)   427     54        0 ( 786/1750- 822/1806) -45.0 E45-E40 N32-N36
      2  564 1675 (N21E53 L=161)   541     99        0 ( 547/1637- 594/1724) -40.6 E57-E48 N18-N24
      3  948 1485 (N19E18 L=196)   486    142        0 ( 897/1425- 983/1550) -29.1 E23-E14 N15-N21
      4 1471 1446 (N33W16 L=230)   569    166        0 (1433/1372-1492/1509) -53.1 W11-W18 N27-N36
      5 1299 1429 (N26W04 L=219)   435    142        0 (1249/1369-1324/1473) -35.5 W00-W07 N22-N27
      6  341 1344 (S04E54 L=161)  2972    486        5 ( 236/1238- 430/1466)  65.0 E67-E46 S11-N04
      7 1675 1389 (N37W34 L=248)   348     58        0 (1655/1356-1690/1415) -63.4 W31-W35 N34-N38
      8 1223 1359 (N19W02 L=216)   229     55        1 (1205/1331-1243/1380)  71.6 W00-W03 N17-N20
      9 1246 1280 (N15W05 L=220)   266     80        0 (1211/1251-1283/1300)  -9.5 W02-W08 N14-N16
     10 1742 1125 (N22W44 L=258)   409     86        0 (1719/1076-1759/1161) -76.0 W40-W45 N19-N24
     11  381 1000 (S23E41 L=173)   790    199        1 ( 308/ 939- 434/1075)  45.0 E47-E37 S27-S17
     12  795 1004 (S12E12 L=202)   211     51        1 ( 761/ 986- 823/1021)  26.6 E14-E10 S13-S10
     13  899  918 (S14E04 L=210)   272     40        0 ( 879/ 894- 920/ 946) -33.7 E05-E02 S15-S12
     14  985  901 (S13W01 L=215)   213     66        1 ( 957/ 877-1012/ 932) -33.7 W00-W03 S14-S11
     15 1069  869 (S12W07 L=221)   599    185        3 (1032/ 802-1104/ 930) -69.4 W04-W09 S16-S08
     16 1156  865 (S10W12 L=226)   386     98        2 (1106/ 837-1208/ 897)  38.7 W09-W14 S12-S07
     17  708  738 (S30E11 L=203)   733    207        1 ( 633/ 707- 805/ 791)   0.0 E17-E04 S31-S27
     18  605  451 (S53E13 L=202)   367     54        0 ( 576/ 427- 636/ 475)  -8.1 E15-E08 S54-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