Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  838 1767 (N55E22 L=141)   547     82        2 ( 795/1737- 890/1788)  -5.2 E27-E16 N51-N56
      2 1019 1731 (N52E02 L=161)   708    140        0 ( 962/1706-1079/1750) -11.3 E07-W03 N49-N53
      3  434 1543 (N35E51 L=112)  2294    321        5 ( 324/1462- 557/1642) -11.3 E81-E36 N29-N40
      4 1533 1558 (N38W42 L=205)   698    104        0 (1481/1532-1579/1580) -12.5 W36-W46 N35-N39
      5 1594 1485 (N32W45 L=208)   328     48        0 (1574/1456-1611/1510) -90.0 W43-W46 N30-N34
      6  820 1444 (N29E15 L=148)  1242    383        2 ( 640/1409-1002/1475)  -4.4 E28-E02 N27-N31
      7  288 1380 (N23E59 L=104)   424     76        0 ( 259/1358- 323/1405)  -6.3 E64-E55 N21-N24
      8 1619 1351 (N23W43 L=206)   370     71        0 (1591/1311-1639/1381) -78.7 W40-W44 N20-N25
      9 1854 1224 (N15W65 L=228)  1294    184        3 (1799/1146-1913/1284) -36.9 W58-W73 N09-N18
     10 1257 1188 (N13W14 L=177)   317    111        0 (1210/1151-1305/1213)  26.6 W11-W17 N11-N15
     11  629 1159 (N11E25 L=137)   314    104        0 ( 588/1129- 670/1196) -26.6 E28-E22 N08-N12
     12 1768 1136 (N10W53 L=216)   458     89        3 (1739/1102-1796/1171) -53.1 W50-W56 N07-N11
     13 1440 1126 (N10W26 L=189)   217     71        0 (1410/1095-1470/1152) -45.0 W24-W28 N07-N11
     14  679  782 (S12E22 L=141)   210     49        2 ( 656/ 759- 704/ 804)   0.0 E23-E19 S13-S10
     15  983  586 (S25E02 L=161)   303    116        0 ( 962/ 541-1001/ 646)  81.9 E03-W00 S27-S20
     16 1266  547 (S27W17 L=180)   690    138        1 (1219/ 516-1325/ 587) -24.0 W13-W22 S29-S24
     17 1066  555 (S27W04 L=167)   527    121        0 ( 998/ 528-1118/ 585)  29.7 E01-W06 S28-S24
     18  733  222 (S57E32 L=131)   669     69        0 ( 699/ 204- 773/ 243) -23.2 E36-E27 S58-S54
  

    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